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Synoptic  monitoring  of  Jupiter's  decameter-wavelength  emissions  over 
many  years  has  yielded  a  considerable  accumulation  of  data.  This 
dissertation  summarizes  an  analysis  of  the  data  for  effects  not  readily 
apparent  from  the  study  of  small  subsets  of  the  observations.  The  data 
base  is  a  merger  of  the  observations  collected  by  the  Universities  of 
Florida  and  Texas.  It  is  contained  on  magnetic  tape  in  punched  card 
format  and  it  has  been  processed  into  an  intermediate  form  more  suitable 
for  direct  analysis.  Information  necessary  for  use  of  the  raw  data  and 
intermediate  data  magnetic  tapes  has  been  provided,  as  well  as  details  of 
the  variety  of  punched  card  formats. 

The  data  in  each  of  ten  frequencies  in  the  range  from  5.6  to  30  MHz 
are  summarized  in  histograms  of  occurrence  probability  for  2°  intervals 
in  System  III  longitude  or  the  position  of  Io  from  superior  geocentric 
conjunction  (phase  of  Io).  Apparent  in  the  high-resolution  longitude 
histograms  are  spectral  changes  in  profile  and  position  of  the  three 
probability  peaks,  often  called  radio  sources  A,  B,  and  C.  Examination 
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of  the  data  comprising  each  source  for  dependence  of  occurrence 
probability  on  the  phase  of  Io  reveals  asymmetries  in  the  dominant  peaks, 
often  accompanied  by  proximate  secondary  maxima.  High-resolution 
occurrence  probability  versus  longitude  histograms  for  each  apparition 
of  the  merged  18,  20,  and  22-MHz  observations  confirm  the  sinusoidal 
drift  of  source  A  with  time,  and  further  studies  emphasize  the  respon- 
sibility of  the  non-Io-related  radiation  for  this  effect.  No  significant 
trend  of  changes  in  the  longitude  positions  of  sources  B  and  C  has  been 
found,  but  the  Io-related  parts  of  these  sources  do  exhibit  a  twelve-year 
cyclic  drift  in  Io  phase,  while  source  A  shows  no  consistent  Io-phase 
motion. 

Histograms  of  occurrence  probability  calculated  for  intervals  in 
both  longitude  and  Io  phase  provide  a  much  deeper  insight  into  the  source 
structure.  Data  at  eight  frequencies  from  10  to  30  MHz  have  been  reduced 
in  this  manner,  four  with  interval  sizes  of  2°  by  2°  and  the  other  four 
with  5°  resolution  in  each  variable.  The  diagrams  reveal  many  frequency- 
dependent  changes  in  the  sources,  such  as  the  varying  appearance  of  two 
projections  on  the  lower  Io-phase  side  of  Io-related  source  B,  a  bridge  of 
radiation  connecting  sources  A  and  C  through  the  source  B  longitudes  at 
lower  frequencies,  and  a  large  modulation  of  the  source  positions  and  shapes 
between  15  and  10  MHz.  Similar  histograms  were  constructed  for  the 
isolated  Io-related  source  regions  from  combined  18,  20,  and  22-MHz 
data  for  individual  apparitions.  The  connection  of  the  motion  of 
sources  B  and  C  in  Io  phase  with  changes  in  their  morphology  is  evident 
from  the  results. 

The  phenomena  uncovered  in  the  histogram  analyses  can  be  compared 
with  theories  of  beaming  of  the  radiation  from  the  emission  region.  Few 


theories  treat  this  subject  in  numerical  detail,  but  correspondence  of 
some  of  the  source  structure  with  predictions  by  McCulloch,  Goertz, 
and  Cavaliere  and  Speranza  is  found.  Calculations  of  the  expected 
Io-phase  and  System  III  longitude  angles  for  which  Io-related  source 
radiation  will  be  received  on  earth  have  been  made  using  the  magnetic 
field  models  of  Acuna  and  Ness  and  Smith  et  al .  in  combination  with  the 
beaming  geometry  of  the  Goldreich  and  Lynden-Bell  theory.  In  general, 
agreement  with  the  observed  source  characteristics  is  good,  but  further 
calculations  indicate  that  more  of  the  high-resolution  phenomena  can  be 
explained  by  occultation  of  the  emission  regions  by  the  dense  Jovian 
atmosphere  near  the  top  of  the  cloud  layer. 


vm 


CHAPTER  I 
INTRODUCTION 

1.1  The  Decametric  Jupiter  Puzzle  in  Plain  Language 

Over  20  years  ago  Burke  and  Franklin  (1955)  first  realized  that  the 
source  of  interference  falling  on  their  new  radio  telescope  was  the 
planet  Jupiter.  Since  that  time  millions  have  been  spent  on  observation 
sites  around  and  above  the  world  to  record  every   facet  of  the  radio 
noise  from  Jupiter. 

Should  the  emissions  be  called  noise?  Noise  implies  a  hodgepodge 
of  data,  a  random  mixture  of  the  electromagnetic  waves.  Although  Jupiter 
is  certainly  not  perfectly  predictable,  it  exhibits  regularities  which 
make  it  one  of  the  most  interesting  paradoxes  in  the  solar  system. 

All  planets  emit  radio  noise  as  a  natural  consequence  of  the  thermal 
motion  of  their  molecules.  Jupiter,  the  sun,  the  earth,  and  possibly 
Saturn  show  emission  far  exceeding  that  explainable  by  temperature  alone. 
On  the  earth  we  know  that  part  of  this  is  due  to  our  communications  net- 
work and  the  rest  to  lightning  storms  and  the  interaction  of  electrons 
and  ions  with  the  earth's  magnetic  field.  Extrapolating  this  to  Jupiter 
implies  the  existence  of  Jovian  magnetic  fields,  lightning  storms, 
and/or  intelligent  life.  Could  the  emissions  be  untranslated  communi- 
cation? Since  planetary  conditions  are  inhospitable  to  life  as  we  know 
it,  theories  based  on  natural  origins  prevail.  Some  of  the  characteristics 
of  the  radiation,  however,  seem  very   unnatural. 


Consider  the  following  facts. 

1)  The  excessive  radio  noise  from  Jupiter  which  is  detected  on 
earth  occurs  mostly  in  the  wavelengths  around  7-60  meters,  referred  to 
as  decameter  wavelengths.  The  corresponding  radio  frequencies  are 
5-40  megahertz  (MHz),  which  are  in  the  shortwave  radio  band.  It  is 
interesting  that  these  are  the  bands  which  are  used  for  our  long 
distance  communications.  They  are  very   prone  to  static  interference  and, 
sometimes,  the  swishing  sound  breaking  up  communications  originates  from 
Jupiter. 

Jupiter's  decametric  radiation  is  not  heard  continuously.  The 
probability  of  receiving  it  depends  on  which  side  of  the  planet  is  toward 
the  earth.  Unfortunately,  the  "side"  of  the  planet  which  faces  the  earth 
is  not  as  easily  defined  as  people  living  on  a  solid  surface  would  expect. 
Gaseous  orbs,  like  Jupiter  and  the  sun,  rotate  at  different  rates  in 
different  latitudes.  In  the  case  of  Jupiter,  a  complete  rotation  takes 
approximately  9  hours  and  50.5  minutes  for  the  equatorial  cloud  band,  and 
9  hours  and  55.5  minutes  for  the  clouds  away  from  the  equator  (see 
Table  1-1).  The  major  amount  of  radio  emission  occurs  regularly  at  a 
period  close  to  9  hours  55.5  minutes,  so  that  at  first,  the  non-equatorial 
regions  seem  to  be  a  likely  point  of  origin. 

The  radio  period  is  a  small  but  significant  amount  different  from 
that  of  the  cloud  bands,  however,  and  very  regular.  Since  the  period 
does  not  match  that  of  the  visible  features,  this  regularity  suggests 
equating  the  period  to  that  of  a  rotating,  solid  object,  perhaps  the 
core  of  Jupiter  if  one  exists.  Carrying  this  further  one  could  speculate 
that  the  radio  signals  may  be  related  to  a  particular  surface  area  or 
"side"  of  the  core  of  Jupiter.  This  is  often  called  the  source  region 
or  just  a  source. 


Table  1-1 
Important  Jovian  Parameters 


(Rj) 
fe) 


Equatorial  Radius 

Ellipticity 

Semimajor  Axis  of  Orbit 

Eccentricity  of  Orbit 

Sidereal  Period 

Synodic  Period 

Inclination  of  Orbit  to  Ecliptic 

Inclination  of  Equator  to  Orbit 

Sidereal  Rotation  Periods 
System  I    Epoch  1897.0  (At  ) 
System  II   Epoch  1897.0  (An  ) 
System  III  Epoch  1957.0  (Attt  ) 
System  III  Epoch  1965  proposed 


71372  km 
0.06466 
5.202803  A.U. 
0.04845 
11.86223  yrs 
398.88  days 
1°  18*  17" 
3°  5' 

9 J  50m  30s,  003 

9"  55m  40s.  632 

9^  55m  29s.  37 

9n  55m  29s.  711 


(1) 

(2) 
(3) 
(3) 
(3) 
(3) 
(3) 
(3) 

(4) 
(4) 
(5) 
(6) 


Satellite  Data 


Satellite 

Sidereal 

Orbit 

No.    Name 

Period 

Incl. 

Ecc. 

Mean  Dist. 

Radius 

(days) 

(°) 

(Rj) 

(km) 

(4) 

1.76986049 

(3) 
0 

(3) 
0.000 

(3) 

(3) 

I    lo 

5.91 

1810 

II   Europa 

3.55409417 

1 

0.000 

9.40 

1480 

III  Ganymede 

7.16638722 

0 

0.001 

14.99 

2600 

IV   Callisto 

16.75355227 

0 

0.007 

26.38 

2360 

V    Amalthea 

0.49823633 

0 

0.003 

2.54 

80 

References: 

(1)  Opp,  1975 

(2)  Anderson  et  al . ,  1974 

(3)  Allen,  1973 

(4)  American  Ephemeris  and  Nautical  Almanac,  1976 

(5)  I.A.U.  Information  Bulletin,  1962 

(6)  Riddle  and  Warwick,  1976 


2)  Receiving  radiation  from  certain  sides  of  Jupiter  is  unusual  but 
not  necessarily  unnatural.  For  example,  the  earth  has  regions  which  are 
more  prone  to  lightning  storms.  This  type  of  radiation  occurs  somewhat 
periodically,  but  it  is  not  nearly  as  regular  as  the  radiation  from 
particular  source  regions  of  Jupiter.  These  are  considerably  enhanced 
at  very  regular  intervals.  Bigg  (1964)  discovered  that  the  intervals 
were  actually  the  times  when  Io,  the  first  Galilean  moon,  was  near  one 
of  two  positions  in  its  orbit  as  observed  from  the  earth.  If  this  is 
not  strange  enough,  the  orbital  positions  at  which  the  enhancements  take 
place  are  not  symmetrical.  If  we  consider  the  time  Io  is  directly 
behind  Jupiter,  as  viewed  from  the  earth,  as  a  starting  point  (i.e., 
superior  geocentric  conjunction)  then  one  of  the  enhancement  positions 
is  one-fourth  of  the  way  around  Io's  orbit  and  the  other  is  two-thirds 
of  the  way.  Dependence  of  radio  events  on  the  orbital  phase  of  a  moon 
sounds  not  only  unnatural,  but  astrological. 

3)  According  to  Desch  and  Carr  (1974)  and  Brown  (1974a),  Jupiter 
is  most  active  at  a  frequency  of  approximately  8  MHz.  The  earth's 
ionosphere  interferes  with  observations  at  frequencies  below  this  value. 
At  higher  frequencies,  Jupiter's  activity  drops  off  at  a  steady  rate- 
until  a  value  of  39.5  MHz  is  reached.  Beyond  this  point  no  abnormal 
decametric  radiation  has  ever  been  received. 

4)  Although  not  as  severe  as  the  earth's  23.5°  angle,  Jupiter's 
rotation  axis  has  a  3°  tilt  so  that  as  Jupiter  orbits  the  sun  the 
northern  and  southern  hemispheres  are  tipped  toward  the  earth  in  12-year 
cycles.  When  the  northern  hemisphere  is  dominant,  radiation  is  received 
more  often,  and  the  area  of  the  sources  appears  to  be  larger.  It  is  as 


though  the  radiation  were  being  beamed  like  a  searchlight  from  the 
northern  hemisphere  and,  when  this  hemisphere  is  tipped  toward  the 
earth,  the  searchlight  is  aimed  more  directly  at  us.  Also,  the  positions 
of  the  radio  sources  shift  slightly  in  longitude  depending  upon  how 
much  the  northern  hemisphere  is  tipped  toward  or  away  from  the  earth. 

5)  The  natural  thermal  radio  noise  emitted  by  a  planet  is  composed 
of  electromagnetic  waves  having  random  vibrational  directions  of  the 
electric  and  magnetic  components.  The  decametric  radiation  from  Jupiter, 
however,  is  polarized  so  that  the  vibrations  of  the  electric  and 
magnetic  waves  appear  to  move  in  an  elliptical  path  as  the  electromagnetic 
wave  propagates  toward  the  observer.  For  most  of  the  radiation  sources 
the  vibrations  appear  to  be  right  elliptically  polarized;  that  is,  they 
appear  to  the  observer  to  travel  counterclockwise  around  the  elliptical 
path.  For  the  other  source  regions  the  predominant  mode  is  in  the  opposite 
sense  or  left  elliptically  polarized. 

6)  The  radio  signals  are  not  confined  to  a  single  frequency.  On 
a  loudspeaker  connected  to  a  single  frequency  receiver  Jupiter  radio 
storms  sound  like  ocean  waves  breaking  on  a  shore.  These  "swishing 
storms"  drift  slowly  upward  or  downward  in  frequency  in  a  manner  which 
depends  upon  the  longitude  of  Jupiter  and  often  the  orbital  position  of 
Io  as  well.  At  times,  events  are  received  which  sound  like  "pops"  on  a 
speaker.  When  slowed  down  about  a  hundred  times  these  pops  sound  more 
like  whistles  of  rapidly  decreasing  pitch,  usually  one  right  after 
another.  These  are  detected  only  when  Io  and  Jupiter  are  in  particular 
orientations  to  each  other. 


Combining  all  of  the  six  sections  together  one  could  argue  for 
Jovian-manned,  frequency-modulated  radio  stations  in  the  northern 
hemisphere  of  Jupiter  which  transmit  narrow  beams  of  polarized  radiation 
when  Io  passes  over  them.  Such  a  situation  is  possible,  but  it  is  part 
of  man's  nature  to  look  for  the  simplest  explanation  first.  Indeed  many 
plausible  theories  have  explained  the  above  phenomena  in  terms  of  the 
natural  ingredients  of  charged  particles  and  magnetic  fields,  although 
so  far  no  one  theory  has  satisfactorily  explained  everything.  The  lack 
of  dominance  for  any  particular  theory  is  due  partially  to  the  indefinite 
characteristics  of  the  decametric  sources.  This  allows  a  wide  range  of 
source  mechanisms  to  yield  the  same  general  results.  The  aim  of  the 
research  presented  in  this  dissertation  is  to  study  the  source  phenomena 
with  higher  resolution  and  so  put  greater  constraints  on  the  theoretical 
models.  The  results  of  this  endeavor  will  be  discussed  in  the  following 
chapters,  but  first,  there  is  a  need  to  step  back  and  see  the  overall 
picture  from  a  more  scientific  viewpoint. 

1 . 2  Decametric  Jupiter:  A  Brief  Technical  Summary 

In  order  to  familiarize  the  reader  with  the  terms  and  generally 
accepted  hypotheses  in  Jovian  decametric  study  the  characteristics  outlined 
in  the  previous  section  will  be  presented  again  in  more  technical 
language.' 

1.2.1  Longitude  Dependence 

The  reception  of  decametric  radiation  (often  referred  to  as  DAM) 
depends  on  which  Jovian  longitude  is  on  the  central  meridian  as  viewed 
from  the  earth.  It  should  be  mentioned  that  measurement  of  longitude 
on  Jupiter  is  based  on  the  definition  of  0°  longitude  coinciding  with  the 


central  meridian  as  viewed  from  earth  at  a  particular  point  in  time  or 
epoch.  Longitude  measurement  differs  from  the  earth's  longitude  system 
in  that  longitude  increases  westward  from  the  0°  meridian  up  to  a  maximum 
of  360°.  The  rotational  dependence  of  activity  is  normally  presented 
in  the  form  of  a  histogram,  as  in  Figure  1-1.  This  figure  shows  the 
probability  of  receiving  DAM  at  22  MHz  for  each  5°  interval  of  longitude 
as  averaged  over  the  years  1962-1967.  The  occurrence  probability  (o.p.) 
is  usually  defined  as  the  ratio  of  the  number  of  activity  events  to  the 
number  of  listening  times  for  each  histogram  interval.  Although  the 
histograms,  change  slightly  in  appearance  with  frequency  and  time, 
Figure  1-1  represents  a  typical  plot  for  the  middle  decametric  frequency 
range.  Labeled  A,  B,  and  C  in  the  order  of  their  prominence,  the  three 
probability  peaks  are  referred  to  as  sources,  although  the  name  implies 
more  individuality  than  it  should.  In  some  of  the  literature  A,  B,  and 
C  are  called  the  main,  early,  and  late  sources.  At  lower  frequencies  a 
fourth  source  sometimes  appears  in  the  20°-100°  longitude  range,  called 
appropriately  D,  or  the  fourth  source. 

As  mentioned  before  there  is  a  period,  different  from  the  cloud 
layer  rotation  periods,  which  best  represents  the  rotation  of  the  radio 
sources.  The  longitude  system  based  on  this  period  is  called  System  III 
longitude  or,  symbolically  (Xm).  Systems  I  and  II  are  the  optically 
observed  periods  of  the  equatorial  and  non-equatorial  region  cloud  bands, 
respectively  (see  Table  1-1).  The  exact  period  of  \         has  been  the 
subject  of  much  discussion  (Carr,  1971;  Alexander,  1975;  Duncan,  1971; 
Lecacheux,  1974).  At  a  recent  meeting  in  Tucson,  many  of  those  concerned 
agreed  upon  a  consensus  value  of  9h  55m  29?711  to  be  submitted  to  the 
International  Astronomical  Union  for  adoption  as  the  official  System  III 
longitude  epoch  1965  to  replace  the  epoch  1957  value  of  9h  55m  29*37. 
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Since  Am  differs  from  Aj  and  X       the  source  of  the  radio  signals 
should  not  be  related  to  the  visible  cloud  surface.  Radio  observations 
in  the  decimetric  (DIM)  frequency  range  (Drake  and  Hvatum,  1959; 
Roberts  and  Ekers,  1966;  Berge,  1966;  Branson,  1968)  indicated  the 
presence  of  radiating  charged  particles.  Changes  in  the  polarization 
and  intensity  of  the  emission  imply  the  presence  of  a  rotating  dipolar 
magnetic  field.  The  period  of  this  rotation  corresponds  to  the  Am 
period  (Carr,  1971).  The  particles  are  apparently  trapped  by  and 
radiate  in  a  Jovian  dipole-like  field,  which  would  mean  the  System  III 
period  represents  the  rotation  of  the  magnetic  field.  Pioneer  10  and  11 
satellite  data  (Smith  et  al . ,  1975;  Van  Allen  et  al . ,  1975)  have  confirmed 
that  periodic  increases  in  energetic  particle  intensity  were  in  synchronism 
with  the  time  of  rotation  for  Jupiter's  magnetic  field.  The  magnetic  field 
of  Jupiter,  like  that  of  the  earth,  is  believed  to  arise  from  electric 
currents  circulating  deep  in  the  planet.  On  Jupiter  the  pressure  is 
sufficient  that  the  currents  may  be  electric  charges  attached  to  a 
rotating  solid  core  of  metallic  hydrogen.  Thus,  in  their  relation  to 
the  magnetic  field,  the  radio  signals  may  yield  information  about  the 
hypothetical  solid  surface  of  Jupiter. 

1.2.2  Io  Dependence 

By  far  the  most  unusual  characteristic  of  DAM  is  its  dependence  on 
the  position  of  the  first  Galilean  satellite,  Io,  as  discovered  by 
Bigg  (1964).  This  dependence  is  best  described  by  adding  a  new  dimension 
to  the  histogram  of  Figure  1-1.  Figure  1-2  is  a  smoothed  contour  plot 
of  1962-1967,  22-MHz  data  computed  from  values  of  occurrence  probability 
for  5°  intervals  in  System  III  longitude  and  the  phase  of  Io  from 
superior  geocentric  conjunction  (<K  ). 
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Figure  1-2.  Smoothed  contour  plot  of  the  number  of  radio  storms  in 
each  5°  interval  of  System  III  Longitude  (65)  and  the 
phase  of  Io  from  superior  geocentric  conjunction.  The 
graph  is  based  on  Florida  22-MHz  data  for  the  years 
1962-1967.  The  contour  interval  is  6. 
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The  features  in  Figure  1-2  are  typical  for  the  20-MHz  frequency 
range.  The  eye  is  naturally  drawn  to  source  B  (60°  <Am  <200°),  which 
is  conspicuous  in  its  Io  dependence.  Only  when  Io  is  near  90°  from 
superior  geocentric  conjunction,  or  90°  west  of  the  central  meridian 
as  seen  from  the  earth,  can  source  B  longitude  region  emission  normally 
be  detected.  The  source  A  region  (200°  <Ani  <290°)  is  quite  different. 
Radiation  is  often  received  regardless  of  the  position  of  Io,  as  can  be 
seen  from  the  vertical  smear  in  the  A  longitudes.  This  is  usually 
labeled  the  non-Io-A  radiation.  There  is  an  enhancement  in  the  region 
180°  <  <J>Io  <  270°  where  the  occurrence  probability  is  a  factor  of  4  to 
5  larger  than  the  non-Io-A.  This  Io-A  region  has  the  same  Io-phase 
dependence  as  the  Io-C  enhancement  just  to  the  right  of  it.  In  the 
C  longitudes  (290°  <  Am  <  360°,  0°<A  n<40°)  there  is  some  non-Io- 
related  emission,  although  it  is  not  nearly  as  prominent  as  non-Io-A. 
Thus,  source  C  is  more  like  B  in  its  isolation.  Source  D  is  only  seen  at 
frequencies  below  18  MHz  in  the  range  of  20°-100°  of  longitude  and  it  is 
isolated  in  the  95°-115°  region.  The  other  sources  also  change  with 
frequency,  as  will  be  seen  in  the  next  section. 

1.2.3  Source  Frequency  and  Sensitivity  Dependence 

The  two-dimensional  histogram  changes  gradually  with  frequency.  When 
frequency  is  increased  the  Io-C  source  slowly  shrinks  until  it  disappears 
around  26  MHz.  The  non-Io-A  vanishes  at  28  MHz,  and  the  Io-A  disappears 
at  approximately  36  MHz.  Finally,  no  decametric  storms  have  ever  been 
found  to  exceed  39.5  MHz  (Warwick,  1970).  A  decrease  in  frequency  brings 
a  decrease  in  the  Io-related  radiation  and  a  growth  in  non-Io-related 
radiation  such  that  a  histogram  like  that  of  Figure  1-2  would  have  contours 
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in  all  regions.  Data  below  10  MHz  are  difficult  to  acquire  because 

of  the  increasing  opacity  of  the  earth's  ionosphere  to  these  frequencies. 

Satellite  studies  are  beginning  to  fill  in  this  gap  (Brown,  1974b; 

Desch  and  Carr,  1974).  There  is  still  insufficient  data  to  say  much 

about  the  dependence  of  emission  on  longitude  and  moon  phase  at  these 

lower  frequencies  although  some  interesting  trends  have  been  found. 

For  further  details  on  the  satellite  dependence  of  the  sources,  the  summary 

by  Carr  and  Desch  (1975)  is  excellent. 

Bozyan  and  Douglas  (1975)  report  that  they  have  discovered  a 
separate  peak  of  the  non-Io-related  radiation  in  the  C  region.  According 
to  them,  there  is  a  distinct  non-Io-C  source  rather  than  a  spillover  of 
the  non-Io-A  radiation  into  the  C  longitudes.  If  sources  A  and  C  each 
have  a  non-Io-related  component  then  maybe  the  other  sources  do  also. 

In  support  of  this  statement,  Desch  and  Carr  (1974)  have  found  a 
spread  of  non-Io-B  related  radiation  using  a  640-element  dipole  array  at 
26.3  MHz.  This  non-Io-B  evidence  suggests  that  each  source  does  indeed 
have  a  non-Io-related  component  which  can  be  detected  with  sufficiently 
sensitive  apparatus  when  source  B  longitudes  face  the  earth. 

1.2.4  Twelve-Year  Periodicities 

The  average  amount  of  Jupiter's  DAM  radiation  received  on  the  earth 
oscillates  with  a  period  between  11  and  12  years.  At  first,  this 
fluctuation  seemed  to  correlate  with  the  eleven  year  solar  activity 
cycle  (Smith  et  a! . ,  1965),  which  was  reasonable  since  the  sun  is  a 
likely  candidate  as  an  energy  source  for  DAM.  Douglas  (1964)  and 
Carr  (1962)  had  already  noticed,  however,  that  the  rhythm  also 
approximated  Jupiter's  sidereal  period  of  11.89  years.  Carr  suggested 
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the  periodicity  could  be  explained  by  the  changing  aspect  of  Jupiter  as 
seen  from  the  earth,  caused  by  the  3°  angle  between  Jupiter's  equatorial 
and  orbital  planes.  From  the  point  of  view  of  Jupiter  this  angle  is  the 
Jovicentric  declination  of  the  earth  as  symbolized  by  Df.     Maximum  or 
minimum  activity  occurs  when  the  Jovian  north  pole  is  pointed  toward  or 
away  from  the  earth  respectively.  The  increased  activity  can  be  noticed 
especially  in  source  A.  During  the  active  years,  the  width  of  the  source  A 
occurrence  probability  peak  in  A...  coordinates  increases,  as  does  the 
overall  reception  probability.  According  to  Carr  (1972)  this  indicates  a 
beaming  of  the  DAM  from  the  northern  hemisphere  such  that  when  this 
hemisphere  is  canted  toward  the  earth  a  greater  portion  of  the  radiation 
beam  is  received.  The  previously-mentioned  suggestion  of  a  pencil  beam 
radio  station  in  the  northern  hemisphere  is  not  too  far-fetched,  although  in 
this  case  the  station  appears  to  have  naturally-occurring  cyclotron  or 
synchroton  radiation  as  its  source.  This  is  inferred  primarily  from  the 
polarization  of  the  DAM. 

1.2.5  Polarization 

The  predominant  sense  of  polarization  in  the  source  A  and  B  longitude 
regions  is  right  elliptical  or  right  circular  (Kennedy,  1969).  This  is 
consistent  with  the  theory  of  cyclotron  or  synchroton  radiation  beamed 
from  the- northern  hemisphere  in  the  extraordinary  mode  (Dulk,  1965). 
The  magnetic  field  of  Jupiter  has  the  opposite  sense  to  that  of  the  earth, 
and  the  field  in  the  northern  hemisphere  will  be  pointing  out  of  the  planet 
and  toward  the  earth-bound  viewer.  In  the  same  manner,  sources  C  and  D, 
which  are  predominantly  left  elliptically  polarized,  are  speculated  to 
originate  in  the  southern  hemisphere.  The  polarization  argument  does 
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not  preclude  the  sources  from  being  in  either  hemisphere  or  even  far 
away  from  Jupiter's  surface,  but  this  evidence,  together  with  the 
magnetic  field  values  to  be  discussed  later,  make  a  convincing  argument 
for  assuming  sources  A  and  B  come  from  the  northern  half  of  Jupiter. 

1-2.6  Dynamic  Spectra 

The  simplest  approach  to  the  study  of  DAM  makes  use  of  narrow-band 
receiving  equipment.  A  single  Jupiter  storm  usually  covers  a  wide 
range  of  frequencies;  however,  and  a  wealth  of  information  can  be  gained 
from  broad-band  antennas  and  receivers.  This  is  especially  true  when 
the  storms  are  examined  with  high  time  resolution.  These  analyses 
generally  involve  very   small  subsets  of  the  data  base.  Since  this 
dissertation  is  concerned  mainly  with  large  volumes  of  narrow-band  data, 
only  a  brief  mention  of  the  field  of  dynamic  spectra  study  will  be  made. 

Jupiter  radiation  can  be  divided  into  the  categories  of  long  and 
short  events  or  L-bursts  and  S-bursts.  The  L-bursts  are  events  of  rela- 
tively long  duration  compared  to  S-bursts,  usually  on  the  order  of  several 
seconds.  These  were  mentioned  previously  as  "swishing"  storms.  The 
swishing  sound  comes  from  the  reception  of  a  succession  of  L-burst  waves 
by  narrow  bandwidth  equipment.  The  appearance  of  modulation  lanes  in 
adjacent  waves  or  envelopes  indicates  that  they  were  once  part  of  a 
continuous  disturbance  (Riihimaa,  1974).  The  separation  into  envelopes 
occurs  in  interplanetary  space  where  solar  plasma  winds  perturb  the  DAM 
waves  (Douglas  and  Smith,  1967).  Although  the  L-bursts  drift  up  and  down 
in  frequency  at  a  rate  of  a  few  MHz/sec, the  modulation  lanes  within 
L-bursts  have  drift  rates  around  100  kHz/sec.  The  rate  and  direction  of 
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frequency  drift  correlate  with  the  longitude  of  the  central  meridian 
and  the  position  of  Io  (Warwick,  1967). 

On  a  loudspeaker  connected  to  a  narrow-band  receiver,  S-bursts 
are  millisecond-length  events  often  described  as  spitting,  although 
percolating  might  be  more  acoustically  accurate.  The  tendency  for 
S-bursts  to  come  in  groups  results  in  a  repetitive  popping  noise.  With 
high  time  resolution  the  pop  is  seen  to  result  from  the  combination  of 
narrow  instantaneous  bandwidth  and  rapid  frequency  drift.  The  drift  is 
always  negative  and  porportional  to  frequency.  The  S-bursts  have  been 
observed  only  during  Io-controlled  storms,  especially  in  the  B  and  C 
longitude  regions.  For  more  detailed  information  on  the  dynamic  spectra 
phenomena  the  summary  by  Carr  and  Desch  (1975)  is  a  good  starting  point. 

Jupiter  radio  astronomy  basics  have  now  been  presented  both 
speculatively  and  technically.  Chapters  III  and  IV  will  expand  the 
previous  discussion  of  Sections  1.2.1-1.2.4,  while  the  material  in 
Sections  1.2.5  and  1.2.6  must  be  considered  in  the  theory  discussion  of 
Chapter  V.  For  the  moment  let  us  digress  from  the  fundamental  facts  to 
the  basis  for  these  facts,  namely  the  data. 


CHAPTER  II 
THE  DATA  BASE 

2.1  Introduction 

Despite  the  accumulation  of  many  astute  theories,  none  of  them  can 
as  yet  satisfactorily  explain  all  of  the  quirks  inherent  in  the  Jovian 
radiation  process.  Part  of  the  problem  in  establishing  a  comprehensive 
theory  of  Jovian  radiation  is  the  intermittent  behavior  of  the  radio 
events.  Their  nature  changes  not  only  with  the  rotation  of  Jupiter  and 
the  phase  of  Io,  but  also  with  the  11.86  year  orbital  period  of  Jupiter. 
Therefore,  two  decades  of  radio  monitoring  have  yielded  only  about  two 
cycles  of  data.  This  is  complicated  by  the  fact  that  ground-based  deca- 
metric  radio  observations  are  favorable  only  during  ionospherically 
quiescent  times,  usually  the  late  nighttime  hours.  The  irregularity  of 
good  conditions  for  reception  together  with  the  periodicities  make  it 
difficult  for  any   single  observing  site  to  get  an  overall  picture  of 
Jupiter's  activity. 

The  obvious  step  of  combining  the  data  from  several  observatories 
fills  in  .some  of  the  gaps  and  gives  a  better  statistical  representation 
of  the  emission  pattern.  Several  institutions  such  as  NASA  and  the 
Universities  of  Texas  and  Florida  have  helped  to  establish  multiple  ob- 
serving sites  for  this  purpose.  Some  individuals  have  accomplished  the 
same  result  by  merging  data  obtained  from  separate  groups.  Duncan 
(1971)  gathered  one  of  the  most  comprehensive  collections  of  the  earliest 
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years  of  Jupiter  decametric  observation.  Although  his  list  of  observing 
groups  is  extensive,  the  amount  of  data  is  small  since  most  of  the  groups 
listed  observed  only  in  the  mid-1950's.  Intensive  synoptic  monitoring 
of  Jupiter  did  not  begin  until  the  late  fifties  and  early  sixties,  during 
the  rising  segment  of  the  Jupiter  activity  cycle. 

The  amount  of  decametric  data  swelled  geometrically  in  the  mid-sixties 
because  of  increased  government  support  as  well  as  the  rise  in  Jovian 
emission  probability.  Exhaustive  collections  of  observations  have, 
therefore,  become  difficult  to  establish.  Duncan  includes  only  small 
subsets  of  the  data  from  groups  observing  in  the  sixties,  since  he  was 
interested  in  studying  only  the  higher  frequencies.  In  most  cases  recent 
publications  have  used  the  data  acquired  by  a  single  group. 

This  dissertation  is  based  upon  a  combination  of  surveillance  gathered 
by  groups  at  the  Universities  of  Florida  and  Texas.  The  University  of 
Florida  has  collected  data  from  sites  in  Maipu  and  Huanta,  Chile,  as  well 
as  from  its  own  observatories  in  Florida.  In  the  interest  of  brevity, 
this  assemblage  will  often  be  referred  to  as  Florida  data.  Likewise, 
the  term  Texas  data  will  mean  the  catalogued  collection  (Bozyan,  et  al . , 
1972)  of  data  from  observatories  at  the  University  of  Texas,  Yale,  and 
Kodaikanal,  India.  Both  Florida  and  Texas  began  synoptically  monitoring 
Jupiter  around  1957.  Both  have  continued  to  accumulate  observations  up 
to  the  present  time  from  many  widely-spaced  sites,  so  that  the  collection 
from  just  these  two  institutions  involves  over  145,000  hours  of  data. 
Figures  II-l  and  1 1  - 2  are  plots  of  the  listening  and  storm  times  through 
all  the  years  and  for  all  frequencies.  The  four-character  codes  along 
the  ordinate  axis  are  the  data  set  names.  These  names  are  listed  in 
Table  II-8  and  are  derived  from  the  card  codes  explained  in  Table  II-l. 


Figure  II— 1  -  Plot  of  listening  and  storm  times  for  all  Florida 

and  Chile  5°  intermediate  data  cards  on  magnetic  tape. 
The  codes  on  the  ordinate  axis  are  interpreted  in 
Table  II— 1  -  The  vertical  lines  are  placed  along  the 
Julian  data  axis  at  the  times  of  Jupiter's  conjunction 
with  the  sun.  See  Figure  1 1-3  for  an  explanation  of 
the  horizontal  bands. 


19 


i 

«\\ 

ii« 

nil 

ii 

N 

g 

'ii 

Iifl 

(Mi 

■*  :: 

iifl 

■j 
j 

i 

Iifl 

ii  ii 

iifl 

II 

H 

iiii 

II 

ll 

111 

ii 

n; 

r 

"  ii 

HI 

HI 

H 

II 

M 

n 

(\ 

-Hi 

i 

II 

i 

"II 

i 

1 

i 

\\ 

1 

Hi 

* 
i 

1! 

II 

Ml 

4 

ii 

1 

1 

i 

I 

-11 

jl 

1! 

ill 

Ml 

mi 

II 

4 

J 

si 

i 

11 

II 

ill 

ill 

M 

■1 

* 

; 

«(i 

II II 

ii 

ii 

•  ii 

::ll 

ill 

III 

*m 

II 

p«! 

1 

•    ii 
■•  ii 

.11 

111 

till 

nil 

;;l 

nil 

iiii 

Hi 

III 

ill 

lid 

ii; 

i 

J 

ill 

iifl 

III 

nil 

iilj 

nil 
»ii 

4| 

ii1 

"V 

H! 

< 

i 

Hi 

ii 

iill 

ill 

iii 

ii 

ill 

ill 

til 

•Ml 

! 

jiJi 

ai[ 

iifl 

nil 

ill 

*! 

ill 

Jl 

ii: 

■ 

"II 

"II 

ii 

=  n 

:i 

u 

;;ii 

=1 

ij 

i! 

I: 

1* 

J 

II 

!1 

ill 

n 

1 

"" 

■I 

Ih 

s 

4 

- 

ill 

II 

> 

Ml 

«l 

! 
"5 

5 

-rl 

~ 

::ii 

1 
■ 

1 

i  1 

i 

f 

3 

T 

= 

1 

? 

5 

1 

1 

S 

s 

s 

| 

5 

Figure  1 1-2.  Plot  of  listening  and  storm  times  for  all  Texas 

intermediate  data  on  magnetic  tape.  The  codes  on  the 
ordinate  axis  are  interpreted  in  Table  II-l.  The 
vertical  lines  are   placed  along  the  Julian  date  axis 
at  the  times  of  Jupiter's  conjunction  with  the  sun. 
See  Figure  1 1-3  for  an  explanation  of  the  horizontal 
axis. 
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The  data  are  plotted  according  to  Julian  date  and  separated  into 
apparitions  by  vertical  lines  representing  Jupiter's  times  of  con- 
junction with  the  sun. 

Figure  1 1 -3  shows  a  blowup  from  Figure  II-l  of  one  apparition  for 
one  frequency.  The  fine  detail  reveals  that  the  graph  is  made  of 
electrostatically  printed  dots  of  0.01  inch  diameter  produced  by  a  Gould 
plotter.  The  four  horizontal  bands  can  be  divided  into  the  two  upper  bands 
which  present  information  on  activity,  and  the  two  lower  bands  concerning 
listening  times.  The  horizontal  axis  is  the  time  axis,  and  in  order  to 
keep  the  plot  dimensions  reasonable,  each  vertical  line  of  dots  represents 
the  maximum  listening  and  storm  times  for  two  consecutive  Julian  dates. 
If  there  was  any  listening  time  on  either  of  the  two  days  there  will  be  a 
vertical  line  of  dots  in  the  lowest  horizontal  band.  The  band  above  it 
represents  the  maximum  listening  time  for  either  of  the  two  days  by  a 
vertical  line  of  from  0  to  20  dots.  The  number  of  dots  is  directly 
proportional  to  the  amount  of  listening  time  from  0  to  600  minutes.  In 
the  same  way,  the  two  activity  bands  give  first  the  presence  of  activity, 
and  then  a  line  of  0  to  20  dots  proportional  to  the  logarithm  of  the 
maximum  activity  from  0  to  600  minutes.  The  logarithm  was  used  since 
the  amount  of  activity  for  a  given  date  varies  greatly  among  the 
frequencies. 

The  rest  of  this  chapter  will  describe  the  initial  processing  and 
characteristics  of  the  data  base.  First,  the  data  gathered  by  the 
University  of  Florida  will  be  examined. 
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2-2  University  of  Florida  Data 

2.2.1  Initial  Processing 

The  means  of  processing  the  Florida  data  has  not  changed  signifi- 
cantly since  the  early  sixties.  All  of  the  observations  from  Florida 
as  well  as  those  from  Chile  have  been  reduced  according  to  the  method 
described  below.  There  are  some  nonuniformities  in  the  processed  data 
which  will  be  listed  later. 

Signals  from  the  radio  receiving  systems  are  recorded  on  roll  charts 
either  directly  or  through  the  intermediary  of  magnetic  tape.  The 
observer  aurally  distinguishes  Jupiter  storms  from  interference  and  marks 
them  as  such.  Upon  reduction,  any  activity  on  the  records  is  divided  into 
five-minute  intervals  to  facilitate  recording  changing  storm  patterns. 
According  to  the  format  of  Table  II-l,  the  time  and  characteristics  of  the 
storm  within  each  interval  are  punched  on  a  computer  card.  Observation 
times  having  no  Jupiter  activity  are  summarized  on  a  single  card.  As  can 
be  seen,  most  of  the  card  is  used  for  calibration  data  in  order  to  obtain 
flux  density  information  on  the  storm. 

In  1973,  a  640-element  phase-steerable  dipole  array  operating  at 
26.3  MHz  was  completed  at  the  University  of  Florida  (Desch  and  Carr,  1974). 
The  many  differences  between  the  new  array  and  the  older  receiving 
equipment- made  a  different  raw  data  card  format  more  convenient.  Table  1 1-2 
explains  the  format  used  for  that  situation. 

The  actual  determination  of  the  longitude  of  Jupiter  and  the  phase 
of  its  satellites  for  the  listening  and  storm  times,  in  addition  to  the 
flux  density  calculation  for  each  storm  interval,  was  done  by  a  program 
developed  by  H.I.  Register  (1968).  The  program  figures  planetary  and 
satellite  orientation  using  information  from  the  appropriate  day's 
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Table 

n-i 

Raw  Data 

Card 

Florida 

Card 

Year 

Month 

Day 

Listening 

beginning  time  -  hours 

" 

ii     M 

-  minutes 

» 

ending    " 

-  hours 

-  minutes 

Activity  beginning  time 

-  hours 

M     ii 

-  minutes 

ending    " 

-  hours 

ii     ii 

-  minutes 

(1) 


Galactic  background  level  -  percent  of  full  scale 
Larqest  Jupiter  pulse  -  percent  of  full  scale 

Quality  of  data  -  0,  certain;  1,  probable;  2,  possible 
Calibration  input  step  No.  1  -  dB 

"  2  - 

"  3  - 

"  4  - 

"  5  - 

"  6  - 

System  response  to  cal.  step  Ho.  1  -  percent  of  full  scale 

2  - 

3  - 

4  - 

5  " 


Number  of  calibration  steps  (2) 

Cal.  step  matching  output  of  standard  noise  diode  -  dB 

Code  identifying  calibrator  in  use  (3) 

Flux  density  matching  standard  noise  diode  -  Iff22  watts/m2-Hz 

Code  indicating  baseline  of  measurement  (4) 

Code  indicating  meaning  of  date  on  card  (5) 

Difference  between  times  on  card  and  Universal  Time  -  Hours 

Card  code  (6) 

FORTRAN  Format  (1 1 I2,2PF2.0,0PF3.2,I1 ,6F3.1  ,2P6F3. 1 ,11 ,0PF3.2, 
I1,F4.2,3I1,A4) 
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Notes  for  Table  II-1: 

(1)  These  fields  may  be  blank,  indicating  the  listening  time  is  the 
same  as  the  corresponding  activity  time.  This  method  is  used 
for  long  continuous  storms,  usually  in  the  later  years. 

(2)  This  field  has  been  neglected  in  the  later  years. 

(3)  For  further  information  on  calibrators  consult  Miller  (1970), 
Register  (1968),  and  Lebo  (1964). 

(4)  A  0  indicates  measurement  from  the  galactic  level  and  a  1 
signifies  use  of  the  recorder  baseline. 

(5)  For  code  1  observation  periods  extending  through  midnight 
cards  are  often  punched  with  the  second  day's  calendar  date. 
Code  0  cards  should  be  split  at  midnight  to  avoid  this  problem. 

(6)  The  card  code  is  a  four  character  field  containing  information 
on  the  observing  site,  frequency,  and  type  of  antenna  according 
to  the  following  format: 

1st  character  -  observing  site 
F  -  Florida 
H  -  Huanta,  Chile 
K  -  Kodaikanal ,  India 
M  -  Maipu,  Chile 
T  -  Texas 

Y  -  Yale 

2nd  and  3rd  characters  -  frequency 
05  or  5  -  5.6  MHz 
10  -  10.0  MHz 

12  -   12.5  MHz 

13  -   13.1   MHz 

15  -   15.0  MHz 

16  -  16.0  MHz  for  Maipu,  Chile  and  16.7  MHz  for 

Texas  or   Yale 

18  -  18.0  MHz 

19  -  19.2  MHz 

20  -  20.0  MHz 

22  -  22.2  MHz 

23  -  23.0  MHz 

25  -  25.2  MHz 

26  -  26.3  MHz 

27  -  27.6  MHz 
30  -  30.0  MHz 

4th  character  -  antenna  type 
A  -  array 
B  -  broadside 
H  -  dual  full-wave  dipoles  in  an  H  pattern, 

essentially  a  broadside 
I  -  interferometer 
P  -  polarimeter 
S  -  south  half  of  interferometer  only 

Y  -  Yagi 
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Table  II-? 

Raw  Data  Card 

Florida  26.3  MHz  Data 


Month 
Day 

Listening  beginning  time  -  hours 

-  minutes 


ending 

M 

-  hours 

>i                   ii 

" 

-  minutes 

Activity  beginning 

time 

-  hours 

M 

-  minutes 

ii 

-  seconds 

ending 

M 

-  hours 

M 

-  minutes 

H                   n 

ii 

-  seconds 

Character  of  storm  (1 ) 

Credibility  -  1,  certain;  2,  probable;  3,  possible 

Largest  Jupiter  pulse  -  dB 


Galactic  background  level 


dB 


FORTRAN  Format  (212, 4(212, 2X) ,211  ,F4.1  , IX, F4. 1 ) 
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Notes  for  Table  1 1 - 2  : 


1)     Character  of  storm: 

1  -  S  bursts;   spitting 

2  -  L  bursts;   swishing 

3  -  scintillating;   period  greater  than  30  sec, 
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ephemeris  card  (Table  II-3).  Each  raw  data  card  is  converted  on  a  one- 
to-one  basis  into  an  intermediate  data  card  structured  according  to 
Table  II-4.  The  position  of  Jupiter  and  its  satellites  is  expressed 
in  zones  rather  than  degrees.  For  the  cards  using  the  format  of 
Table  II-4  each  zone  represents  a  5°  segment  of  orbit  or  longitude. 
In  later  years,  and  in  reprocessing  some  of  the  data  from  earlier  years, 
the  card  format  was  modified  to  that  of  Table  1 1-5 ,  which  eliminates  some 
unimportant  fields  and  adds  information  on  the  newly-proposed  System  II 
longitude  epoch  1965  (see  Table  1-1).  With  5°  zones,  orientation  can  be 
expressed  in  two  digits,  allowing  more  room  for  other  information.  At 
first,  with  a  small  amount  of  data,  accuracy  within  5°  was  sufficient. 
With  time  the  storm  counts  in  each  zone  accumulated  to  the  point  where 
smaller  zones  would  be  statistically  significant  and  therefore 
desirable. 

Leacock  (1971)  modified  the  Register  program  to  replace  the  5°  zones 
with  2°  zones.  The  modified  intermediate  data  card  (Table  1 1-6)  contains 
less  information  on  the  outer  Galilean  satellites,  but  the  increased 
positional  accuracy  led  to  several  discoveries  which  will  be  discussed 
later.  Why  weren't  1°  or  even  smaller  zones  used,  since  three  card 
columns  would  be  needed  in  each  case?  As  mentioned  before,  storm 
characteristics  are  recorded  in  five  minute  intervals.  In  particular, 
the  flux  density  calculations  are  an  average  for  the  storms  in  each 
interval.  Jupiter  rotates  through  2°  in  approximately  3.3  minutes.  The 
accuracy  of  the  data  reduction,  therefore,  put  a  natural  limit  on  zone 
size,  in  the  range  of  2°.  This  number  is  used  as  a  compromise  giving 
good  positional  information  without  stretching  the  accuracy  of  the 
reduction  process  too  far. 


31 


Table   1 1-3 
Ephemeris  Card 


Year 
Month 


R^Day  (1) 


*h  „m 


Julian  Date  at  0"  0     Universal   Time   (2,3) 


h  nm 


System  III  Longitude  (57.0)  at  0n  0  Universal  Time  (3,4) 


h  nm 


Position  of  Io  from  s.g.c.   at  0     0     Universal   Time   (3,5) 

Position  of  Europa  from  s.g.c.    at  0h  0m  Universal   Time   (3,5) 

Position  of  Ganymede  from  s.g.c.   at  0h  0m  Universal  Time   (3, 5] 
■Light  time  from  Jupiter  to  Earth  -  minutes   (3,6) 
Earth-Sun  angle  as  viewed  from  Jupiter  (3,6) 


IZ^>~ Declination  of  Jupiter  as  viewed  from  Earth  -  degrees   (7) 
; Z^P  -  minutes 

f°p?~Jupiter  transit  -   hours        (7) 
a»vi? —  "  -  minutes 


Florida  sunrise  -  hours 

"    -  minutes 

Chile         -  hours 

-  minutes 


FORTRAN  Format  (3I2,3X,F6. 1  ,1 X ,F4. 1 ,1X ,3(F7. 3, IX) , 12, IX, 
F6. 2, IX, 13, 312, 2(11, 12)) 
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Notes  for  Table  1 1-3: 

(1)  Invalid  day  numbers  such  as  00  and  32  were  included  in  the  early 
years  for  extrapolation  purposes. 

(2)  The  first  two  digits  of  the  Julian  date  are  always  24  and  are 
excluded  in  the  FORTRAN  format. 

(3)  These  fields  are  often  filled  with  zeroes  in  the  early  years 
at  times  near  the  conjunction  of  Jupiter  and  the  sun  or  when 
no  observations  were  made. 

(4)  Prior  to  1961  System  II  longitude  was  used. 

(5)  The  abbreviation  s.g.c.  stands  for  superior  geocentric  conjunction 
These  values  are  uncorrected  for  light  time  through  1967. 

(6)  After  1967  these  fields  are   not  used. 

(7)  These  data  were  included  starting  in  1970.  Times  are  expressed 
according  to  Eastern  Standard  Time. 

The  data  set  containing  all  of  the  ephemeris  cards  (currently  up  to 

1975)  is  at  the  end  of  tape  volumes  JRID1  and  JRID2.  The  data  set 
name  is  EPH  and  the  number  is  48. 
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Table  1 1-4 

Intermediate  Data  Card 

Florida  5°  Zones  -  Early  Years 


■Card  Code  (1 

•Year 
Month 
■Day 


•Julian  Date  -  24 

■System  III  Longitude  (57.0)  -  listening  beginning  zone 

ending 

activity  beginning  " 

ending  " 
Position  of  Io  from  s.g.c.  -  listening  beginning  zone  (2) 

"     "  "    ii    it         M        ..  „   „' 

ending    '   " 

activity  beginning   "   " 
„  jy         "  "      ending    "   " 

~*°S     Sub-Io  System  III  Longitude  (57.0)  -  listening  beginning  zone 
"  °/  _   "     "     ii     ii      M        ,i        j-     I, 

ending    " 

activity  beginning  " 
ending    " 
Position  of  Europa  from  s.q.c.  -  listening  beqinninq  zone  (2,3) 

ending    "    " 
activity  beginning  "    " 

"       II      II        ll        I,  I!  j-     3   ,,       „ 

ending    "    " 

Sub-Europa  System  III  Longitude  (57.0) -  listening  beginninq  zone  (3' 
11        "     m     n       ii       M        i  •     ii 

ending       " 

"    activity  beoinning  "    " 
"     "     n       ii       n        j.     n    „ 

ending       " 

Position  of  Ganymede  from  s.g.c.  -  listening  beqinninq  zone  (2) 

ending    "    " 

activity  beginning  "    " 

ll        M       ii  n        ,,  iii. 

ending    "    " 
Length  of  activity  in  minutes  (4) 


Maximum  flux  density  for  activity  -  10~22  watts/m2-H; 


Quality  of  data  -  0,  certain;  1,  probable;  2,  possible 

Basis  for  angular  measurements  -  0,  geocentric;  1,  heliocentric  (5) 


FORTRAN  Format  (A4, 312, 15, 2412, IX, II ,1X,F7.2,5X,2I1 ) 
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Notes  for  Table  1 1-4: 

(1)  For  card  code,  see  Table  II-l. 

(2)  The  abbreviation  s.g.c.  stands  for  superior  geocentric 
conjunction. 

(3)  In  some  years  data  for  Jovian  moon  V,  Amalthea,  were  substituted 
in  the  Europa  fields.  Table  1 1-8  lists  the  exceptions. 

(4)  Activity  greater  than  9  minutes  often  occurs  in  the  early 
years  or  as  a  result  of  errors  in  the  data.  The  result  is  an 
asterisk  in  column  65.  In  some  of  the  earlier  apparitions  the 
asterisks  have  been  replaced  by  zeroes. 

(5)  All  of  the  data  listed  in  Table  1 1 -8  have  been  reduced  geo- 
centrically. 
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Table  1 1-5 

Intermediate  Data  Card 

Florida  5°  Zones  -  Later  Years 


t-pr-  Year 
*C- Month- 


Card  Code  (1) 


Julian  Date  -  24 

•System  III  Longitude  (57.0)  -  listening  beginning  zone 
|  '^~7~     "  "       "      ending    " 

"i^~7-  "  "     activity  beginning  " 

3  ~y~  ending    ' 

-f^~r- Position  of  Io  from  s.g.c.  -  listening  beginning  zone  (2) 

ending 
activity  beginning   "    " 
ending    "    " 

Sub-Io  System  III  Longitude  (57.0)  -  listening  beginning  zone 

ending    " 
activity  beginning  " 
-J/7-   "  "  "      ending    " 

"U^— Position  of  Europa  from  s.g.c.  -  listening  beginning  zone  (2) 
I  ~y  "     "  "      ending    "    " 

:°V~   "        "     "    "     activity  beginning  " 

t-    o/  II  II  II  II  II  ||  1  •  II  || 

=  <J/  ending 

;°V  Sub-Europa  System  III  Longitude  (57.0)  -  listening  beginning  zone 

I    es/  II  II  II  II  II  II  _J  •  II 

ending 
activity  beginning  " 
ending 

listening  beginning  zone  (2) 

"      ending  "    " 

activity  beginning  "    " 

ending  " 


Position  of  Ganymede  from  s.g.c 


—  8  "'  'I       II       ii        ii       n 

S-f-s^r- Length  of  activity  in  minutes  (3) 

—  System  III  Longitude  (65)  -  listening  beginning  zone 


-22       ^ 
Maximum  flux  density  for  activity  -  10    watts/m  -Hz  (4] 

System  III  Longitude  (65)  -  listening  ending  zone 

activity  beginning  zone 
"     ending    " 


FORTRAN  Format  (A4, 312, 15, 2412, IX, II ,12, F6. 2, 312) 
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Notes  for  Table  II-5: 

(1)  For  card  code,  see  Table  II-l. 

(2)  The  abbreviation  s.g.c.  stands  for  superior  geocentric  conjunction. 

(3)  Activity  greater  than  9  minutes  sometimes  occurs  resulting 
in  an  asterisk  in  column  65. 

(4)  The  flux  density  field  is  smaller  than  the  field  used  in  the 
early  years  format.  A  greater  frequency  of  field  overflows, 
and  therefore  asterisks,  is  to  be  expected. 
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Table  1 1 -6 

Intermediate  Data  Card 

Florida  2°  Zones 


■Card  Code  (1) 

■Year 
-Month 

•Day 

•Julian  Date  -  24 


■System  III  Longitude  -  listening  beginning  zone  (2) 

"        "  "  ii  j  •  ii      n 

ending 

activity  beginning    "   " 
ending     "   " 

■Position  of  Io  from  s.g.c.  -  listening  beginning  zone  (3) 

»     ii   ii   n     n       ii        ,  .     ,,    „ 

ending    "    " 

activity  beginning  "    " 
ending        " 

Sub-Io  System  III  Longitude  -  listening  beginning  zone  (2) 

ending    "    " 

activity  beginning  "    " 

ii     n     n      ii         ,i        ,  . 

ending 


ii    n 


Length  of  activity  in  minutes 

Maximum  flux  density  for  activity  -  10~22  watts/m2-Hz  (4) 


Quality  of  data  -  0,  certain;  1,  probable;  2,  possible 

Basis  for  angular  measurement  -  0,  geocentric;  1,  heliocentric  (4,5] 


FORTRAN  Format  (A4, 312, 15, 1213, 11X, 13,1 X,F7.2,5X,2I1 ) 
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Notes  for  Table  1 1-6 : 

(1)  For  card  code,  see  Table  II-l. 

(2)  These  fields  may  be  in  epoch  1957.0  or  epoch  65  coordinates 
depending  on  the  data  set  used.  See  Table  II-8  for  listing  of 
data  sets. 

(3)  The  abbreviation  s.g.c.  stands  for  superior  geocentric  conjunction. 

(4)  These  fields  are  not  used  in  epoch  65  data. 

(5)  All  of  the  data  listed  in  Table  1 1-8  are  reduced  geocentrically. 
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Desch  developed  the  program  to  produce  the  26-MHz  array  intermediate 
cards  from  the  raw  data  cards.  Again,  the  characteristics  of  the  array 
make  more  desirable  an  intermediate  card  format  differing  from  that  of 
the  other  frequencies  (Table  II-7).  For  example,  the  sensitivity  of 
the  array  can  result  in  storm  times  greater  than  100  minutes,  requiring 
three  card  columns  for  storm  length.  Also,  increased  accuracy  in 
data  reduction  permitted  the  use  of  1°  zones. 

2.2.2  Data  Characteristics 

Almost  20  years  of  observation  have  yielded  over  80,000  raw  data 
cards  and  an  even  larger  number  of  intermediate  cards.  Such  a  profusion 
of  cards  made  multi-year  analyses  quite  difficult.  In  1973,  the  obvious 
step  of  organizing  and  loading  the  data  onto  magnetic  tape  was  completed. 
The  raw  data  cards  have  been  placed  on  tape  volume  JRRD1  (Jupiter  Radio 
Raw  Data  #1).  JRID1  (Jupiter  Radio  Intermediate  Data  #1)  was  the  name 
selected  for  the  intermediate  card  tape  volume.  JRID2  is  an  exact  copy 
of  JRID1. 

The  cards  had  been  produced  by  many  people  over  the  long  span  of 
years.  As  can  be  expected,  they  are  not  uniform  in  their  content  since 
the  information  they  contain  reflects  what  was  considered  important  at 
the  time  of  their  generation.  Table  1 1 -8  summarizes  the  characteristics 
of  the  data  contained  on  the  tapes.  An  overall  summary  for  each  frequency 
is  shown  in  Table  1 1 -9. 

To  access  any  of  the  data  sets  listed  in  Table  1 1 -8  it  is 
necessary  to  know  the  local  tape  file  management  system.  For  the 
IBM  370/165  computer  at  the  Northeast  Regional  Data  Center  (NERDC),  the 
following  three  data  definition  cards  are  suggested  for  reading  from  any 
of  the  three  tape  volumes: 
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Table  I I -7 

Intermediate  Data  Card 

Florida  26.3  MHz  Data  -  1°  Zones 


Card  code  (1 ) 

Year 

Month 

Day 

Julian  Date  -  24 

■System  III  Longitude  (57.0)  -  listening  beginning  zone 

ending    ' 

activity  beginning  " 

"     "     »      ii        ii        .  •     I. 

ending    " 

•Position  of  Io  from  s.g.c.  -  listening  beginning  zone  (2) 

ending 

activity  beginning  "    " 
ending       " 

■Position  of  Europa  from  s.g.c.  -  listening  beginning  zone  (2) 

"     ii    'I     ii     ii        n        .  •      ,,   „ 

ending     "   " 

activity  beginning   "   " 

ending     "   " 

Position  of  Ganymede  from  s.g.c.  -  listening  beginning  zone  (2) 

H  'I  ii  ii  n  n  .  •  „         ,, 

ending         "         " 

activity     beginning     " 

"  "  ii  ii  n  ii  j  •  M         ii 

ending  " 

Length  of  activity  in  minutes 

Maximum  flux  density  for  activity  -  Janskys 

Character  of  storm  (3) 

Credibility  -  1,  certain;  2,  probable;  3,  possible 

Basis  for  angular  measurement  -  0,  geocentric;  1,  heliocentric  (4) 


FORTRAN   Format   (A4, 312, 15,1 713,1 X,F6.0,4X, 311 ) 
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Notes  for  Table  1 1-7: 

(1)  The^card  code  is  F26A.  See  Table  II-l  for  an  explanation  of  the 

(2)  The  abbreviation  s.g.c.  stands  for  superior  geocentric  conjunction. 

(3)  Character  of  storm: 

1  -  S  bursts;  spitting 

2  -  L  bursts;  swishing 

3  -  scintillating;  period  greater  than  30  seconds 

(4)  All  of  the  26.3-MHz  data  listed  in  Table  II-8  are  reduced 
geocentrically. 
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Table  1 1-9 

Summary  of  Frequency  Statistics 

Florida  Data 


DSN 

#  Cards 

Total  List 
(minutes) 

(1) 

(2) 

(3) 

M05 

460 

28815 

M10 

3063 

135994 

H10I 

100 

7363 

F10P 

492 

64958 

F12H 

2163 

74621 

M13P 

680 

56321 

M15 

2211 

121274 

F15Y 

5778 

369284 

M16 

2701 

223281 

M18Y 

9503 

750464 

H18Y 

2291 

91823 

F18Y 

15425 

878970 

M20 

462 

120425 

F20Y 

7300 

658540 

M22Y 

6162 

720419 

F22Y 

10217 

959295 

F25Y 

1222 

130716 

F26A 

936 

59541 

M27Y 

4049 

685555 

F27Y 

4907 

881709 

Totals 

80122 

7019368 

Total  Act. 

Prob 

(minutes) 

(3) 

(4) 

1326 

.046 

11490 

.084 

194 

.026 

686 

.106 

8855 

.119 

2447 

.043 

8191 

.068 

20037 

.054 

8059 

.036 

30780 

.041 

9499 

.103 

56827 

.065 

489 

.004 

23442 

.036 

15159 

.021 

32623 

.034 

2059 

.016 

9494 

.159 

5166 

.008 

7605 

.009 

254428 


.036 


Notes  for  Table  1 1-9: 

(1)  DSN  is  the  name  of  the  data  set  (see  Table  1 1-8)  which  is  a 
generalization  of  the  card  codes  as  explained  in  Table  II-l. 
The  5°  names  are  used  since  the  rest  of  the  table  is  based  on 
these  data. 

(2)  The  number  of  cards  in  the  5°  data  sets  is  listed.  The 
number  of  cards  in  the  corresponding  raw  data  and  2°  data 
sets  will  sometimes  vary  because  of  lost  cards. 

(3)  The  activity  and  listening  times  are  computed  from  the  5°- 
zone  intermediate  cards. 

(4)  The  probability  shown  is  the  average  occurrence  probability  or 
simply  the  ratio  of  total  activity  time  to  total  listening  time. 
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//filename  DD  UNIT=(TAPE9, , DEFER) ,DISP=0LD,V0L=( , RETAIN, SER=vname) , 
//         DCB=(RECFM=FB,LRECL=80,BLKSIZE=3200), 
//        LABEL=(ds#,,,IN),DSN=dsn 

The  appropriate  tape  volume  name  is  substituted  for  vname  and  ds#  and 
dsn  are  listed  in  Table  II-8. 

Naturally,  there  is  a  need  for  adding  new  data  to  the  tape  as  it 
becomes  available.  The  suggested  system  for  doing  this  is  shown  in 
Figure  1 1-4.  JRID2  is  the  volume  and  serial  number  of  the  main  tape 
kept  in  the  NERDC  tape  library.  JRID1  is  the  backup  tape  to  be  kept  by 
the  person  in  charge  of  the  data.  The  program  STLMY  performs  the  conversion 
from  raw  to  intermediate  data  in  either  2°  or  5°  format.  The  new  inter- 
mediate data  are  generated  on  a  temporary  disk  data  set  and  listed  on  the 
printer.  The  disk  is  merged  with  the  old  tape,  JRID2,  to  yield  an  updated 
JRID1  tape  volume.  After  the  resulting  tape  has  been  carefully  verified 
for  accuracy,  it  is  copied  onto  volume  JRID2.  Note  that  JRID1  and  JRID2 
currently  contain  ephemeris  cards  for  1957-1975,  which  are  necessary  for 
processing  raw  data  (see  Table  1 1-3) .  If  ephemeris  cards  later  than  1975 
are  needed,  either  the  ephemeris  file  on  tape  will  have  to  be  augmented  or 
the  ephemeris  cards  should  be  read  in  directly.  In  the  same  way,  raw  data 
may  be  put  in  on  cards  or  merged  with  tape  volume  JRRD1  before  use. 

2.3  University  of  Texas  Data 

Thanks  to  the  generosity  of  the  decametric  radio  observers  at  the 
University  of  Texas,  the  University  of  Florida  received  a  catalog  of  the 
Texas  data  (Bozyan,  et  al . ,  1972)  at  the  same  time  the  Florida  data  were 
being  organized  for  tape  loading.  The  texas  observing  frequencies  were 
often  close  to  or  equal  to  those  of  Florida,  so  that  a  merger  of  the 
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Figure  1 1-4.  System  flowchart  for  update  of  intermediate  data 
tapes. 
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two  data  sets  was  obviously  advantageous.  Raw  data  cards  punched  according 
to  Table  11-10  were  converted  to  intermediate  cards  of  the  form  shown  in 
Table  11-11  by  a  modification  of  the  Florida  reduction  method.  Note 
that  the  accuracy  of  the  storm  and  listening  times  was  deemed  sufficient 
to  justify  expressing  the  positions  of  Jupiter  and  its  satellites  in  1° 
zones. 

Besides  the  Florida  data,  the  magnetic  tapes  JRID1  and  JRID2  also 
contain  the  Texas  intermediate  cards.  However,  raw  Texas  data  have  not 
been  loaded  onto  the  JRRD1  tape.  Some  overall  characteristics  of  the 
Texas  observations  are  found  in  Table  11-12.  For  more  details  the  reader 
should  consult  the  catalog. 

In  conclusion,  the  Texas  and  Florida  intermediate  cards  are  the 
foundation  for  this  dissertation.  The  large  amount  of  data  is  especially 
important  for  the  higher  resolution  studies  to  be  presented  in  the  next 
few  chapters. 
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Table  11-10 

Raw  Data  Card 

Texas  Data 


Month 

Day 

Listening  beginning  time 

-  hours 

-  minutes 

ending    " 

-  hours 

-  minutes 

Activity  beginning  time 

-  hours 

-  minutes 

ending    " 

-  hours 

-  minutes 

Storm  type  (1) 

Storm  strength  (2) 

Storm  credibility  (3) 

FORTRAN  Format  (1012,311) 
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Notes  for  Table  11-10: 

(1 )  Storm  type: 

1  -  Main  storm  in  observing  time 

2  -  Precursor  in  observing  time 

3  -  Post-cursor  in  observing  time 

4  -  Main  storm  outside  observing  time 

5  -  Precursor  outside  observing  time 

6  -  Post-cursor  outside  observing  time 

(2)  Storm  strength: 

Blank  -  Not  specified 

1  -  Strong  (peak  flux  greater  than  20  times  Cygnus  A) 

2  -  Moderate  (peak  flux  10-20  times  Cygnus  A) 

3  -  Weak  (less  than  10  times  Cygnus  A) 

(3)  Storm  credibil ity: 

1  -  Certain  (probability  event  is  Jovian  99%) 

2  -  Probable  (    "       "    "    "   90%) 

3  -  Possible  (    "       "    "    "   50%) 

4  -  As  in  1,  but  may  be  solar  emission 

5  -  "  "  2,  " 

6  -  "  "  3,  ' 


11  11 


All  of  the  above  are  explained  in  greater  detail  in  the  cataloq  by 
Bozyan  et  a! .  (1972). 


Table  11-11 

Intermediate  Data  Card 

Texas  -  1°  zones 


^N-Card  code  (1) 


Year 

Month 

Day 


Julian  Date  -  24 


System  III  Longitude  (57.0)  -  listening  beginning  zone 

"     "     "      "        »       j  •      ii 

ending     " 

activity  beginning   " 

"     "     "      »        ii       j  •      ii 

ending    " 

Position  of  Io  from  s.g.c.  -  listening  beginning  zone  (2) 

"     "  "    "    "        H        j  •      ii    ii 

ending    "    " 

activity  beginning  "    " 

ending    "    " 

Position  of  Europa  from  s.g.c.  -  listening  beginning  zone  (2) 

"     "   "     "     "         n       j  •     i, 

ending    "    " 

activity  beginning     "         " 

"  "       "  "  ii  ii  j  •  M         ii 

ending         "         " 

Position  of  Ganymede  from  s.g.c.   -  listening  beginning  zone  (2] 

"  "       "  "  ii  n  j  •  ii         ii 

ending         "         " 

activity  beginning  "    " 


ending 


H         ii 


Storm  type  (3) 

Storm  strength  (3) 

Storm  credibility  (3) 

Activity  beginning  time  in  Universal  Time  -  hours  and  minutes 

ending    "   "     "      "     "     <>         « 
Length  of  activity  in  minutes 


FORTRAN  Format  (A4, 312, 15, 1613, 311 ,214,13) 
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Notes  for  Table  11-11: 

(1)  For  explanation  of  card  code,  see  Table  II-l. 

(2)  The  abbreviation  s.g.c.  stands  for  superior  geocentric 
conjunction. 


(3)  For  explanation  of  these  fields,  see  Table  II-l 
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Table  11-12 

Summary  of  Frequency  Statistics 

Texas  Data 


DSN 

DS# 

#  Cards 

Total  List, 
(minutes) 

Total  Act. 
(minutes) 

Prob. 

(1) 

(2) 

(3) 

(4) 

(4) 

(5) 

Y10I 

1 

169 

22638 

3795 

.168 

T10I 

2 

73 

4672 

147 

.031 

Y16I 

3 

330 

41815 

7001 

.167 

T16I 

4 

724 

78167 

7542 

.096 

Y19I 

5 

117 

9592 

514 

.054 

Y20I 

6 

2509 

437106 

43727 

.100 

T20I 

7 

896 

107993 

7388 

.068 

Y22I 

8 

2755 

570155 

37945 

.067 

K22I 

9 

782 

125931 

8285 

.066 

T22I 

10 

1219 

164561 

7427 

.045 

Y23I 

11 

263 

28931 

1067 

.037 

T30I 

12 

1307 

207296 

1883 

.009 

Totals 


11144 


1798857 


126721 


.070 


Notes  for  Table  11-12: 

(.1)  DSN  is  the  data  set  name  to  be  used  in  the  DSN=  parameter  of 
the  data  definition  cards  for  magnetic  tape  volumes  JRID1  and 
JRID2.  It  is  a  generalization  of  the  card  code  which  is  explained 
in  Table  1 1  - 1 . 

(2)  The  data  set  number  in  this  column  is  used  in  the  LABEL= 
parameter  of  the  data  definition  cards  for  magnetic  tape 
volumes  JRID1  and  JRID2. 

(3)  Listed  are  the  number  of  cards  in  the  data  set  on  magnetic 
tape  volumes  JRID1  and  JRID2. 

(4)  The  activity  and  listening  times  are  computed  from  the 
intermediate  cards  on  magnetic  tape. 

(5)  The  probability  shown  is  the  average  occurrence  probability 
or  simply  the  ratio  of  total  activity  time  to  total  listening 
time. 


CHAPTER  III 
ONE-DIMENSIONAL  STUDIES 

3.1  Histogram  Accuracy 

3.1.1  Statistical  Significance 

The  usual  method  of  categorizing  Jupiter's  discrete  radio  events 
has  been  the  histogram,  such  as  that  of  Figure  1-1,  in  which  the  data 
are  related  to  a  single  variable.  In  this  sense,  the  radio  storms  are 
studied  for  single-dimensional  dependence.  Generally,  the  rotation  of 
Jupiter  or  the  orbit  of  its  satellites  is  divided  into  intervals,  and 
the  number  of  events  or  the  occurrence  probability  (the  number  of  events 
divided  by  the  total  number  of  observations)  is  plotted  for  each  zone. 
The  representation  is  straightforward,  but  conclusions  based  on  it  must 
consider  the  statistical  significance  of  the  method. 

Let  us  look  at  the  typical  occurrence  probability  histogram.  The 

usual  method  of  computation  for  Jupiter's  decametric  probability 

histograms  is  to  set 

n 
p(x,Ax)  =  I     a.(x,Ax)/n  (III-l) 

i=l  1 

where  x  is  the  System  III  longitude,  Galilean  moon  phase,  declination  of 
the  earth,  or  any  other  appropriate  independent  variable.  The  probability 
in  the  interval  from  x  to  x+Ax  is  symbolized  by  p(x,Ax).  A  value  of  1 
or  0  for  the  a^x.Ax)  represents,  respectively,  the  presence  or  absence 
or  a  radio  storm  for  each  of  n  observations  within  the  interval.  The 
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main  question  is;  Do  the  histograms  closely  approximate  the  functional 
curve  expected  from  observing  for  an  infinite  amount  of  time  and 
computing  the  probability  for  infinitely  small  intervals?  Restating 
the  question  mathematically,  if 

n 

p(x)  =  lim  I   a.(x,Ax)/n  (III-2) 

n-x»  i=l 

Ax->-o 

then  is  p(x)=p(x,Ax)?  There  are  two  problems  in  answering  this  question. 
First,  in  any  interval,  how  much  statistical  variation  can  be  expected 
with  a  number  of  observations  substantially  less  than  infinite?  Second, 
assuming  no  statistical  fluctuation,  is  Jupiter's  behavior  sufficiently 
repetitive  that  the  same  curve  will  be  continually  retraced? 

Consider  the  first  problem.  A  multitude  of  factors  affect  the 
number  of  radio  storms  in  any  interval.  Several  critical  examples  are  the 
sensitivity  of  the  radio  telescope,  the  transparency  of  the  earth's 
ionosphere,  and  the  separation  of  Jupiter  emissions  from  non-Jovian 
interference.  There  will  also  be  minor  influences  such  as  the  relative 
state  of  repair  of  the  detection  equipment  and  the  human  error  involved 
in  identifying  the  beginning  and  ending  of  storms  and  in  codifying  the 
data.  Of  course,  the  exact  influence  of  all  of  these  factors  cannot  be 
calculated.  The  assumption  will  be  made  that  their  contribution  is  a 
random  fluctuation  about  the  limiting  functional  curve.  Thus,  if  the 
number  of  independent  observations  is  large,  the  standard  deviation  can 
be  estimated  as  the  square  root  of  the  number  of  events  in  any  interval: 


0  %L   ai(x'Ax>  (II 1-3) 

Ap(x,Ax)  =  a/n  (II 1-4) 
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The  second  equation  gives  the  corresponding  variation  in  the  probability. 
The  credibility  of  this  estimate  is  directly  related  to  the  amount  of 
data. 

It  is  possible  that  statistical  variation  may  be  blamed  for  what 
is  actually  a  natural  variation  in  the  functional  curve  with  time.  To 
what  extent  are  Jovian  emission  events  predictably  repeatable?  This 
requires  a  knowledge  of  all  independent  variables  which  influence  the 
emission.  It  is  known  that  at  middle  and  high  decametric  frequencies 
occurrence  probability  versus  System  III  longitude  histograms  would 
change  drastically  for  different  ranges  of  the  position  of  Io  from 
superior  geocentric  conjunction  (hereafter  called  the  "phase"  of  Io  and 
symbolized  by  <t>T0)-  Combining  radiation  from  all  of  the  positions  of  Io 
brings  a  measure  of  stability  to  the  source  probabilities,  but  they  are 
still  altered  in  a  timespan  of  years  due  to  the  changing  declination  of 
the  earth  as  viewed  from  Jupiter  (Dp).  Therefore,  the  probability  of 

receiving  a  radio  storm  depends  on  at  least  three  variables,  XTTT,  d>T  , 

III  rIo 

and  Dp.  Since  only  about  two  cycles  of  the  variable  Dp  have  been 
observed,  there  is  incomplete  coverage  of  all  of  their  possible  com- 
binations. 

Bozyan  and  Douglas  (1975)  consider  this  problem  in  presenting  the 
results  of  their  decametric  studies.  They  assume  the  probability  to  be 
a  function  of  the  three  variables  such  that  at  any  time,  t, 
p(t)=p(Ajjj,c|>Io,Dp).  To  get  a  good  estimate  of  this  probability  function 
Bozyan  and  Douglas  state  the  need  for  simultaneous  observations  at  many 
terrestrial  longitudes,  and  the  longest  possible  runs  of  synoptic 
monitoring.  Obviously,  simultaneous  observations  will  improve  the 
calculation  of  the  occurrence  probability  for  a  given  instant  of  time. 
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To  assume  that  continuous  monitoring  is  beneficial  assumes  also  that 
the  probability  curve  is  cyclic.  How  can  one  be  certain  that  the 
occurrence  probability  is  a  function  of  xin»   <)>,  ,  and  D£  only,  such 
that  several  sets  of  observations  of  a  given  range  of  functional 
values  will  yield  the  same  probability?  The  chances  are  that  the 
probability  does  depend  on  other  variables,  but  the  data  so  far  indicate 
first-order  dependence  on  only  \jjji  <|>l0.  and  DV.  The  simplest  model  is, 
then,  only  dependent  on  these  three. 

This  dissertation  is  based  on  a  search  for  high-resolution,  second- 
order  effects.  The  term  "second-order"  is  used  loosely  to  classify 
those  phenomena  which  are  not  immediately  obvious  upon  inspection  of 
the  data,  and  which  are  barely  separable  from  statistical  fluctuation 
"noise."  Although  these  phenomena  may  be  more  difficult  to  detect,  the 
search  for  them  should  not  be  abandoned.  More  care  must  be  taken  in 
their  identification.  They  may  be  related  to  the  variables  already 
mentioned,  in  which  case  they  would  recur  with  the  appropriate  periods, 
or  they  may  have  entirely  different  periods  or  no  period  at  all.  Since 
repetition  of  effects  cannot  be  anticipated,  obtaining  the-maximum  number 
of  independent  observations  for  each  instant  of  time  becomes  important. 
Simultaneous  observations  from  several  terrestrial  longitudes  have  already 
been  mentioned.  This  dissertation  will  also  merge  proximate  frequency 
data  obtained  at  the  same  observing  site  from  different  radio  telescopes. 
Neither  of  these  methods  yields  truly  independent  observations,  since 
changes  in  the  propagating  medium  at  the  earth,  Jupiter,  or  anywhere 
between  may  affect  all  simultaneous  observations.  In  general,  the  second- 
order  effects  will  be  assumed  to  persist  from  several  months  to  several 
apparitions,  so  that  with  time  the  short-term  fluctuations  will  be 
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eliminated  by  averaging.  In  cases  where  the  effects  are  periodic, 
long-time-span  analyses  will  be  beneficial. 

3.1.2  Time  Vs.  Event  Histograms 

In  computing  the  occurrence  probability  for  a  histogram  interval 

according  to  Equation  III-l,  accuracy  is  sometimes  lost.  Within  a 

particular  interval  the  listening  time  often  exceeds  the  duration  of 

activity,  yet  the  occurrence  probability  will  still  be  computed  as  unity. 

To  avoid  this,  the  time  of  activity  should  be  divided  by  the  listening 

time  for  each  histogram  zone. 

n  m 

p(t,x,Ax)  =  I     a.(t,x,Ax)  /  I   1  (t,x,Ax)  ( 1 1 1-5) 

1=1  1         j=l  J 

Here  a.(t,x,Ax)  represents  the  storm  time  for  the  ith  event  within  the 
histogram  zone  x  to  x+Ax.  The  listening  time  for  the  jth  observation 
of  the  zone  is  represented  by  1  (t,x,Ax).  Since  more  than  one  storm  can 
occur  during  a  span  of  listening,  these  variables  cannot  be  included  in 
the  same  summation.  For  a  large  number  of  observations  occurrence 
probability  histograms  computed  from  either  Equation  III-l  or 
Equation  II 1-5  appear  similar,  although  the  individual  probabilities 
computed  from  the  former  equation  will  be  proportionally  higher.  For 
accurate  probabilities  and  when  data  are  scarce,  Equation  1 1 1-5  is  more 
advantageous.  This  dissertation  will  use  histograms  based  on  time 
calculations  exclusively. 

Time-based  nistograms  present  some  new  difficulties.  Usually, 
the  beginning  and  ending  listening  times  for  a  particular  night  are  known 
only  approximately.  Since  only  the  beginning  and  ending  listening  zones 
are  given  on  the  intermediate  data  card  (see  Table  II-4),  for  ease  of 


63 


computation  the  total  computed  listening  time  is  expanded  to  correspond 
to  the  time  between  entering  the  first  zone  and  exiting  the  last. 
This  will  slightly  lower  the  probabilities  in  the  endpoint  intervals, 
but  the  error  decreases  as  zone  size  shrinks,  and  it  affects  only  two 
zones  for  each  observation  period.  For  activity,  the  intermediate  data 
card  usually  has  information  only  on  the  storm  length  and  the  beginning 
and  ending  zones.  The  storm  time  is  adjusted  so  that  the  activity  extends 
equally  into  the  beginning  and  ending  intervals.  The  error  caused  by 
this  should  average  out  as  other  storms  are  added  to  the  zones,  and  it 
decreases  with  a  reduction  in  interval  size.  For  the  Florida  data 
continuous  radio  storms  are  usually  split  into  five-minute  increments. 
These  divisions  are  rejoined  before  they  are  processed  in  the  manner 
described  above. 

For  time-based  histograms  the  error  estimate  described  in 
Equation  III -4  no  longer  applies.   It  is  difficult  to  figure  any  error 
without  knowledge  of  the  average  time  contribution  per  storm  to  each 
zone.  Consequently,  the  number  of  events  in  each  interval  has  still 
been  recorded.  With  this  information,  the  standard  deviation  of  the 
probability  can  be  computed  as: 

n  n  n 


o  -  a  I     a.(t,x,Ax)/  I   a.(x,Ax)  =  J  a .  (t,x,Ax)/ "1  7  a  .  (x,Ax)  (II 1-6) 
z         1=1  ^  j=l  J       i=l  }  [j=l  J 

m 
Ap(t,x,Ax)  =  oj   I   1.  (t,x,Ax)  (III-7) 

1   k=l  k 

where  a  is  the  standard  deviation  from  Equation  1 1 1 -3  and  the  second 
equation  is  the  error  estimate  for  the  probability.  The  error  estimate 
of  Equation  1 1 1 -  7  will  be  lower  than  that  defined  in  Equation  III-4,but 
the  time-based  probability  will  also  be  lower  than  its  event-based 
equivalent. 
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3.1.3  Decreasing  the  Histogram  Interval  Size 

The  previous  chapter  introduced  the  obvious  method  for  increasing 
histogram  resolution  in  presenting  the  smaller  2°- histogram- interval 
Florida  data.  The  discoveries  made  by  Leacock  (1971)  proved  the  value 
of  reprocessing  the  Florida  raw  data  into  2°  zones.  These  findings 
will  be  presented  later,  but  it  will  suffice  to  say  that  a  subdivision 
of  the  histogram  interval  is  necessary  for  observation  of  some  second- 
order  effects. 

How  much  subdivision  can  be  tolerated  before  the  data  within  each 
zone  are  too  sparse  to  make  the  results  believable?  This  is  a  matter  of 
opinion.  One  measure  of  the  accuracy  is  the  uniformity  of  listening  time 
in  the  histogram  intervals.  Given  a  reasonable  number  of  observation 
periods  at  random  observing  times,  the  amount  of  listening  time  in  each 
of  the  intervals  should  be  approximately  equal.  The  percent  standard 
deviation  of  the  group  of  interval  listening  times  has  been  monitored 
as  a  test  of  the  uniformity. 

3.2  Decametric  Spectral  Analysis 

3.2.1  System  III  Longitude 

The  easiest  way  to  display  some  important  aspects  of  Jupiter 
emission  throughout  the  decametric  wavelength  range  is  to  construct  a 
three-dimensional  plot  of  occurrence  probability  for  all  System  III 
longitudes  over  the  appropriate  span  of  frequencies.  Figure  III-l  is 
an  attempt  to  simulate  a  diagram  of  this  nature  by  combining  System  III 
longitude  histogram  data  for  the  individually-observed  frequencies  having 
more  than  100,000  minutes  of  observing  time.  Ninety  probability  values, 
obtained  by  averaging  occurrence  probabilities  for  odd-even  pairs  of 
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zones  in  180-zone  histograms,  are  plotted  for  each  frequency.  The 
5.0-MHz  data  are  below  the  listening  time  cutoff  and  span  only  a 
limited  number  of  years,  but  they  are  included  as  an  indication  of  the 
lowest  frequency  results  obtainable  without  satellites.  Values  of 
occurrence  probability  can  be  determined  by  comparing  the  vertical 
scale  at  the  right  of  the  figure  with  the  height  of  the  curve  above  the 
zero  plane.  It  should  also  be  noted  that  the  Florida  16.0-MHz  and 
Texas  16.7-MHz  data  have  been  combined  and  labeled  as  16.0. 

In  this  figure  it  is  obvious  that  sources  A,  B,  and  C,  as  delineated 
in  Figure  1-1,  are  well-defined  in  the  middle  decametric  frequencies. 
At  first  glance  it  appears  that  at  15  MHz  and  above  the  plots  are  like 
the  profile  of  an  island  which  is  gradually  sinking  beneath  the  sea  level 
of  the  zero  plane.  This  is  consistent  with  the  fact  that  for  similar 
receiving  equipment  Jovian  activity  decreases  as  frequency  increases 
(Carr  et  al . ,  1964).  Closer  inspection  discloses  a  disproportional 
decrease  in  the  size  of  source  C  above  15  MHz.  The  same  is  true  for 
source  A  above  22  MHz.  Some  process  more  complicated  than  an  overall 
decrease  in  activity  must  be  involved.  The  source-B  region  appears  to 
decrease  gradually  and  then  increase  at  30  MHz.  This  upward  trend  is 
an  artifact  resulting  from  data  acquired  solely  from  the  Texas-based 
interferometer  used  at  that  frequency.  All  of  the  other  frequencies 
involve  single  Yagi  antennas,  or  a  mixture  of  Yagis  and  interferometers, 
providing  lower,  more  comparable  sensitivities. 

Below  15  MHz  the  sources  no  longer  maintain  their  integrity.  Two 
peaks  are  prominent  at  the  10-MHz  level,  and  the  5-MHz  data  are  too 
erratic  for  determination  of  any  trend.  The  exact  metamorphosis  of 
the  plot  from  15  MHz  to  10  MHz  and  beyond  would  be  enlightening,  but 
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lack  of  sufficient  data  in  this  void  makes  interpolation  necessary. 
Figure  II 1-2  suggests  how  the  longitude  positions  of  the  source  pro- 
bability peaks  change  between  10  and  30  MHz.  Similar  graphs  have  been 
made  before  (Smith,  et  al . ,  1965;  McCulloch  and  Ellis,  1966).  Data 
at  5  MHz  strongly  influenced  these  earlier  investigators,  but  they 
have  been  rejected  from  Figure  1 1 1-2  as  too  sparse  to  be  given  credence. 
McCulloch  and  Ellis'  results  were  based  on  relative  flux-density  plots 
and  imply  that  sources  A  and  B  fade  away  below  10  MHz  while  a  new  source, 
positioned  around  180°,  appears  near  15  MHz  and  continues,  as  does 
source  C  at  approximately  350°,  through  5  MHz  and  even  lower.  The 
occurrence-probability-based  plots  of  Smith  et  a! .  suggest  that  at  5  MHz 
sources  B  and  C  combine  into  a  single  source  located  near  270°  and  source 
A  moves  to  about  120°.  Two  modes  of  evolution  into  this  lower-frequency 
configuration  were  proposed.  Both  groups  of  investigators  used  only  two 
or  three  years  of  observations  from  the  early  sixties.  In  Figure  III-2, 
the  lines  drawn  through  the  source  peaks  give  another  interpretation 
based  on  (a)  updated  values  for  the  locations  of  the  maxima,  (b)  polarization 
information  (Kennedy,  1969),  and  (c)  two-dimensional  plots  to  be  presented 
in  the  next  chapter.  The  10-MHz  peak  located  at  70°  corresponds  to  the 
appearance  of  source  D.  This,  too,  will  be  more  evident  in  Chapter  IV. 
What  causes  the  changes  in  peak  locations  and  the  suggested  merger 
of  sources  A  and  B  at  the  lower  frequencies?  This  and  other  theoretical 
questions  will  be  saved  until  Chapter  V,  when  the  results  of  Chapters  III 
and  IV  can  be  treated  together. 
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3.2.2  Phase  of  Io 

Bigg  (1964)  proved  the  control  of  Io  on  certain  aspects  of 
Jupiter's  emissions  by  plotting  occurrence  probability  histograms 
for  the  phase  of  Io.  This  revealed  the  now-famous  intensification  of 
occurrence  probability  at  90°  and  240°  phase  values.  Additional  two- 
dimensional  maps,  like  that  of  Figure  1-2,  disclosed  the  responsibility 
of  source  B  for  the  90°  peak  and  the  combination  of  sources  A  and  C 
in  the  240°  enhancement. 

To  study  the  effects  of  Io  on  an  individual  source,  one  must  isolate 
the  storm  events  or  parts  thereof  occurring  in  the  longitude  range  defined 
for  that  source.  The  longitudinal  limits  for  the  sources  have  been 
defined  as  follows:  source  B,  65°-200°;  source  A,  201°-290°;  and 
source  C,  the  remainder.  The  storms  associated  with  each  source  can  then 
be  analyzed  in  histograms  of  occurrence  probability  vs.  Io  phase. 
Figures  III-3,  1 1 1-4,  and  1 1 1 -5  add  information  on  the  variation  of  Io's 
influence  with  frequency.  Like  Figure  III-l,  Figures  III-3,  III-4,  and 
III-5  are  pseudo-three-dimensional,  smoothed  plots  derived  from  histogram 
data. 

Figure  1 1 1 -3  emphasizes  the  strength  of  the  Io  control  of  source  B. 
The  maximum  occurrence  probability,  as  seen  on  the  scale  to  the  right, 
is  0.43  which  more  than  doubles  the  maximum  of  Figure  III-l.  The  90° 
peaks  predominate  through  all  frequencies  except  5  MHz,  which  has  no 
activity  in  the  90°  region.  Again,  the  lack  of  data  (52  per  cent 
standard  deviation  of  the  listening  times  per  zone)  and  limited  number 
of  years  of  observation  at  this  frequency  lessens  the  significance  of  this 
finding.  The  peaks  remain  very   stable  in  position  through  all  of  the 
frequencies  at  a  value  usually  quoted  as  90°,  although  the  average  peak 
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position  for  the  nine  frequencies  was  93°  with  a  standard  deviation 
of  2°.  Notice  also  the  increased  non-Io-related  radiation  as  the 
frequency  decreases. 

Figure  III-4,  a  diagram  of  the  isolated  A  region,  displays  much 
more  prominent  non-Io-related  radiation,  but  the  Io-related  peaks  are 
still  obvious.  Io-A  is  more  frequency  dependent  than  B.  The  maximum 
Io-related  probability  occurs  around  20  MHz,  and  seems  to  taper  off  in 
both  frequency  directions  until  it  is  no  longer  evident  at  10  and 
30  MHz.  Between  these  limits,  the  peaks  have  an  average  position  of 
247°  with  a  standard  deviation  of  6°.  Presumably,  the  large  mixture  of 
non-Io-related  radiation,  especially  at  lower  frequencies  where  the  peak 
position  variation  is  greater,  accounts  for  the  increase  in  standard 
deviation  over  source  B. 

The  quick  disappearance  of  the  Io  influence  on  source  C  above 
22  MHz  is  displayed  in  Figure  III-5.  Unlike  sources  A  and  B,  however, 
Io  control  still  seems  to  dominate  down  to  5  MHz  and  perhaps  even  lower, 
although  the  lack  of  data  at  5  MHz  makes  that  peak  suspicious.  The 
average  Io  phase  for  these  maxima  is  the  same  as  that  of  source  A,  247°. 
The  smaller  standard  deviation  of  3°  follows  from  the  lower  "noise" 
level  of  events  unrelated  to  Io. 

The  frequency  dependence  of  the  maximum  occurrence  probabilities 
in  Figure  III-l  and  Figures  1 1 1 -  3  to  III-5  may  be  more  precisely  examined 
in  Figure  III-6.  The  error  in  the  ordinate  values  is  difficult  to 
determine,  since  data  obtained  by  equipment  with  different  sensitivities 
have  often  been  merged.  Except  for  the  30-MHz  data,  the  Florida 
observations  outweigh  those  from  Texas  and  the  sensitivities  seem  to 
agree  well  enough  for  general  statements  based  on  maximum  probability 
tendencies  to  be  made. 
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Figure  II 1-6.     Spectral  distribution  of  the  maximum  occurrence 

probabilities  in  the  Io-phase  histograms  for  isolated 
source-B,  A,  and  C  regions  and  in  the  Ajjj   histograms 

for  all   sources.     (See  Figures  III-l   and  III-3  to   III-5.) 
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As  expected,  the  maximum  occurrence  probabilities  in  Figures  III-l 
and  1 1 1-4  change  in  a  similar  manner  with  a  change  in  frequency  since 
the  major  contribution  to  the  peaks  in  each  case  comes  mainly  from 
the  same  source-A  Io-related  radiation.  From  10  to  20  MHz  both  sets 
of  maxima  retain  nearly  constant  or  slowly  increasing  levels.  Beyond 
20  MHz  the  probabilities  decrease  at  an  almost  steady  rate  except  for 
the  30-MHz  values  which  are  artificially  high.  Io-dependent  source-B 
maxima  follow  the  same  pattern  except  for  the  slower  rate  of  decline  at 
higher  frequencies.  The  isolated  source-C  region  exhibits  interesting 
differences.  Its  decline  begins  around  15  MHz,  and  the  source  is  just 
barely  noticeable  above  22  MHz.  The  labeling  of  the  barely-perceptible 
storms  at  these  higher  frequencies  as  Io-related  is  questionable.  Note 
that  the  5-MHz  peak  probability  has  been  included  for  source  C  since 
the  peak  is  similar  to  and  agrees  well  with  those  at  higher  frequencies. 
From  the  three  peaks  at  15  MHz  and  below  the  maxima  can  be  construed  as 
approximately  constant  to,  and  possibly  through,  5  MHz.  Dominant 
Io-related-C  radiation  at  frequencies  below  5  MHz  would  be  an  interesting 
discovery.  To  date  no  one  has  published  evidence  positively  linking 
low-frequency  Io-related  radiation  to  a  particular  source. 

3.3  Source  Chronology 


3.3.1  System  III  Longitude  Source  Motion 

Modification  of  the  sources  with  frequency  has  been  discussed, 
but  their  time  evolution  must  also  be  considered.  The  discovery  that 
Jupiter's  activity,  and  especially  the  source-A  characteristics,  change 
over  a  period  of  years  was  made  soon  after  regular  Jupiter  monitoring 
began  (Gallet,  1961;  Douglas  and  Smith,  1963).  The  two  effects  most 
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quickly  recognized  were  a  slow  change  in  the  source-A  occurrence 
probability  and  a  gradual  translation  of  its  A„,  centroid.  Gulkis 
and  Carr  (1966)  first  recognized  the  shifting  of  the  position  of  source 
A  as  a  linear  drift,  resulting  from  a  slight  error  in  the  definition  of 
the  XTTT  period,  combined  with  a  sinusoidal  offset.  The  pure  sinusoidal 
effect  can  be  observed  once  the  data  have  been  reduced  using  the  corrected 
period  defined  by  AJTT(65)  and  the  centroid  of  the  source  is  plotted 
for  each  apparition. 

Calculation  of  the  centroid  of  a  source  is  \/ery   susceptible  to 
error  from  random  statistical  fluctuation  in  the  histogram  probabilities. 
For  this  reason,  the  results  presented  in  the  rest  of  this  chapter  will 
be  based  upon  the  largest  subset  of  the  data  base  considered  meaningful: 
the  merged  observations  for  18,  20,  and  22  MHz.  Even  with  this  pool  of 
more  than  90,000  hours  of  data,  the  lower-activity  years  show  suscepti- 
bility to  statistical  error.  An  estimation  of  the  maximum  error  in  the 
computed  centroid  value  will  be  made  by  the  following  method.  Once  the 
centroid  has  been  computed  from  a  particular  range  of  histogram 
occurrence  probabilities,  the  probabilities  to  the  left  of  the  centroid 
zone  will  be  decreased  by  an  amount  determined  by  Equation  III-7.  The 
probabilities  to  the  right  of  the  centroid  are  increased  in  the  same 
manner  and  a  new  centroid  is  computed.  The  new  value  is  used  as  an 
upper  error  limit  for  the  original  centroid.  The  reverse  process  yields 
a  lower  limit. 

With  this  method  an  updated  study  has  been  made  of  the  source-A 
sinusoidal  motion.  The  gradual  change  in  the  source  is  apparent  from 
Figure  III-7,  which  compares  the  XTTT  histograms  from  each  of  the 
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Figure  1 1 1-7.  Individual-apparition  occurrence  probability  histograms 
for  2°  increments  in  System  III  lonqitude  (65)  based  on 
merged  18,   20,  and  22-MHz  data. 
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apparitions.  Register  and  Smith  (1969)  and  Goertz  (1971)  found  that 
the  variation  was  more  prominent  in  the  Io-dependent  part  of  the 
source-A  radiation,  and  Figure  III -8  confirms  that  this  continues 
through  the  1974-75  apparition.  The  source-A  centroid  positions  and 
accompanying  errors  have  been  plotted  for  each  apparition  for  the  Io- 
dependent  radiation  region,  the  non-Io-related  region,  and  the  combination 
of  the  two.  All  source-A  radiation  occurring  when  the  phase  of  Io  was 
between  180°  and  270°  was  defined  to  be  Io-related.  The  centroids  were 
computed  from  data  within  particular  boundaries  defined  for  each  of  the 
Io-dependent,  Io-independent,  and  overall  regions.  The  boundaries  were 
determined  from  histograms  for  each  region  combining  all  apparitions. 
The  minima  on  each  side  of  the  source-A  peak  were  used  as  the  limits, 
yielding  values  of  188°-288°  for  the  Io-related  region,  204°-296°  for 
the  non-Io  part,  and  192°-292°  for  the  overall  mixture.  It  may  seem 
wiser  to  define  separate  limits  for  each  apparition's  A-.j  histogram, 
since  the  minima  positions  vary  greatly  from  year  to  year.  This  procedure 
is  difficult,  especially  in  low  activity  periods,  since  the  exact  position 
of  the  minimum  dividing  source  A  from  a  neighboring  source  is  often 
hard  to  determine.  Nevertheless,  an  attempt  was  made  to  carry  out 
this  process.  The  results  were  not  significantly  different  from  those 
of  Figure  III-8. 

The  dominance  of  the  Io-independent  radiation  in  Figure  III -8  in 
the  overall  source-A  motion  if  obvious.  Some  periodic  displacement 
appears  in  the  Io-related  sector,  but  it  must  be  remembered  that  this 
region  still  contains  some  non-Io-related  radiation.  Note  that  the 
inactive  years  in  this  segment,  years  when  Io-dependent  radiation 
dominates,  have  nearly  constant  centroid  positions.  Register  and  Smith 
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concluded  that  the  Io-dependent  radiation  retained  a  constant  rotational 
period  while  a  variable-amplitude,  variable-period,  non-Io-related 
radiation  was  superimposed  on  it.  Although  the  rotational  period  of 
the  Io-independent  source  A  may  indeed  be  variable,  the  oscillation  of 
the  centroid  can  be  explained  by  causes  other  than  a  complete  shift 
of  the  source  in  A....  coordinates.  This  will  be  discussed  again  in 
Chapter  V. 

As  mentioned  in  Chapter  I,  Carr  explained  the  sinusoidal  source 
motion  as  a  beaming  effect  related  to  the  11.86-year  cycle  of  the 
declination  of  the  earth  (Dp)  as  viewed  from  Jupiter,  resulting  mainly 
from  the  3°  difference  between  the  Jovian  orbital  and  rotational  axes. 
Donivan  and  Carr  (1969)  graphically  displayed  the  correlation  of  centroid 
position  with  Dp  in  the  manner  of  Figure  1 1 1-9.  This  figure  continues 
the  method  of  Donivan  and  Carr  up  to  1975,  using  the  data  of  Figure  1 1 1-8. 
The  centroid  for  each  apparition  is  plotted  according  to  the  mean  Dp 
value  obtained  by  averaging  declination  values  for  three  months  before 
and  after  opposition.  Weighted,  least-squares  lines  fitted  to  the  data 
measure  the  amount  of  correlation.  The  equations  for  the  non-Io-A, 
Io-A,  and  the  overall-A  lines,  respectively,  are  as  follows: 

A   (65)  -  238.9°  ±   .8°  +  (1.8+.2)DE 

Am(65)  =  251.0°  +  .8°  +  (2.7±.2)D£ 

A  n(65)  =  244.1°  i  .5°  +  (2.7+.l)TL 

where  the  errors  are  standard  deviations.  All  three  lines  fit  the  data 
quite  well,  but  the  larger  non-Io-A  slope  and  the  matching  overall-A 
slope  indicate  the  dominance  of  the  non-Io  radiation  in  the  periodic  shift. 

Donivan  and  Carr  also  looked  for  a  drift  in  the  source-B  18-MHz 
radiation,  finding  a  possible  displacement  in  the  same  phase  as  that  of 
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source  A.  With  the  larger  data  base  in  the  present  study  no  significant 
motion  was  found  for  either  source  B  or  source  C  (see  Figure  111-10). 
The  boundaries  for  source  B  were  defined  to  be  62°  <  A    <  188°,  and 
300°  <  Ain  <  30°  for  source  C.  Perhaps,  as  in  source  A,  the  variation 
occurs  in  the  non-Io-related  segment.  Unfortunately,  Io-independent-B 
storms  are  scarce  even  in  the  larger  data  base.  Non-Io-C  storms,  though 
more  numerous,  are  difficult  to  separate  from  their  non-Io-A  counterparts, 
and  the  position  of  this  boundary  is  critical  in  the  calculations.  The 
apparition  centroids  of  these  regions  were  computed,  however,  and  some 
minor  trends  could  be  suggested,  but  overall  the  results  were  too  erratic 
and  error-prone  to  be  conclusive. 

3.3.2  Io  Phase  Drifts 

In  the  process  of  generating  System  III  longitude  histograms  for 
each  apparition,  Io-phase  histograms  of  the  individual  source  regions 
were  also  generated,  as  shown  in  Figures  III-ll,  111-12,  and  111-13. 
The  method  of  isolating  the  sources  was  the  same  as  used  in  Figures  III-3, 
III-4,  and  III-5.  It  should  be  mentioned  that  the  earliest  apparition 
shows  artificially  large  probabilities  and  a  rugged  appearance  as  a 
result  of  insufficient  observation  time,  since,  by  chance,  the  small 
number  of  observing  periods  in  each  zone  may  each  have  contributed  a  storm. 
Close  examination  of  the  figures  reveals  that  the  prominent,  Io-dependent 
peaks  are  unstable  in  position.  This  could  result  from  random  fluctuations 
in  data  acquisition  and  combination,  but  a  plot  of  the  centroids  of  these 
peaks  suggest  differently.  Figure  1 1 1-14  compares  the  changing  centroid 
values  with  the  now-familiar  A    shift  of  source  A.  The  Io-phase 
motions  of  sources  B  and  C  seem  to  be  correlated  with  each  other  and 
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Figure  III-ll.  Individual-apparition  histograms  of  occurrence 

probability  for  2°  increments  in  the  phase  of  Io.  All 
18,  20,  and  22-MHz  storms  occurring  in  the  source-B 
System  III  longitude  (65)  region  are  plotted. 
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Figure  II 1-12.      Individual-apparition  histograms  of  occurrence 

probability  for  2°  increments   in  the  phase  of  Io.     All 
18,   20,  and  22-MHz  storms  occurring  in  the  source-A 
System  III   longitude  (65)  region  are  plotted. 
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Figure  111-13.  Individual-apparition  histograms  of  occurrence  probability 

for  2°  increments  in  the  phase  of  Io.  All  18,  20,  and  22-MHz 
storms  occurring  in  the  source-C  System  III  longitude  (65) 
region  are  plotted. 
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anticorrelated  with  the  A  .  shift  of  source  A,  that  is,  an  increase 

in  the  A  .,  position  of  A  corresponds  with  a  decrease  in  the  Io  phase 

of  B  and  C  and  vice  versa.  Io-related  A  exhibits  a  possible  Io-phase 

drift  in  the  early  apparitions,  where  larger  errors  are  apparent,  but 

later  years  show  no  significant  trends.  The  strange  appearance  of  the 

source-C  1957-58  and  1958-59  centroids  results  from  the  effects  of  very 

small  amounts  of  data  on  the  method  of  computing  errors. 

!   It  is  interesting  that  sources  B  and  C  appear  to  oscillate  in 

Io-phase  coordinates  but  show  no  consistent  motion  in  System  III 

longitude,  whereas  source  A  behaves  in  the  opposite  manner.  One  obvious 

difference  between  sources  B  and  C  and  source  A  is  the  ratio  of  non-Io- 

related  to  Io-related  radiation  involved  in  the  centroid  calculations. 

Since  the  data  used  to  compute  the  Io-phase  B  and  C  centroids  were 

confined  to  the  regions  65°  <  ATTT  <  200°,  48°  <  <t>T  <  126°  for  source  B, 
^  111  Io 

and  291°  <  ATTT  <  64°,  186°  <  <h  <  264°  for  source  C,  the  Io-related 
III  Io 

storms  should  greatly  outweigh  any  satellite-independent  radiation.  The 
Io-dependent  source-A  region  (201°  <  X...  <  290°,  184°  <  <K  <  286°), 
which  contains  a  larger  fraction  of  Io-independent  storms,  has  no 
oscillatory  motion,  although  the  larger  variability  of  the  activity  in 
this  region  makes  centroid  computation  difficult.  When  the  A,.,  centroids 
were  calculated  without  limits  on  the  phase  of  Io,  source  A  showed  the 
only  significant  motion  and  had  by  far  the  largest  ratio  of  Io-independent 
to  Io-dependent  radiation.  Further  discussion  of  the  motions  will  be 
delayed  until  after  they  have  been  studied  in  two  dimensions. 

The  correlation  of  the  Io-phase  drifts  with  the  averaged  declination 
of  the  earth  is  seen  in  Figure  111-15.  The  weighted,  least-squares 
analysis  yielded  the  following  equations  for  sources  B,  A,  and  C 


89 


4->  4-> 

■a      s_ 

c     ra  o 

re     <d  +-> 


"a 
x> 


-M      I    >* 


<c  -a  a> 


o 

s_ 

•  4- 

C 

o  c: 
•r-   cu 

•r-  (0 
S-  4-> 
fO 

a.  oj 

CL   S- 

ra    ra 


90 


respectively: 

♦Io  =  91.1°  ±  0.5°  +  (-l.liO.UlL 
*Io  =  236.9°  f  2.4°  +  (-0.4±0.5)DE 
♦Io  =  241.1°  +  1.0°  +  (-2.9±0.2)D 

The  slope  of  the  middle  equation,  which  is  within  a  standard  deviation 
of  being  zero,  corroborates  the  lack  of  motion  of  source  A.  The  other 
two  sources  show  good  correlation  with  D£,   although  the  smaller  drift 
in  source  B  results  in  a  slope  of  less  than  half  the  value  of  that  of 
source  C.  It  is  interesting  that  the  fitted  values  and  standard  deviations 
for  source  C  are  in  the  same  range  as  those  previously  computed  for  the 
Xjjj  movement  of  source  A. 

3.3.3  Asymmetries  and  Secondary  Peaks 

The  individual -apparition  histograms  were  also  useful  for  a 
year-by-year  study  of  Leacock's  discovery  (1971)  of  a  bifurcation  of 
the  source-B  peak  in  the  Io-phase  histograms.  It  is  easier  to  see  these 
secondary  peaks  in  Figure  111-16,  which  is  the  smoothed,  pseudo-three- 
dimensional  counterpart  of  Figure  III-ll.  All  source-B  peaks  except 
those  in  1957-58  and  1970-71,  which  are  low-data  years,  are  preceded  in 
Io  phase  by  one  or  more  secondary  peaks.  Quantitative  examination  of 
the  minor  maxima  is  difficult  because  of  their  low  storm  levels,  but 
general  trends  can  be  distinguished.  The  second-order  peaks  seem  to  be 
statistically  significant  since  they  are  often  more  than  2a  above  their 
surrounding  minima  as  measured  by  Equation  III-6.  Close  inspection 
reveals  that  they  do  not  retain  a  fixed  Io-phase  position,  but  move  in 
approximate  synchronization  with  the  major  peak  motion.  In  the  active 
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years  the  separation  of  the  most  prominent  secondary  peak  from  the 
main  peak  remains  in  the  range  of  20°-25°.  Secondary  peaks  following 
the  highest  peak  can  sometimes  be  seen  in  Figure  111-16  and  more 
readily  in  Figure  III-ll,  near  an  Io  phase  of  115°.  This  is  a  con- 
tribution of  the  18-MHz  data  in  the  mixture,  and  it  is  characteristic 
of  the  frequencies  below  20  MHz,  as  is  shown  in  Figure  111-17,  a 
different  view  of  Figure  III-3.  The  following  peak  is  apparent  at  18, 
16,  and  15  MHz,  but  it  diminishes  at  10  MHz.  Above  18  MHz,  only  pre- 
ceding secondary  maxima  seem  to  be  permitted.  Even  where  these  are 
missing,  a  more  persistent  characteristic  is  noticeable.  The  main  peak 
is  nearly  always  skewed  toward  smaller  values  of  Io  phase,  that  is,  the 
occurrence  probabilities  increase  slowly  with  Io  phase,  reach  a  maximum 
and  then  decrease  more  rapidly.  This  is  difficult  to  see  in 
Figure  111-17,  but  Figures  III-ll  and  II 1-16  provide  unobstructed 
viewing  angles. 

The  slow  onset  and  rapid  cutoff  of  Io-related  occurrence  probability 
is  manifest  in  more  than  just  the  B  region.  Figures  111-18  and  111-19 
taken  from  Figures  111-12  and  111-13,  reveal  similar  asymmetries  for  the 
A  and  C  regions.  The  non-Io-related  radiation  makes  the  source-A 
picture  somewhat  confusing,  but  the  asymmetry  is  still  evident.  Also 
reminiscent  of  source  B  are  the  secondary  peaks  which  often  occur  in 
the  rising  incidence  of  storms  preceding  the  A  and  C  peaks  in  Io  phase. 
Source  A  may  have  several  of  these,  while  there  is  normally  only  a 
single  bifurcation  of  source  C.  They  have  no  preferred  position, 
although  the  enhancement  in  the  C  region  occurs  fairly  often  in  the 
200°-210°  range.  Secondary  peaks  and  asymmetries  can  be  observed  at 
other  frequencies  above  10  MHz  in  the  plots  of  Figures  1 1 1-4  and  III-5. 
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Since  all  apparitions  have  been  averaged  for  each  of  these  frequencies, 
the  asymmetries  and  secondary  maxima  which  are  present  have  some 
stability  through  the  years. 

Skewness  is  a  quality  inherent  also  in  the  System  III  longitude 
picture  of  the  sources.  The  slow,  almost  linear  increase  of  source  B 
with  increasing  \„,   coordinates,  followed  by  a  more  rapid  descent,  is 
characteristic  above  10  MHz  (see  Figure  III-l).  Source  C  behaves  in 
just  the  opposite  manner,  so  that  sources  B  and  C  seem  to  bracket  and 
build  toward  source  A  until  an  interruption  gives  them  their  separate 
identities.  Source  A  is  usually  almost  symmetrical,  with  the  direction 
of  skew,  if  any,  changing  from  one  apparition  to  the  next  according  to 
the  rhythm  of  the  centroid  motion,  as  seen  in  Figure  III-7. 

Thus  far,  the  source  frequency  behavior,  centroid  motions, 
asymmetries,  and  bifurcations  have  been  treated  only  in  their  relation 
to  System  III  longitude  or  the  phase  of  Io.  Many  of  these  phenomena 
indicate  simultaneous  longitude  and  satellite  influence.  The  next 
chapter  will  study  the  second-order  effects  in  their  dependence  upon 
both  Ajjj  and  $.  ,  and  reveal  other  results  discovered  inthe  process. 


CHAPTER  IV 
TWO-DIMENSIONAL  STUDIES 

4 . 1  Projection  Methods  and  Problems 

The  best  means  of  examining  the  simultaneous  dependence  of  the 
decametric  radiation  on  two  independent  variables  is  a  three-dimensional 
graph,  which  is  difficult  to  project  onto  a  two-dimensional  page.  His- 
torically contour  plots,  like  Figure  1-2,  have  been  the  choice  in  the 
research  of  many  authors  (i.e.,  Alexander  et  al . ,  1974;  Wilson  et  al . , 
1968),  including  the  original  work  by  Bigg  (1964).  Three  difficulties 
limit  the  usefulness  of  these  plots:  the  resolution  of  the  underlying 
data  is  hard  to  ascertain;  fine  structure  is  often  lost  between  contour 
lines  unless  a  confusingly  large  number  of  lines  is  used;  and  the  sense 
of  relief  is  unsatisfactory  for  the  average  observer.  Another  mode  of 
presentation  is  the  hidden-line,  perspective  view  like  Figure  III— 1 . 
This  certainly  improves  comprehension  of  relief,  but  the  hidden-line 
technique  can  conceal  important  detail,  and  deciding  exact  coordinate 
locations  of  a  particular  feature  is  troublesome. 

A  third  method  has  been  tried  which  emphasizes  the  two-dimensional 
histogram  basis  for  studying  the  data.  In  this  process,  the  two 
independent  variables  are  divided  into  intervals,  which  in  turn  segments 
the  coordinate  plane  into  rectangular  "bins."  Since  the  number  of  bins 
increases  geometrically  when  both  variables  are  subdivided,  compactness 
of  display  of  the  results  in  each  two-dimensional  zone  is  helpful. 


97 


Assignment  of  a  single  line-printer  character  for  each  bin  was  the 
least  expensive  and  most  convenient  solution.  Figure  IV-1  demonstrates 
the  technique  for  an  occurrence  probability  histogram  of  18-MHz  data 
from  the  years  1963-65.  In  analogy  with  digital  imaging  systems  each 
of  the  characters  could  be  considered  as  a  code  for  the  "brightness"  of 
a  picture  element  or  pixel,  but  the  lack  of  imaging  equipment  makes  it 
necessary  to  display  only  the  symbols  themselves.  The  interval  sizes  used 
for  the  System  III  longitude  and  the  phase  of  Io  axes  were  5°  each.  The 
print  characters  are  valued  according  to  the  order  shown  at  the  bottom  of 
the  picture.  A  total  of  52  characters  are  used.  A  blank  means  no 
activity  occurred  in  that  bin,  and  a  zero  represents  the  maximum 
occurrence  probability  zone.  In  this  case  the  small  amount  of  data  resulted 
in  several  zones  having  the  maximum  probability  of  1.0.  The  other 
50  characters  symbolize  equally-divided  percentiles  of  the  maximum 
probability,  beginning  with  the  "?"  representing  0-2  per  cent,  and 
ending  with  the  "1"  representing  98-99.9  per  cent. 

This  method  displays  nearly  all  of  the  detail  attainable  from  the 
data,  but  again  there  are  problens.  Comprehension  of  relief  is  even 
more  difficult  than  in  the  contour  plots.  Only  the  contrast  between  the 
characters  and  the  blank  areas  provides  any  immediate  sense  of  the  third 
dimension.  Many  times,  the  data  are  so  diffuse  that  few  blank  areas 
appear.  Plots  of  this  nature  can  be  improved  by  eliminating  the. less- 
important  lower  levels  of  probability  by  covering  them  with  low-lying 
"clouds."  This  has  been  done  in  Figure  IV-2,  in  which  all  characters 
below  the  10  per  cent  level  have  been  replaced  with  blanks.  The  observer 
still  has  trouble  distinguishing  the  shape  of  the  peaks  penetrating 
this  cloud  layer.  Apparently  no  single  plotting  method  is  free  of 
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Figure  IV-1.  Histogram  of  relative  occurrence  probability  for  5°  intervals 
in  System  III  Longitude  (65)  and  the  phase  of  Io,  based  on 
18-MHz  data  from  the  apparitions  of  1953-65.  The  characters 
represent  percentiles  of  the  maximum  occurrence  probability. 
Their  number  and  rank  from  lowest  to  highest  are  listed  at 
the  bottom. 
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Figure  IV-2.  Histogram  of  relative  occurrence  probability  for  5°  intervals 
in  System  III  longitude  (65)  and  the  phase  of  Io,  based  on 
18-MHz  data  from  the  apparitions  of  1963-65.  The  characters 
represent  percentiles  of  the  maximum  occurrence  probability. 
Their  number  and  rank  from  lowest  to  highest  are  listed  at  the 
bottom.  All  probabilities  below  0.100  have  been  replaced  by 
blanks. 
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defects.     Consequently,  a  combination  of  two  techniques  will  generally 
be  used  to  illustrate  the  results  in  the  next  section.     The  line- 
printer-character  display  provides  the  resolution  needed  for  observation 
of  second-order  effects,   and  a  perspective  plot  of  the  same  data  should 
provide  a  perception  of  the  third  dimension. 

Figures  IV-1   and  I V- 2  also   illustrate  a  problem  inherent  in  the 
data  which  all  methods  of  displaying  two-dimensional   histograms  inherit. 
In  low-storm-level    regions  of  the  plot  linear,   downward-sloping  structure 
often  stands  out.      For  illustration,  see  the  arrangement  of  characters  in 
Figure  IV-1   beginning  at  an  Io  phase  around  300°  and  extending  from 
0°-100°  in  System  III   longitude.     This  type  of  feature  represents  the 
contribution  of  a  single  lengthy  radio  storm  or  the  superposition  of 
several   storms.     The  slope  is   characteristic  of  the  family  of  lines  which 
are  described  by  the  elimination  of  the  time  parameter  from  the  parametric 
equations  of  the  time  evolution  of  Am  and  $  The  magnitude  of  the 

slope  is   .234,   the  ratio  of  the  rotational   period  of  Jupiter  to  the 
orbital   period  of  Io.     Thus,  any  two  storms  having  nearly  equal    UTTT,<j>T    ) 
coordinates  at  some  point  in  their  duration  will    lie  along,  the  same 
lines  or  bands  and  will   be  collectively  summed  where  their  tracks     in  the 
histogram  overlap.     The  downward-sloping  streaks  can  occur  in  low- 
probability  areas  of  even  the  largest  accumulations  of  data.     This  is  due 
to  a  selection  effect  caused  when  most  of  the  storms  in  these  regions  take 
place  during  the  active  years.     Only  one  complete  period  of  active  years 
has  been  observed,  and  the  incomplete  coverage  of  the  L..-L     plane 
prejudices  the  data.      It  is  expected  that  inclusion  of  other  active  periods 
will   gradually  smooth  out  the  ridges  and  obviate  this  effect.     The 
high-probability,   Io-related  source  regions  will   also  have  some  bias 
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toward  more  active  years,  but  in  their  case  the  inactive  years  contribute 
a  larger  number  of  storms  and,  therefore,  the  coverage  is  more  continuous. 
Care  must  be  taken  to  recognize  the  bias  toward  inclined  structure  so 
that  it  is  not  mistaken  for  phenomena  independent  of  a  low  probability 
level . 

4.2  Spectral  Analysis 

Chapter  III  began  the  study  of  source  metamorphosis  through  the 
frequencies  by  simultaneous  display  of  the  one-dimensional  histograms 
for  heavily-observed  wavelengths.  Since  simultaneous  comparison  of 
two-dimensional  histograms  for  all  frequencies  is  not  possible  without 
a  fourth  dimension,  each  frequency  will  be  treated  in  turn  starting  at 
the  high  end  and  working  down. 

Distribution  of  the  data  into  two-dimensional  zones  stretches  even 
the  largest  accumulation  of  data  to  the  breaking  point.  Consequently, 
the  frequencies  of  5.6  and  25.2  MHz  have  been  dropped  from  consideration 
in  the  present  chapter  for  insufficient  data.  All  of  the  other  frequencies 
have  over  200,000  minutes  of  observing  time.  Even  at  these  levels  the 
histograms  for  the  lesser-observed  frequencies  cannot  be  meaningfully 
subdivided  beyond  5°  intervals  for  the  two  independent  variables.  For 
accuracy  of  occurrence  probabilities,  the  calculations  have  been  based 
on  listening  and  activity  times,  although  the  number  of  events  in  the  zones 
is  still  monitored.  The  uniformity  of  the  data  can  still  be  measured 
from  the  average  and  standard  deviation  of  the  collection  of  the  listening 
times  in  the  zones. 

Three  figures,  IV-3  to  IV-5,  display  the  30-MHz  results.  This  end 
of  the  decametric  spectrum  has  been  chosen  as  a  starting  point  since 
future  diagrams  will  become  more  complicated  as  frequency  decreases. 
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Figure  IV-3.  Histogram  of  relative  occurrence  probability  for  5°  intervals 
in  System  III  longitude  (65)  and  the  phase  of  Io,  based  on  all 
30-MHz  data.  The  characters  represent  percentiles  of  the 
maximum  occurrence  probability.  Their  number  and  rank  from 
lowest  to  hiqhest  are  listed  at  the  bottom.  The  arrow  indicates 
the  direction  of  view  for  Figure  IV-4. 
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Figures  IV-3  and  IV-4  are  the  detailed  presentation  and  perspective 

view  of  the  Am  -  ^  histogram.  Zones  of  5°  by  5°  were  chosen  since 

division  of  the  210,000  minutes  of  observing  time  into  5184  parts  gives 

an  average  of  40  minutes  and  6  observations  per  zone,  numbers  considered 

too  small  for  further  subdivision.  Percentiles  of  relative  occurrence 

probability  are  represented  by  the  same  character  sequence  as  in 

Figures  IV-1  and  IV-2.  The  sequence  will  remain  the  same  in  all  of 

the  similar  plots  in  this  chapter.  The  arrow  in  Figure  IV-3  indicates 

the  direction  of  observation  used  in  Figure  IV-4.  This  angle  was  chosen 

as  a  standard  for  all  of  the  perspective  plots,  since  it  is  the  most 

beneficial  for  revealing  details  to  be  discussed  later. 

The  outstanding  and  almost  the  only  feature  in  either  figure  is 

the  Io-related,  source-B  radiation.  Though  surrounded  by  emptiness, 

the  source-B  island  still  exhibits  a  high  degree  of  probability,  since 

the  maximum  occurrence  probability  is  0.95.  Its  shape  in  the  ATTT  -  A 

III   *lo 

plane  of  Figure  IV-3  is  typical.  There  is  an  elongated  smear  of  data 
from  75°  to  180°  in  \m   and  from  85°  to  110°  in  <frI(j.  This  is  augmented 
by  a  projection  extending  to  the  upper  left.  While  the  upper  edge  if 
irregular  in  shape,  the  lower  edge  cuts  off  sharply  at  <(..  of  110°,  which 
explains  the  slow  rise  and  sharp  cutoff  of  source  B  in  the  one-dimensional 
histogram. 

In  the  traditional  locations  for  sources  A  and  C  (<j>.  =  240°  and 
AIII  =  240°  and  330°'  respectively)  there  are  slight  indications  of  storm 
activity.  These  are  not  significantly  different  from  other  storms 
scattered  elsewhere.  The  dispersion  at  first  appears  random,  but  no 
storms  occur  in  what  is  called  the  "null  zone"  to  the  left  side  of 
source  B  in  Figure  IV-3.  In  this  Am  range,  no  Io-related  source  appears 
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until  source  D  emerges  around  16  MHz.  Since  the  other  Io-related  sources 
persist  to  much  higher  frequencies,  this  suggests  that  the  scattered 
radiation  may  be  harbingers  of  these  sources,  while  storms  preceding 
source  D's  appearance  will  not  appear  until  lower  frequencies  are 
encountered.  The  small  number  of  outlying  storms  makes  this  statement 
rather  speculative. 

The  graphs  in  Figure  IV-5  are  2°-interval  histograms  for  System  III 
longitude  and  for  the  phase  of  Io  in  isolated  source  B,  A,  and  C  regions. 
Although  these  plots  did  not  specifically  appear  in  the  previous  chapter, 
they  served  as  the  basis  for  the  30-MHz  parts  of  Figures  III-l  and 
III-3  to  III-5.  Figure  IV-5  and  similar  figures  in  this  chapter  will  serve 
to  facilitate  comparison  of  features  between  the  one-  and  two-dimensional 
histograms.  In  the  30-MHz  case  the  source-B  peaks  are  obvious  in  the 
upper  two  histograms  and  the  lack  of  concentrated  storms  in  the  isolated 
A  and  C  areas  makes  the  lower  histograms  seem  erratic.  The  System  III 
longitude  histogram  reinforces  the  identification  of  the  scattered 
radiation  as  precursors  of  sources  A  and  C,  since  secondary  peaks  appear 
at  the  appropriate  longitudes. 

Between  30  and  27  MHz  Io-related  source  A  appears,  as  is  shown  in 
Figures  IV-6  and  IV-7.  Figure  IV-6  has  more  detail  than  its  30-MHz 
counterpart  since  the  larger  amount  of  observing  and  activity  time  at 
27  MHz  permits  division  of  the  data  into  2°  by  2°  zones.  Even  with  the 
smaller  intervals,  there  is  still  an  average  of  49  minutes  of  observation 
and  18  separate  listening  periods  for  each  zone.  The  individual 
characters  are  difficult  to  see,  but  for  cursory  analysis  the  more 
important  feature  is  the  relationship  between  the  blank  and  lettered  areas. 


Figure  IV-6.  Histogram  of  relative  occurrence  probability  for  2° 

intervals  in  System  III  longitude  (65)  and  the  phase  of 
Io,  based  on  all  27-MHz  data.  The  characters  represent 
percentiles  of  the  maximum  occurrence  probability.  Their 
number  and  rank  from  lowest  to  highest  are  listed  at  the 
bottom. 


109 


o 

roO 


*0 


60 


|°  120         150  180        210         240        270        300 

SYSTEM    m    LONGITUDE     (65) 


330        360 


no 


UJ  >- 
o  LU  -J 

-i  tr  <, 

3  03 

o  o 

o  CC 
o  a 


•i-  CD 

cn>— 

o  o 
1 —  (/) 


E  -c 

CD  -t-> 
4->  -r- 

1/5  CD 

to 

S-  A3 

o  -a 

4- 


.—  CD 

■I-  > 

-Q  O 


1 

al 

r-~. 

Q 

>, 

Cvl 

F 

+J 

CI 

4- 

C   ) 

o 

XI 

« 

«3 

CD 

O  XI 

o 

> 

5-                1 

4- 

CD- 

Q> 

O 

> 

CD 

CD  -SZ 

+-> 

in 

+-> 

O 

01 

CD  -C 

</) 

O- 

o--o 

1/1 

<- 

a; 

CD 

CD 

-c 

Q- 

4-> 

>■> 

Ill 


Source  B  has  changed  slightly  from  its  30-MHz  appearance.  Not  only 
is  there  a  projection  to  the  upper  left,  but  a  low-probability  area  has 
appeared  at  the  upper  right.  The  combination  of  these  two  regions  explains 
the  bifurcated  peak  at  70°   in  the  Io-phase  histogram  of   source  B  in 
Figure  IV-8.  Although  source  A  has  just  appeared  and  is  fairly  weak, 
the  corresponding  Io-phase  histogram  already  exhibits  the  secondary  peaks 
and  asymmetry  discussed  in  the  previous  chapter.  The  slow 'rise  of  source 
A  with  increasing  Io  phase  corresponds  to  the  slope  leading  up  to  the 
peak  in  Figure  IV-7,  and  it  extends  over  the  entire  30°  longitudinal  width 
of  the  source.  The  maximum  probability  is  not  centered  in  this  range,  but 
occurs  at  a  slightly  lower  longitude.  It  is  difficult  to  say  whether  the 
secondary  peaks  come  from  bumps  in  the  slope  or  ridges  extending  across 
the  width.  Assuming  that  the  cause  for  the  secondary  peaks  is  the  same  at 
lower  frequencies,  their  origins  will  be  easier  to  determine  in  the  next 
set  of  figures. 

The  most-observed  frequency  is  22  MHz,  with  an  average  of  80  minutes 
and  30  observations  per  2°  by  2°  zone.  This  seems  poetically  just,  since 
the  original  discovery  of  Jupiter's  radio  emission  was  made  at  this 
frequency.  The  increase  in  monitoring  time  from  27  to  22  MHz  is  also 
justified  by  the  larger  amount  of  activity,  displayed  in  Figures  IV-9  and 
IV-10.  To  unclutter  Figure  IV-9  (and  other  similar  figures  at  lower 
frequencies)  all  zones  having  probabilities  less  than  the  average  occurrence 
probability,  defined  as  the  total  activity  time  divided  by  the  total 
listening  time,  have  been  blanked  out.  In  this  case  the  cutoff  value  is 
0.041. 

Between  27  and  22  MHz,  Io-related  source  C  has  emerged  as  well  as 
profuse  Io-independent  radiation  in  the  A  and  C  regions.  The  non-Io-related 
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Figure  IV-9.  Histogram  of  relative  occurrence  probability  for  2° 

intervals  in  System  III  longitude  (65)  and  the  phase  of 
Io,  based  on  all  22-MHz  data.  The  characters  represent 
percentiles  of  the  maximum  occurrence  probability.  Their 
number  and  rank  from  lowest  to  highest  are  listed  at  the 
bottom.  All  probabilities  below  0.041  have  been  replaced 
by  blanks. 
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radiation  shows  the  downward-sloping  structure  characteristic  of  low- 
probability  regions.  The  Io-phase  secondary  peaks  preceding  sources  A 
and  C  in  Figure  IV-11  coincide  with  ridges  usually  extending  across  the 
entire  longitudinal  width  of  the  source.  Whether  they  are  artifacts  of 
the  active  years  or  persistent  through  the  other  periods  will  be  studied 
in  the  chronology  section  of  this  chapter. 

As  is  true  at  27  MHz,  the  bifurcation  and  asymmetry  of  the  Io-phase, 
source-B  peak  in  Figure  IV-11  are  products  of  the  projections  extending 
to  the  upper  right  and  upper  left  in  Figure  IV-9.  Unlike  the  situation 
at  the  higher  frequency,  the  main  body  of  the  source  extends  to  lower 
values  of  Io  phase,  but  the  sharp  cutoff  at  the  higher  Io-phase  edge 
remains  at  the  same  position.  The  shape  of  Io-related  source  C  should 
also  be  mentioned.  In  this  case  an  interesting  bulge  extends  to  lower 
values  of  Io  phase,  ending  around  220°. 

There  is  not  too  much  change  between  the  22-MHz  picture  and  the 
20-MHz  data  presented  in  Figures  IV-12  to  IV-14.  The  increased  probability 
of  the  non-Io-related  radiation  from  sources  A  and  C  makes  observation  of 
Io-related  source  C  in  Figure  IV-13  more  confusing,  but  its  growth  is  still 
evident  from  the  higher  probabilities  and  more  prominent  extension  past 
360°  and  into  the  low  longitudes.  The  increased  activity  has  also  raised 
the  average  occurrence  probability  cutoff  to  0.059,  but  the  number  of 
observations  and  listening  time  per  zone  have  dropped  to  16  and  41  minutes 
respectively.  The  projection  to  the  upper  left  of  source  B  in  Figure  IV-12 
has  increased  in  size,  while  the  indentation  between  the  two  bulges  has 
been  partially  filled  with  a  small  number  of  storms.  This  helps  to 
explain  the  lack  of  bifurcation  in  Io-phase  of  source  B  in  Figure  IV-14, 
but  conspicuous  secondary  peaks  are  absent  from  the  other  sources  as  well, 
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Figure  IV-12.  Histogram  of  relative  occurrence  probability  for  2° 
intervals  in  System  III  longitude  (65)  and  the  phase 
of  Io,  based  on  all  20-MHz  data.  The  characters 
represent  percentiles  of  the  maximum  occurrence 
probability.  Their  number  and  rank  from  lowest  to 
highest  are  listed  at  the  bottom.  All  probabilities 
below  0.059  have  been  replaced  by  blanks. 
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although  the  asymmetry  is  still  unmistakeable.  This  argues  for  the 
random  superposition  origin  of  the  secondary  maxima. 

The  18-MHz  data  have  a  still  higher  level  of  activity,  elevating 
the  average  occurrence  probability  to  0.065,  and  it  has  been  monitored 
slightly  more  than  its  20-MHz  predecessor  with  an  average  of  20  observations 
and  53  minutes  of  listening  in  each  2°-square  zone.  Previous  figures  have 
shown  some  signs  of  non-Io-related  radiation  in  the  source-'B  longitudes, 
but  Figure  IV-15  has  a  consistent  low  level  of  radiation  scattered  in 
the  120°-180°  longitude  region.  There  is  a  definite  lack  of  Io-independent 
radiation  in  the  70°-120°  region,  even  though  the  Io-related  part  of 
source  B  extends  over  to  the  70°  longitude.  One  major  change  has  taken 
place  between  20  and  18  MHz.  The  Io-related  part  of  B  no  longer  cuts 
off  sharply  at  an  Io  phase  of  approximately  100°  but  continues  on  to  120°. 
As  a  result  the  Io-phase,  source-B  peak  in  Figure  IV-17  is  now  more 
symmetric  and  trifurcated.  The  source  A  and  C  regions  remain  virtually 
the  same  as  before.  Source  C  has  increased  in  probability  and  extended 
farther  into  the  low-longitude  region.  Some  secondary  peaks  have  reappeared 
in  Figure  IV-17  preceding  Io-related  source  A,  but  the  surrounding 
fluctuation  of  Io-independent  radiation  makes  them  undistinguished. 

The  shorter  total  monitoring  time  at  16  MHz  necessitates  returning 
to  5°  by  5°  zones  for  the  histograms  of  Figures  IV-18  and  IV-19.  The  time 
has  been  divided  into  an  average  of  68  minutes  per  zone  contributed  by 
10  observations.  The  average  occurrence  probability  cutoff  for 
Figure  IV-18  is  0.067.  Figures  IV-18  to  IV-20  do  not  reveal  much  change 
in  sources  A  and  C,  although  the  non-Io-A  radiation  has  higher  probabilities 
and  a  more  erratic  appearance  than  at  18  MHz.  This  can  be  blamed  on  the 
lower  data  levels.  The  alteration  of  source  B  is  more  noticeable.  The 


Figure  IV-15.  Histogram  of  relative  occurrence  probability  for  2° 

intervals  in  System  III  longitude  (65)  and  the  phase  of 
Io,  based  on  all  18-MHz  data.  The  characters  represent 
percentiles  of  the  maximum  occurrence  probability.  Their 
number  and  rank  from  lowest  to  highest  are  listed  at  the 
bottom.  All  probabilities  below  0.065  have  been  replaced 
by  blanks. 
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Figure  IV-18.  Histogram  of  relative  occurrence  probability  for  5°  intervals 
in  System  III  longitude  (65)  and  the  phase  of  Io,  based  on 
all  16-MHz  data.  The  characters  represent  percentiles  of 
the  maximum  occurrence  probability.  Their  number  and  rank 
from  lowest  to  highest  are  listed  at  the  bottom.  All 
probabilities  below  0.067  have  been  replaced  by  blanks. 
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upper  left  bulge,  characteristic  of  the  source  at  higher  frequencies, 
has  all  but  disappeared  in  Figure  IV-18,  and  the  lower  edge  of  the  source 
extends  to  an  Io  phase  of  120°  as  it  did  at  18  MHz.  The  Io-related 
source-B  peak  in  Figure  IV-20  is  now  skewed  in  the  opposite  manner  from 
its  behavior  at  frequencies  above  18  MHz.  Also,  the  non-Io-related 
source-B  radiation  has  resumed  an  intermittent  nature  similar  to  that  at 
frequencies  above  18  MHz. 

The  amount  of  data  increases  at  15  MHz  to  an  average  of  94  minutes 
and  14  observations  per  5°  by  5°  zone.  The  average  occurrence 
probability  has  also  risen  slightly,  to  0.075.  Two  major  differences 
make  Figure  IV-21  distinct  from  Figure  IV-18.  Between  16  and  15  MHz, 
source  D  has  become  visible  between  the  longitudes  of  15°  and  100°  at  an 
Io  phase  above  100°.  The  source  has  a  surprisingly  high  degree  of  proba- 
bility, considering  that  it  was  not  noticeable  at  a  frequency  one  MHz 
higher.  There  also  seems  to  be  an  arching  bridge  of  radiation  between 
the  edge  of  source  C  at  60°  longitude  and  the  200°  edge  of  source  A. 
Figure  IV-23  exhibits  the  same  skewing  of  source  B  toward  higher  values 
of  Io  phase  and  the  same  bifurcation  that  was  present  at  16  MHz.  Some 
interesting  secondary  peaks  are  present  in  the  non-Io-related  part  of  the 
Io-phase  source-B  histogram.  The  two  at  105°  and  240°  relate  to  source  D 
and  the  storms  linking  sources  A  and  C  through  the  B  region.  Other 
secondary  peaks  in  all  of  the  Io-phase  histograms  may  be  due  to  the 
superposition  of  storms  in  the  active  years,  although  the  peak  just 
preceding  source  C  in  Io  phase  seems  to  be  persistent  through  many 
frequencies. 

A  drastic  change  occurs  between  15  and  10  MHz.  Figure  IV-24  only 
vaguely  resembles  Figure  IV-21  because  non-Io-related  radiation  has 
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Figure  IV-21 .  Histogram  of  relative  occurrence  probability  for  5°  intervals 
in  System  III  longitude  (65)  and  the  phase  of  Io,  based  on  all 
15-MHz  data.  The  characters  represent  percentiles  of  the 
maximum  occurrence  probability.  Their  number  and  rank  from 
lowest  to  highest  are  listed  at  the  bottom.  All  probabilities 
below  0.075  have  been  replaced  by  blanks. 
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Figure  IV- 24.  Histogram  of  relative  occurrence  probability  for  5°  intervals 
in  System  III  longitude  (65)  and  the  phase  of  Io,  based  on 
all  10-MHz  data.  The  characters  represent  percentiles  of  the 
maximum  occurrence  probability.  Their  number  and  rank  from 
lowest  to  highest  are  listed  at  the  bottom.  All  probabilities 
below  0.121  have  been  replaced  by  blanks. 
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spread  throughout  the  *m-<i>Io  plane.  The  Io-related  source  regions 
no  longer  dominate  the  picture  with  higher  probabilities,  and  so  it  is 
difficult  to  pick  out  any  features  in  Figure  IV-25.  The  most  noticeable 
part  of  Figure  IV-24  is  the  downward-sloping  structure  indicative  of 
the  fact  that  this  frequency  is  one  of  the  lesser-observed,  with  an 
average  of  47  minutes  of  observation  accumulated  in  7  sessions  for  each 
of  the  5°  by  5°  zones.  Even  with  the  elimination  of  all  characters 
representing  values  below  the  average  occurrence  probability  of  0.121, 
the  picture  remains  cluttered. 

The  biggest  aid  in  sorting  out  the  features  is  Figure  IV-26. 
According  to  the  longitude  plot  in  this  figure,  the  probability  peaks  at 
210°  and  350°.  The  radiation  near  210°  is  fairly  evenly  spread  over 
all  values  of  Io  phase.  This  is  revealed  in  the  one-dimensional, 
source-A,  Io-phase  histogram  and  in  close  examination  of  Figure  IV-24. 
The  Io-related  part  of  source  A  in  Figure  IV-24  is  no  longer  detectable 
above  its  surroundings.  Io-B  is  still  noticeable,  but  the  dominant  part 
of  it  occurs  at  a  higher  longitude  than  in  the  15-MHz  case.  Upward 
projections  from  the  source  toward  lower  Io  phase  are  no  longer 
distinguishable  from  the  background  radiation,  if  they  are  even  present. 
It  is  difficult  to  separate  the  Io-related  parts  of  source  B  and  source  D, 
since  B  has  narrowed  in  its  Io-phase  extent.  A  gap  at  approximately  120° 
longitude  seems  to  separate  D  from  B.  If  this  is  true,  then  source  B 
begins  at  a  higher  longitude  and  has  its  maximum  at  higher  values  than 
in  the  upper-frequency  diagrams,  suggesting  that  the  whole  source  has  moved 
upward.  Also,  the  null  zone  of  Io-independent  radiation  has  changed  from 
the  50°-100°  longitudes  to  the  100°-170°  band.  Further  evidence  of 
migration  is  found  in  Kennedy's  polarization  work  (1969).  The  source  A 
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and  B  regions  are  known  to  be  polarized  in  the  right-handed  sense, 
whereas  sources  C  and  D  generally  have  the  opposite  sense,  though  they 
are  more  of  a  mixture.  Kennedy  finds  that  the  transitions  between 
dominance  by  one  sense  to  dominance  by  the  other  occur  at  approximately 
100°  and  330°  of  longitude  in  the  mid-decametric  frequencies.  At  10  MHz 
the  upper  boundary  remains  nearly  the  same,  but  the  lower  one  has  moved  to 
around  130°,  indicating  a  higher-longitude  beginning  for  source  B. 
Careful  comparison  of  Figures  IV-24  and  IV-21  suggests  that  between  15 
and  10  MHz  the  non-Io-related  source-A  storms  have  split  away  from  their 
source-C  counterparts  and  moved  to  lower  longitudes.  It  seems  reasonable 
to  say  that  the  longitude  probability  maxima  of  sources  A  and  B  have  moved 
toward  each  other  and  blended  together  to  form  the  longitude  histogram 
peak  at  210° as  implied  earlier  by  the  lines  in  Figure  III-2. 

The  longitude  peak  at  350°  in  Figure  IV-26  is  obviously  formed  from 
the  combination  of  Io-related  and  Io-dependent  source  C.  Secondary  peaks 
in  the  60°-100°  longitude  range  align  with  the  Io-related  source-D 
radiation  located  at  an  Io  phase  of  90°.  Another  secondary  peak  at  150° 
seems  to  be  caused  by  the  increased  number  of  storms  near  .an  Io  phase  of 
240°.  These  storms  are  in  the  same  position  as  the  "bridge"  at  15  MHz 
connecting  sources  C  and  A.  In  comparing  Figure  IV-24  to  the  15-MHz 
equivalent,  traces  of  a  bridge  can  still  be  seen.  Little  can  be  said 
about  the  asymmetry  of  the  Io-related  sources  still  present,  since  the 
pervasive  Io-independent  storms  obscure  the  fine  details. 

4.3  Io-Related  Source  Chronology 

Subdivision  of  the  two-dimensional  histograms  into  shorter  time 
periods  will  again  require  the  use  of  the  largest  possible  set  of  data. 
As  in  the  previous  chapter,  the  merged  18,  20,  and  22-MHz  data  will  be 
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used,  since  they  have  been  extensively  observed  and  their  histograms 
closely  resemble  each  other.  Figures  IV-27  to  IV-29  display  the 
properties  of  the  merged-frequency  data  for  all  years.  Some  individual 
characteristics  of  the  frequencies  are  still  apparent,  such  as  the 
radiation  on  the  higher  Io-phase  side  of  source  B  contributed  by  the 
18-MHz  data.  For  Figure  IV-27  the  2°  by  2°  zones  now  represent  an 
average  of  175  minutes  of  observation  and  65  observing  periods.  All 
characters  symbolizing  occurrence  probabilities  below  the  average  of 
0.053  have  been  replaced  by  blanks.  It  will  be  noticed  that  the  downward- 
slanting  structure  indicative  of  low  storm  levels  is  still  visible  in  the 
low-probability  regions.  Attempts  to  separate  the  data  in  these  areas  into 
different  time  periods  would  reduce  them  to  a  near-meaningless  collection 
of  slanted  lines.  Effort  has  been  concentrated  instead  on  the  high- 
probability,  Io-related  parts  of  the  sources.  Even  for  these  areas 
subdivision  of  the  data  into  single  apparitions  would  display  too  small 
a  number  of  storms  for  the  years  of  low  activity.  For  this  reason  the 
apparitions  during  the  times  of  negative  declination  of  the  earth,  D  , 
will  be  grouped  together. 

The  changes  in  source  B  through  the  years  can  be  seen  in  Figure  IV-30. 
Nearly  all  of  the  histograms  show  some  storm-evolution  streaks,  which  are 
difficult  to  eliminate  completely  from  short  spans  of  data.  The  same 
problem  causes  artificially  high  probability  concentrations  such  as  the 
large  number  of  l's  in  the  1962-63  diagram.  The  general  structure  of  the 
source  can  still  be  discerned  despite  these  problems.  The  process  of 
blanking  out  all  characters  below  the  average  occurrence  probability  for 
the  region  has  been  carried  out  for  each  diagram. 


Figure  IV-27.  Histogram  of  relative  occurrence  probability  for  2° 

intervals  in  System  III  longitude  (65)  and  the  phase  of 

Io,  based  on  combined  18,  20,  and  22-MHz  data.  The  characters 

represent  percentiles  of  the  maximum  occurrence  probability. 

Their  number  and  rank  from  lowest  to  highest  are  listed 

at  the  bottom.  All  probabilities  below  0.053  have  been 

replaced  by  blanks. 
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Figure  IV-30.  Comparison  of  histograms  of  relative  occurrence  probability 
for  2°  intervals  in  System  III  longitude  (65)  and  the  phase 
of  Io,  centered  on  the  Io-related  source-B  region.  Combined 
18,  20,  and  22-MHz  data  have  been  split  into  individual 
apparitions  for  years  of  positive  D£  and  groups  of 
apparitions  spanning  negative  Dr  years.  The  number  and 
rank  of  the  characters  remain  trie  same  as  in  Figure  IV-1 . 
All  probabilities  below  the  average  occurrence  probability 
for  each  diagram  have  been  replaced  by  blanks. 
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The  figure  suggests  that  the  motion  of  source  B  in  Io  phase  is 
caused  by  the  entire  source  shifting  its  Io  phase  boundaries  with  time, 
rather  than  having  an  oscillation  of  a  smaller  but  heavily-weighted 
segment  within  the  source.  For  example,  the  lower  edge  (side  toward 
higher  Io  phase)  remains  near  an  Io  phase  of  100°  in  the  active  years, 
but  in  the  combined  negative  D£   years  this  extends  to  between  105°  and 
110°.  The  more  interesting  part  is  the  opposite  edge.  The  two  upward 
projections  and  the  intervening  blank  area  vary  in  appearance  over  the 
years.  Some  of  the  change  can  be  attributed  to  the  incomplete  coverage 
and  low  storm  levels  of  this  region.  Overall  examination  shows  a  build-up 
of  the  upper  structure  in  the  active  part  of  the  cycle.  At  those  times 
the  projections  extend  to  around  55°  of  Io  phase  on  the  left  and  60°-65° 
on  the  right.  In  the  middle,  there  is  generally  an   indentation  to  around 
75°  at  a  longitude  near  125°.  The  secondary  peaks  leading  the  main  peak 
in  Io  phase,  as  discussed  in  the  previous  chapter,  are  formed  from  the 
combined  protuberances,  and  so  they  are  more  probable  during  the  active 
years.  The  years  of  lower  probability  seem  to  have  a  more  uniform  upper 
edge  with  projections  absent  and  the  indentation  filled.  One  other 
characteristic  that  can  be  seen  in  many  apparitions  may  be  important. 
The  high-longitude  side  of  the  source  is  usually  almost  linear  in  nature, 
sloping  downward  at  an  angle  which  differs  from  the  slant  characteristic 
of  a  storm's  continuance  in  time.  The  implications  of  all  of  these 
phenomena  will  be  discussed  further  in  Chapter  V. 

The  study  of  source  A  in  Figure  IV-31  has  a  heavily-streaked 
appearance.  Source  A  is  not  as  probable  as  source  B  and  it  is  bracketed 
by  even  lower-probability,  non-Io-related  radiation  on  the  top  and  bottom. 
This  makes  it  difficult  to  acquire  enough  data  in  a  limited  time  span  to 


Figure  IV-31.  Comparison  of  histograms  of  relative  occurrence 

probability  for  2°  intervals  in  System  III  longitude  (65) 
and  the  phase  of  Io,  centered  on  the  Io-related  source-A 
region.  Combined  18,  20,  and  22-MHz  data  have  been  split 
into  individual  apparitions  for  years  of  positive  Df  and 
groups  of  apparitions  spanning  negative  [)£   years.  The 
number  and  rank  of  the  characters  remain  the  same  as  in 
Figure  IV-1 .  All  probabilities  below  the  average 
occurrence  probability  for  each  diagram  have  been 
replaced  by  blanks. 
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fill  in  the  gaps.  The  confused  scattering  of  radiation  complicates 
feature  recognition.  One  persistent  trait  is  the  rounded  bulge  extending 
to  the  left  of  the  source  to  a  longitude  of  approximately  190°.  This 
changes  slightly  in  appearance  over  the  years,  but  not  enough  to  be 
considered  significant.  From  the  diagram  it  is  easy  to  see  that  secondary 
peaks  preceding  the  main  peak  in  the  histograms  of  the  previous  chapter 
are  caused  by  the  intermittent  storm-line  ridges.  Although  the  limited 
time  span  of  a  single  apparition  gives  incomplete  coverage  of  the  plane 
and,  therefore,  a  streaked  appearance,  some  of  the  streaks  may  be  enhanced 
by  passing,  through  a  naturally  higher-probability  area.  Thus,  the 
secondary  peaks  may  hint  of  an  underlying  feature  and  be  more  than  just 
a  random  superposition.  For  example,  the  two  histograms  combining  the 
negative  DE  years,  fairly  lengthy  time  spans,  show  similar  islands  of 
emission  preceding  the  main,  Io-related  part  of  the  source  in  Io  phase, 
although  their  Io-phase  boundaries  differ.  Some  other  years  have  storm 
bands  passing  through  the  same  area,  around  210°  of  Io  phase,  but 
1966-67  is  definitely  barren  at  that  location.  Separation  of  real  from 
artificial  enhancements  will  require  more  years  of  observation  and/or 
additional  data  from  the  past  filling  in  the  gaps  in  the  times  of 
monitoring. 

The  lower  level  of  non-Io-related  radiation  in  the  source-C  region 
permits  a  closer  examination  of  the  Io-related  qualities.  First 
impressions  of  Figure  IV-32  suggest  that  source  C  is  simply  a  fairly 
wide  swath  of  radiation  beginning  around  an  Io  phase  of  240°  and  sloping 
downward  in  the  normal  manner  for  storm  evolution.  The  end  of  this  swath, 
however,  usually  levels  out  and  sometimes  bends  back  upward,  as  in  1967-74. 
This  often  gives  the  bottom  edge  a  rounded  appearance.  The  overall  picture 


Figure  IV-32.  Comparison  of  histograms  of  relative  occurrence 

probability  for  2°  intervals  in  System  III  longitude  (65) 
and  the  phase  of  Io,  centered  on  the  Io-related  source-C 
region.  Combined  18,  20,  and  22-MHz  data  have  been  split 
into  individual  apparitions  for  years  of  positive  D£   and 
groups  of  apparitions  spanning  negative  D£  years.  The  number 
and  rank  of  the  characters  remain  the  same  as  in  Figure  IV-1 . 
All  probabilities  below  the  average  occurrence  probability 
for  each  diagram  have  been  replaced  by  blanks. 
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of  source  C  in  Figure  IV-27  suggests  another  feature  which  can  be  picked 
out  in  careful  observation.  There  is  a  bulge  at  the  upper  edge  of  the 
source  at  a  longitude  of  about  340°  and  an  Io  phase  near  230°.  This 
enhanced  area  is  especially  noticeable  in  the  active  years  from  1963-66. 
Part  of  the  responsibility  for  the  motion  of  the  source  in  Io  phase 
can  be  attributed  to  the  appearance  of  this  projection.  Also,  as  in 
source  B,  source  C's  lower  edge  extends  to  higher  values  of  Io  phase 
in  the  years  of  negative  D£.  The  secondary  peaks  in  the  one-dimensional 
histograms  around  200°  to  210°  of  Io  phase  are  seen  to  be  bands  of  storms 
differing  little  from  the  other  non-Io-related  radiation  in  C's  longitude 
range.  Again,  whether  the  secondary  peaks  suggest  significant  phenomena 
or  merely  result  from  the  limited  observational  cycle  is  hard  to  distinguish. 

If  the  radio  emission  from  Jupiter  is  dependent  upon  the  variables 
^III'  *Io'  and  DE  on1y'  then  apparitions  which  are  twelve  years  apart 
should  appear  similar,  since  they  have  nearly  equal  values  of  DV. 
Sufficient  coverage  of  the  ^tjj-4>Io  plane  during  each  apparition  is 
assumed.  In  Figures  IV-30  to  IV-32  this  corresponds  to  a  statement  that 
the  top  and  bottom  rows  of  histograms  should  match.  In  each  of  the  pairs 
some  vague  resemblances  can  be  picked  out  and  have  been  mentioned  as 
peculiar  to  the  low-activity  years  which  they  represent.  Overall  they  are 
far  from  exact  copies.  The  large  differences  in  the  amount  of  data  and  in 
the  observational  technique  before  and  after  a  12-year  interlude  probably 
explain  the  modification,  but  there  is  still  a  possibility  that  another 
variable  may  be  needed  to  account  for  all  discrepancies.  The  active  part 
of  the  cycle,  presently  taking  place,  should  make  it  possible  to  compare 
single  apparitions  having  relatively  high  storm  levels  in  the  near  future. 
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Many  first-  and  second-order  effects  have  been  revealed  in  this 
and  the  previous  chapter,  with  little  supposition  regarding  their 
relationship  to  the  decametric  emission  process.  This  speculation  has 
been  saved  until  now  so  that  the  basic  theories  and  the  implications  of 
the  effects  on  these  theories  might  be  presented  as  a  unit.  The  next 
chapter  is  dedicated  to  this  task. 


CHAPTER  V 
THEORETICAL  IMPLICATIONS 


5.1  Introduction 


The  toil  of  a  multitude  of  observers  scanning  Jupiter  over  a  score 
of  years  can  be  matched  only  by  the  man-hours  expended  in  explanation. 
A  pile  of  pertinent  theoretical  publications,  excluding  superficial 
opinions  proffered  in  the  midst  of  experimental  outcome,  is  daunting 
in  stature.  Comprehensive  review  of  said  articles  would  be  unduly  long 
and  unnecessary.  Fortunately,  the  field  has  been  blessed  by  periodic 
resumes  of  exceptional  value  (Warwick,  1964,  1967,  1970;  Ellis,  1965; 
Carr  and  Gulkis,  1969;  Smith,  1976).  On  page  1147  of  the  latest  of  these 
is  a  declaration  which  summarizes  the  inherent  worth  of  the  accumulation 
of  analyses:  ".  .  .we  believe  that  in  the  paradigmatic  nature  of  science 
a  successful  theory  will  emerge  from  the  gradual  modification  and  conso- 
lidation of  previous  work,  rather  than  from  radically  new  hypotheses." 
Assuming  that  a  completely  independent  solution  has  not  occurred  and  is 
nor  forthcoming,  it  would  seem  well  to  concentrate  on  comparing  results 
with  the  newer  theories,  which  have  profited  from  critiques  of  prior 
postulates  and  have  been  honed  by  the  expanding  assemblage  of  evidence. 
Indeed,  the  latest  offerings  are  more  given  to  numerical  statements  of 
expectations,  permitting  direct  testing  of  their  validity,  especially 
in  light  of  the  post-Pioneer  parameter  boundaries. 

In  the  course  of  Smith's  review,  the  theories  are  divided  according 
to  their  treatment  of  the  coupling  mechanism  which  must  provide  a 
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trigger  for  Io  to  stimulate  radiation;  the  instability  mechanism  and 
corresponding  amplification  process  required  to  explain  the  power  for 
both  Io-related  and  Io-independent  radiation;  and  the  subsequent 
propagation  of  the  radio  waves  through  the  Jovian  plasmasphere.  Emphasis 
in  this  dissertation  has  been  on  the  morphology  of  the  sources,  which 
is  most  directly  affected  by  the  last  category,  although  the  previous 
sections  certainly  have  some  influence.  Propagation  theory  is  the  least 
discussed  aspect  of  decametric  radiation,  making  direct  comparison  of 
prediction  with  observation  applicable  to  only  a  few  theoretical  models. 
Some  of  the  second-order  effects  uncovered  in  the  previous  chapters  relate 
in  an  indirect  fashion  to  other  theories,  and  this  too  will  be  considered 
wherever  possible.  Let  us  examine  the  results  in  the  context  of  several 
hypotheses. 

5.2  The  Source  B  Anomaly 

Evident  throughout  Chapter  IV  is  the  consistency  and  longevity  of 
the  Io-related  source-B  shape.  Its  uniformly  high  degree  of  probability 
and  control  by  Io  make  it  a  natural  starting  point  for  studying  the  nature 
of  Io's  influence.  Figure  V-l  summarizes  all  of  the  activity  in  the  Io-B 
region  for  the  frequencies  18,  20,  and  22  MHz.  In  the  histogram  the 
zones  are  2°  by  2°,  the  character  sequence  is  the  same  as  for  similar  plots 
in  Chapter  IV,  and  all  symbols  representing  probabilities  below  0.114 
have  been  replaced  by  blanks.  In  a  way,  the  shape  resembles  a  profile 
view  of  a  cellophane  tape  dispenser.  The  two  upper  projections  responsible 
for  this  morphology  are  the  interesting  part.  The  line  drawn  through  the 
diagram  illustrates  one  of  the  family  of  lines  of  equal  slope  which  trace 
the  change  of  longitude  and  Io  phase  with  time.  The  positioning  of 
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this  particular  line  shows  that  both  bulges  can  be  sampled  in  a  single 
night  of  monitoring.  Examination  of  the  histogram  occurrence  pro- 
babilities along  or  near  this  line  reveals  a  probability  of  0.13  at 
point  1  on  the  line,  a  zero  probability  at  point  2,  and  a  0.32  value 
at  the  third  position.  These  probabilities  seem  insignificant  when 
compared  with  the  maxima  nearby  (see  the  perspective  plots  of  Chapter  IV), 
but  there  is  still  a  definite  modulation  effect.  The  mechanism  producing 
the  anomaly  seems  to  allow  storms  to  occur,  then  shuts  them  off  completely, 
to  be  followed  by  an  even  more  probable  occurrence  of  storms— all  within 
the  space  of  an  hour  and  a  half. 

Does  the  radiation  actually  follow  this  scenario,  or  are  pre-  and 
post-null  probabilities  built  up  from  separate  storms  on  different  nights? 
Even  in  the  merged-frequency  data  presented  in  the  figure,  there  are  few 
nights  in  which  the  plotted  track  of  the  period  of  monitoring  would  lie 
close  to  the  depicted  line  and  span  both  enhancements  during  a  time  of 
reasonably  radio-quiet  conditions.  In  a  study  of  many  apparitions  of 
data,  25  nights  were  found  to  have  storms  contributing  to  zones  along  or 
near  the  line  in  question.  Of  these,  7  of  the  cases  exhibited  activity 
both  preceding  and  following  the  null  zone  at  point  2.  Thus,  in  some 
instances  the  radiation  does  turn  on  and  off,  and  then  on  again  from 
the  earth's  viewpoint.  On  none  of  the  nights  was  this  behavior  observed 
at  more  than  one  frequency.  This  could  reflect  the  difficulty  of 
receiving  the  low-probability  events,  or  it  could  mean  they  are  narrow- 
band by  nature.  It  would  be  interesting  to  examine  the  traces  of  these 
events  as  recorded  by  the  spectrograph  at  Boulder,  Colorado.  Unfortunately, 
according  to  the  catalog  of  Warwick  et  al.  (1975)  only  two  of  the 
incidents  were  received  in  Boulder.  Nevertheless,  the  study  of  the 
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Boulder  data  from  all  appropriate  nights  would  be  beneficial.  One  other 
comment,  which  will  be  discussed  in  more  detail  later,  is  that  all  but 
one  of  the  sets  of  pre-  and  post-null  storms  occurred  in  the  first  half 
of  an  apparition. 

The  two  bulges  of  source  B  are  mentioned  in  the  literature  indirectly. 
The  upper  left  protuberance  and  the  radiation  in  the  main  part  of  the 
source  directly  below  and  to  the  right  of  it  often  combine  to  form  a 
secondary  peak  at  approximately  125°  longitude  in  one-dimensional 
histograms  (cf.  Figures  III-l  and  1 1 1  -7 )  or  two-dimensional  contour 
plots  of  source  B.  This  has  prompted  some  researchers  (Smith,  et  al . , 
1962;  Goertz,  1971)  to  consider  the  two  peaks  as  separate  sources, 
labeled  B^  and  B2,  which  overlap  in  their  longitudinal  widths  making 
them  appear  as  one  source.  The  approximate  longitudinal  range  of  each 
source  is  indicated  in  Figure  V-l .  Goertz  reverses  the  order  of  labeling, 
but  in  this  work  B^  will  refer  to  the  dominant  peak  in  the  longitude 
histograms.  Several  arguments  can  be  made  for  and  against  the  division 
of  source  B. 

Close  scrutiny  of  Figure  V-l  reveals  two  nearly  equa.l  maxima  which 
have  been  circled.  The  first  is  symbolized  by  a  "0"  at  127°  in  longitude 
and  93°  in  Io  phase,  and  a  "1"  represents  the  other  at  151°  in  longitude 
and  89°  in  Io  phase.  The  dual  maxima  were  also  apparent  in  some  of  the 
contour  plots  mentioned  previously,  but  the  difference  in  Io  phase  between 
the  two  peaks  was  not  as  obvious.  In  general,  the  surrounding  zones 
emphasize  the  change  of  Io  phase  resulting  in  a  slight  upward  slant  of 
the  probability  contours  across  the  main  body  of  the  source.  This,  in 
turn,  causes  a  slight  indentation  along  the  bottom  edge  at  longitudes 
near  B-j .  If  B1  and  B2  are  separate  sources  they  seem  to  vary  slightly 
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in  Io-phase  position.  It  is  interesting  that  Bp,  the  secondary  peak  in 
the  one-dimensional  longitude  histograms,  has  an  Io-phase  position  equal 
to  the  one-dimensional  Io-phase  maximum. 

Division  of  source  B  has  also  been  indicated  by  the  spectral  changes 
in  the  B]  and  B2  bulges,  if  they  may  be  so  labeled,  brought  out  in  the 
previous  chapter.  The  B2  extension  is  prominent  at  30  and  27  MHz,  B, 
dominates  at  15  and  16  MHz,  and  both  are  conspicuous  for  frequencies  in 
between.  On  the  other  hand,  the  one-dimensional  longitude  histograms  of 
the  entire  source  remain  nearly  unchanged  by  frequency  differences. 

The  previous  chapter  has  also  discussed  the  occurrence  of  non-Io-B 
radiation,  especially  at  the  lower  frequencies,  and  its  usual  confinement 
to  longitudes  between  120°  and  180°.  This  means  the  non-Io-related  storms 
are  centered  on  the  B^  longitude  range  of  the  source.  Few  storms  are  found 
in,  or  extend  to,  the  lower  B2  longitudes.  One  exception  to  this  appears 
in  the  work  by  Desch  and  Carr  (1974)  at  26  MHz,  in  which  the  non-Io-B 
band  seems  to  be  centered  on  a  lower  longitude  between  the  B,  and  B„  peaks. 
It  is  possible  that  the  sensitivity  of  their  antenna  array  has  allowed 
detection  of  low-intensity  non-Io-B2  radiation.  In  combination  with 
Io-independent-B-j  storms  this  would  explain  the  decrease  in  the  central 
longitude. 

At  this  point  it  will  be  beneficial  to  consider  the  theory  of 
Cavaliere  and  Speranza  (1971),  which  assumes  the  existence  of  non-Io-B 
radiation.  For  explanation  of  Io-related  sources  the  theory  espouses  a 
now-common  interpretation  of  the  Io  influence  first  proposed  by  Dulk  (1965), 
in  which  the  Io-related  sources  originate  near  the  point  where  the 
magnetic  field  lines  passing  through  or  near  Io  intersect  the  Jovian  cloud 
tops  (the  "foot"  of  Io's  flux  tube),  and  the  radiation  is  emitted  nearly 
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perpendicular  to  these  lines.  The  longitude  beneath  Io  at  the  time  of 
storm  onset  and  the  direction  to  the  earth  determine  which  of  the  sources, 
if  any,  are   seen.  This  is  shown  in  Figure  V-2,  a  modification  of  a 
diagram  in  the  article  by  Cavaliere  and  Speranza.  The  figure  looks  down 
on  the  north  pole  of  the  axis  of  rotation  and  Io  is  placed  at  two  positions 
of  interest.  At  the  foot  of  the  Io  flux  tube,  the  curved  line  connecting 
Io  to  Jupiter,  the  direction  of  emission  of  Io-related  radiation  is 
represented  by  arrows  nearly  orthogonal  to  the  flux  lines.  If  the  emissions 
reach  the  earth,  their  identification  as  a  particular  Io-related  source 
depends  on  the  position  of  Io  and  the  longitude  facing  the  earth  at  the 
time  of  storm  generation.  The  two  examples  of  Io  position  were  chosen  for 
their  angular  distances  from  the  north  magnetic  pole  (NMP),  assumed  to 
point  toward  a  longitude  of  200°  as  represented  by  the  dashed  line.  They 
yield  the  Io-related  sources  indicated  at  the  tips  of  the  arrows  at  the 
flux  tube  feet.  Thus,  Io-A  and  Io-Bp  storms  may  have  the  same  point  of 
origin,  but  they  are  distinguished  by  their  direction  of  propagation.  For 
small  changes  in  the  given  positions  of  Io,  the  identification  of  the  sources 
will  remain  the  same.  Larger  changes  will  bring  a  cutoff,  of  the  source  by 
some  type  of  modulation  mechanism. 

According  to  Cavaliere  and  Speranza,  the  non-Io  sources  are  turbulent 
plasma  regions  beaming  radiation  radially  outward  from  the  planet.  The 
theory's  modulation  mechanism  permits  occurrence  of  Io-related  storms  only 
when  Io  is  over  a  longitude  corresponding  to  the  non-Io-A  and  non-Io-B 
sources.  The  non-Io-A  and  non-Io-B  sources  surround  the  magnetic  pole  as 
depicted  by  the  shaded  regions  and  labeled  arrows  in  Figure  V-2.  The 
Io-related  storms  are,  therefore,  stronger  emissions  generated  in  the 
same  turbulent  plasma  regions  which  produce  non-Io  radiation.  The  major 
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Figure  V-2.  View  of  Jupiter  from  directly  over  the  north  pole  showing 

sample  beaming  directions  of  the  decametric  sources  according 
to  an  extension  of  the  theory  of  Cavaliere  and  Speranza  (1971) 
The  longitude  of  the  north  magnetic  pole  (NMP)  is  assumed  to 
be  ^200°.  The  geometry  of  source  emission  for  the  Io-related 
sources  applies  to  many  other  theories  as  well. 
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drawback  to  this  hypothesis  is  the  failure  to  account  for  non-Io-C 
radiation.  Extension  of  the  theory  may  explain  this  source  and  imply 
existence  of  non-Io-B2.  Cavaliere  and  Speranza  state  that  the  Io-dependent 
A  and  B  sources  are  positioned  symmetrically  about  the  longitude  of  the 
north  magnetic  pole  ( % 200°±50° ) .  An  extension  of  the  theory  to  include 
a  separate  non-Io-C  source  at  a  longitude  centered  around  310°  would 
suggest  by  symmetry  a  counterpart  near  90°.  This  latter  source  could  be 
interpreted  as  non-Io-B2  or  perhaps  non-Io-D,  if  they  exist.  Some  evidence 
for  a  band  of  Io-independent  storms  centered  on  the  90°  longitude  can  be 
found  in  Figure  IV-24,  the  two-dimensional  histogram  for  10-MHz  data.  The 
problem  with  the  extension  of  the  theory  is  that  one  would  expect  Io-related 
radiation  to  occur  in  the  non-Io-C  and  non-Io-B?  regions.  The  angular 
distance  of  these  sources  from  the  north  magnetic  pole  may  weaken  the 
plasma  turbulence  and  prevent  this  emission.  In  general,  the  modulation 
mechanism  proposed  by  Cavaliere  and  Speranza  works  well  for  explanation  of 
approximate  source  positions,  but  the  added  assumptions  necessary  to  explain 
non-Io-C  radiation  and  the  exact  source  boundaries  would  stretch  its 
credibility. 

In  Figure  V-2  the  perpendicular-beaming  model  suggests  that  Io-related 
sources  B2  and  A  originate  on  the  opposite  side  of  the  north  magnetic  pole 
from  sources  B]  and  C.  The  reason  for  this  is  clearer  in  the  next  figure. 
Figure  V-3  replots  the  22-MHz  data  of  Figure  IV-9  in  terms  of  System  III 
longitude  versus  the  longitude  beneath  Io  at  the  time  of  emission.  The 
sources  are  altered  in  shape,  but  still  recognizable,  and  they  are  more 
closely  grouped.  A  line  drawn  along  a  sub-Io  longitude  near  225°  separates 
not  only  sources  A  and  C,  but  also  B]  and  B2.  For  purposes  of  radiation- 
beaming  symmetry  most  theories  incorporate  some  mechanism  which  places  the 


Figure  V-3.  Histogram  of  relative  occurrence  probability  for  2° 

intervals  in  System  III  longitude  (65)  and  the  sub-Io 
System  III  longitude  (65)  for  all  22-MHz  data.  The 
characters  represent  percentiles  of  the  maximum  occurrence 
probability.  Their  number  and  rank  are  listed  at  the 
bottom. 
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origin  of  Io-related  storms  at  a  longitude  approximately  15°  lower  than 
Io's  position.  This  is  equivalent  to  moving  the  sources  upward  by  15° 
in  Figure  V-3,  and  it  changes  the  separation  line  to  210°,  a  longitude 
which  is  close  to  that  of  the  north  magnetic  pole. 

The  symmetry  about  the  dividing  line  can  be  taken  one  step  further 
by  a  comparison  of  Figure  V-4,  an  enlarged  view  of  the  Io-related 
source-C  region  with  the  B,  side  of  source  B  in  Figure  V-l.'  An  overlay 
of  the  two,  with  a  longitude  of  185°  and  an  Io  phase  of  90°  in  the 
source-B  picture  aligned  with  360°  and  240°,  respectively,  in  the 
source-C  diagram,  reveals  many  similarities.  The  B,  bulge  and  the  pro- 
tuberance above  source  C  coincide,  as  does  the  downward  slant  of  the  edges 
just  to  the  right  of  the  bulges.  Ignoring  the  scattered  non-Io-C 
data,  the  extent  of  each  source  in  Io  phase  is  comparable.  The  sole 
difference  of  major  import  is  the  "tail"  on  the  right  side  of  C.  Source 
B  sometimes  has  tails,  extending  into  the  non-Io-A  longitudes,  but  they 
never  have  an  upward  turn  like  that  of  source  C. 

In  the  comparison  of  source  B,  with  source  C,  one  is  reminded  of 
the  Io-phase  sinusoidal  drifts  of  sources  B  and  C  discussed  in  Chapter  III. 
Source  B  varies  its  Io-phase  peak  in  the  same  direction,  but  less  than 
half  as  much  as  source  C.  Perhaps  this  is  caused  by  a  large  motion  of 
B1  like  that  of  source  C,  coupled  with  a  lack  of  motion  of  B?  like  that 
of  A.  Alignment  of  the  changing  source  positions  by  correcting  the  data 
for  the  sinusoidal  drifts  should  reveal  such  a  situation  by  a  concentration 
of  the  maximum  probability  areas  at  a  single  point.  In  accordance  with 
the  results  presented  in  Chapter  III,  all  18,  20,  and  22-MHz  storms 
were  modified  in  observed  Io-phase  by  subtracting  2.9xDp  degrees.  Dp 
was  determined  from  the  date  of  the  storm,  and  the  factor  2.9  is 
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the  proportionality  constant  determined  in  Chapter  III  relating  DF  to 
Io  phase  for  source  C.  At  the  same  time,  to  study  the  longitudinal 
motion  exhibited  by  source  A,  the  longitudes  of  all  storms  were 
corrected  by  -2.7xDE  degrees,  a  factor  ascertained  in  Chapter  III  from 
the  non-Io-A  drift.  The  results  are  shown  in  Figure  V-5  and  can  be 
compared  with  Figure  IV-27.  Changes  between  the  two  figures  are  nearly 
insignificant.  The  non-Io-related  radiation  seems  to  be  more  concentrated 
within  each  source's  longitude  range,  implying  that  each  source  does  have 
its  own  Io-independent  segment  and  that  all  non-Io  sources  drift  in 
longitude  in  a  similar  manner.  If  the  drifts  are  all  similar,  they  can 
be  explained  by  a  variable  period,  but  this  is  difficult  to  reconcile  with 
DE  dependence.  Explanation  in  terms  of  beaming  geometry  seems  more 
likely,  but  similar  changes  in  the  geometry  for  all  non-Io  sources  are  hard 
to  justify.  Although  the  Io-related  source-C  radiation  is  more  concentrated, 
the  B1  counterpart  is  not.  This  suggests  source  B  moves  as  a  unit  in  its 
Io-phase  drift,  whether  or  not  it  is  composed  of  separate  sources. 

One  theory  proposes  the  existence  of  separate  B,  and  B?  sources. 
Goertz  (1973)  hypothesized  the  triggering  of  Io-dependent. storms  by 
magnetohydrodynamic  (MHD)  waves  traveling  down  magnetic  flux  lines  which 
have  been  perturbed  by  Io's  passage.  The  waves  encounter  dense  plasma 
in  the  lower  parts  of  the  Jovian  ionosphere  and  stimulate  coherent 
cyclotron  radiation.  The  radiation  is  beamed  nearly  perpendicular  to  the 
flux  tube  so  that  the  propagation  geometry  for  Io-related  sources  is 
equivalent  to  the  configuration  shown  in  Figure  V-2.  In  this  case,  the 
modulation  mechanism  limiting  the  extent  of  the  sources  involves  the 
amount  of  amplification  of  the  radiation.  The  amplification  factor  first 
increases  and  then  decreases  rapidly  as  Io  moves  away  from  the  longitude 


Figure  V-5.  Histogram  of  relative  occurrence  probability  for  2° 

intervals  in  System  III  longitude  (65)  and  Io  phase  for 
all  18,  20,  and  22-MHz  data.  The  storms  have  been 
corrected  for  longitudinal  and  Io-phase  drifts  according 
to  the  method  described  in  the  text.  The  characters 
represent  percentiles  of  the  maximum  occurrence  probability. 
Their  number  and  rank  are  listed  at  the  bottom. 
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of  the  north  magnetic  pole,  since  plasma  characteristics  change  with 
angular  distance  from  the  pole.  The  radiation  is  lessened  when  Io  is 
over  the  pole  because  amplification  depends  inversely  upon  the  area  of 
spread  of  the  MHD  wave  and  this  is  greatest  at  the  pole.  The  result  of 
these  conditions  is  a  maximum  of  amplification  on  each  side  of  the  pole, 
and  the  Io-related  sources  are  again  beamed  as  in  Figure  V-2.  The  theory 
may  also  explain  the  longitudinal  merging  of  sources  A  and  B  at  10  MHz 
observed  in  earlier  chapters,  since  this  radiation  should  be  emitted  farther 
from  the  planet  where  amplification  may  be  limited  to  points  closer  to 
the  pole..  Goertz  predicts  source  positions  in  a  contour  plot  of  reception 
probability  for  Ajn  versus  <J>Iq.  Peaks  of  probability  occur  near  the 
normal  positions  for  sources  A,  C,  B]  and  B2  and  Goertz  labels  them  as 
such. 

The  same  pair  of  amplification  maxima  produce  sources  B,  and  B?  when 
the  phase  of  Io  is  90°,  and  sources  A  and  C  for  a  240°  Io  phase,  and  yet 
sources  A  and  C  have  greater  separation  in  longitude  than  B,  and  B?.  The 
theory  accounts  for  this  by  stating  that  sources  A  and  C  are  produced 
when  Io  is  on  the  hotter,  afternoon  side  of  Jupiter.  Heating  of  the 
ionosphere  makes  it  thicker,  and  this  in  turn  causes  the  amplification 
factor  to  peak  farther  from  the  pole.  A  simple  test  can  be  used  to  study 
this  hypothesis. 

Assume  that  the  separation  of  the  amplification  maxima  increases 
linearly  as  the  Jovian  day  proceeds.  The  approximate  angular  distance 
in  longitude  between  sources  A  and  C  has  been  measured  from  the  different 
apparitions  of  the  18,  20,  and  22-MHz  data,  and  it  averages  around  70°. 
Measurement  of  the  average  B^Bg  separation  is  more  difficult,  but  it  is 
roughly  30°.  This  means  there  is  a  40°  increase  in  separation  over  the 
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150°  change  in  Io's  position  from  the  dawn  side  to  near  dusk.  Although 
Io-related  sources  A  and  C  are  received  on  earth  when  Io  is  near  240° 
from  superior  geocentric  conjunction,  Io's  position  from  superior  helio- 
centric conjunction  at  those  times  changes  with  the  sun-Jupiter-earth 
angle.  Let  the  change  in  this  angle  from  one  point  on  the  earth's  orbit 
to  another  be  represented  by  the  angle  B.  For  example,  the  maximum  change 
occurs  between  the  two  points  of  quadrature  on  the  earth's  orbit,  giving  a 
23°  value  of  3,  as  illustrated  by  Figure  V-6. 

If  the  earth's  orbit  is  divided  into  the  first  and  second  halves 
of  a  Jovian  apparition,  then  the  average  change  in  the  sun-Jupiter-earth 
angle  between  the  positions  of  earth  during  the  first  half  and  the  locations 
in  the  second  half  is  15°.  This  means  that  for  a  given  geocentric  position 
Io  is  an  average  of  15°  farther  from  superior  heliocentric  conjunction  in 
the  first  half  of  an  apparition  than  in  the  second  half.  If  the  linear- 
increase-in-separation  theory  were  strictly  true,  sources  A  and  C  should  be 
approximately  4°  farther  apart  in  the  histogram  from  the  first  half  of 
an  apparition  than  in  that  from  the  second  half.  ,In  studying  the  18, 
20,  and  22-MHz  source  centroids  computed  before  and  after. opposition,  the 
average  difference  was  1.8°.  This  result  is  inconclusive  since  the  large 
number  of  subtractions  involved  in  the  computation  causes  an  error  of  the 
same  magnitude  as  the  result.  The  small er-than-expected  difference  is 
more  realistic  in  the  sense  that  the  temperature  of  the  ionosphere  will 
not  continually  increase,  but  will  start  to  decline  at  some  time  in  the 
Jovian  afternoon.  Also,  conditions  are  more  favorable  for  radio  reception 
near  opposition  and  this  will  effectively  decrease  the  average  value  of  6. 

As  a  side  note,  in  the  course  of  the  study  outlined  above  it  was  found 
that  the  average  occurrence  probability  decreased  from  a  value  of  0.055 
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Figure  V-6.  Illustration  of  the  maximum  value  of  the  angular  change 
in  the  sun-Jupiter-earth  angle  from  the  first  half  of  an 
apparition  to  the  second  half.  When  the  earth  is  in  the 
position  shown  during  the  first  half  of  an  apparition,  Io 
must  be  in  position  1  for  Io-related  source-B  radiation  to 
be  received.  The  second  position  of  Io  is  necessary  for 
reception  of  Io-B  radiation  at  the  earth's  indicated  second- 
half  position. 
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in  the  first  half  of  all  apparitions  to  0.050  for  the  second  half.  This 
is  not  surprising  since  Jupiter  is  positioned  in  the  early  morning  sky 
during  the  first  half  of  an  apparition,  and  the  conditions  for  radio 
reception  are  generally  better  at  that  time.  The  change  in  the  probability 
does,  however,  provide  initial  confirmation  of  an  effect  observed  by 
Gruber  (1965).  He  divided  the  apparitions  into  15-day  periods  and 
discovered  a  large  decrease  in  the  reception  probability  just  after 
opposition.  Analysis  of  the  Florida-Texas  data  base  in  a  similar  vein 
might  be  a  worthwhile  future  project. 

One  major  difficulty  with  the  Goertz  model  seems  to  be  the  polarization 
He  predicts  a  larger  ratio  of  left-hand  polarization  as  the  sub-Io 
longitude  nears  the  north  magnetic  pole.  This  is  based  on  the  assumption 
that  the  radiation  is  emitted  in  the  extraordinary  mode  with  right-hand 
polarization  emanating  from  the  northern  hemisphere  and  left-hand 
polarization  from  the  southern  hemisphere.  The  reverse  seems  to  be  true, 
because  source  C,  the  only  major  source  of  left-hand  polarization  that 
Goertz  has  included,  is  farther  from  the  pole  in  sub-Io  longitude  than 
the  dominant  right-handed  B]  source.  If  the  left-hand  polarization  of 
source  D  were  added  to  the  picture  the  symmetry  of  polarization  about  the 
pole  longitude  would  be  more  complete  since  Io-related  D  lies  at  sub-Io 
longitudes  farther  from  the  pole  than  sources  A  and  B„.  In  general,  the 
source  Goertz  has  considered  to  be  B2  might  be  better  labeled  as  source 
D,  or  a  combination  of  the  two.  Indeed,  the  majority  of  theoreticians 
have  called  the  radiation  which  has  the  beaming  characteristics  of  B?, 
as  illustrated  in  Figure  V-2,  source  D. 

In  summary,  the  separation  of  B]  and  B2,  if  such  a  division  exists, 
is  probably  a  phenomenon  associated  with  the  location  of  the  north 
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magnetic  pole  and  its  effect  on  the  surrounding  plasma.  Further  arguments 
concerning  the  link  between  sources  Bj,  By,   and  D  will  be  presented 
shortly,  but  first  let  us  consider  some  of  the  other  theories  in  regard 
to  wave  propagation. 

5.3  Beaming  of  the  Sources 

Hypotheses  concerning  the  exact  mechanism  responsible  for  production 
of  Io-related  storms  span  a  wide  range  of  possibilities.  We  are  more 
concerned  with  the  subsequent  propagation  of  the  radiation  away  from  the 
point  of  origin.  In  discussion  of  emission  geometry  theories  fall  into 
two  categories— those  in  which  the  beaming  is  perpendicular,  or  nearly 
perpendicular,  to  the  magnetic  flux  lines  passing  through  Io,  and  those 
involving  beaming  in  or  close  to  the  direction  of  the  magnetic  flux  tube. 
The  former  category  is  by  far  the  more  popular.  One  notable  exception 
to  these  two  classes  is  the  model  of  Gledhill  (1967).  In  this  case  the 
point  of  origin  was  assumed  to  be  near  Io,  as  the  satellite  perturbed  a 
disk  of  plasma  while  passing  through  it.  This  theory  is  generally 
disregarded  today  since  it  requires  an  inordinately  high  plasma  density 
in  Io's  vicinity,  which  is  inconsistent  with  Pioneer  spacecraft  findings. 
Through  the  years,  Ellis  (1965)  and  McCulloch  (1971)  have  developed 
a  theory  which  represents  the  case  of  emission  aligned  with  the  magnetic 
field.  In  this  model  Io  is  assumed  to  generate  whistler  radiation  and, 
as  described  by  McCulloch,  the  waves  propagate  forward  and  backward  to 
higher  and  lower  longitudes  than  the  sub-Io  meridian.  The  finite  time  of 
travel  in  combination  with  the  rotation  of  Jupiter  under  Io  causes  the  rays 
to  intersect  the  planet  at  asymmetrical  angular  distances  from  the  sub-Io 
longitude.  It  is  also  assumed  that  the  rays  must  intersect  the  planet  at 
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magnetic  latitudes  between  75°  and  80°  so  that  the  radiation  generated 
will  be  beamed  toward  earth.  The  most  efficient  time  for  this  to  happen 
is  when  Io  passes  the  longitude  of  the  north  magnetic  pole.  McCulloch 
finds  that  the  sum  of  these  conditions  yields  two  peaks  of  radiation 
intensity  corresponding  to  the  positions  of  Io-related  sources  A  and  B. 
A  Ain  versus  $Jo   contour  plot  of  the  predicted  relative  intensity  near 
30  MHz  closely  resembles  the  observed  shape  of  source  B.  There  is  a 
projection  similar  to  the  B2  bulge  found  at  30  and  27  MHz.  McCulloch 
also  computes  a  one-dimensional  Io-phase  histogram  of  relative  intensity. 
The  source-B  peak  follows  observed  asymmetry  since  it  rises  slowly  and  falls 
off  more  rapidly  with  increasing  Io  phase.  Unfortunately,  the  predictions 
for  source  A  are  not  as  successful.  According  to  the  theory,  A  behaves 
as  a  mirror  image  of  B  about  the  sub-Io  longitude  of  the  pole.  The  two- 
dimensional  contour  plot  shows  a  projection  to  the  high  Io-phase,  high 
longitude  side  of  source  A  that  does  not  match  observations.  The  one- 
dimensional  Io-phase  peak  rises  rapidly  and  has  a  slower  descent,  which 
also  opposes  reality.  It  can  be  argued  that  the  plots  are  of  relative 
intensity  rather  than  occurrence  probability,  but  the  probability  should 
roughly  follow  the  intensity  profiles.  McCulloch  does  not  discuss 
sources  C  and  D,  but  addition  of  southern  hemisphere  effects  would  probably 
be  necessary  to  explain  them.  The  model  has  its  shortcomings,  but  its 
strengths  should  be  kept  in  mind  as  an  alternate  method  of  beaming. 

The  great  majority  of  thenewer  theories  subscribe  to  the  perpendicular 
beaming  patterns  illustrated  in  Figure  V-2.  This  requires  some  explanation 
of  the  asymmetry  of  the  position  of  Io  about  the  earth-Jupiter  line  at  the 
time  of  Io-controlled  storms.  Most  theories  invoke  some  system  which  allows 
the  point  of  origin  to  occur   at  a  longitude  15°  lower  than  that  under  Io, 
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since  this  would  symmetrize  the  starting  positions.  As  an  example, 
we  have  already  discussed  the  Goertz  model  in  which  Io  perturbs  the 
magnetic  field  lines  by  its  passage.  MHD  waves  change  the  orientation 
of  the  flux  lines  and  cause  the  disturbance  created  by  Io  to  penetrate 
the  ionosphere  at  a  longitude  15°  lower  than  the  sub-Io  longitude.  Another 
explanation  involves  the  assumption  of  high  conductivity  for  Io  as 
suggested  by  Goldreich  and  Lynden-Ball  (1969)  or  Piddington  and  Drake  (1968) 
In  this  case,  charged  particles  are  accelerated  along  the  magnetic  flux 
tube  moving  with  Io.  The  motion  of  the  flux  tube  gives  rise  to  a  magnetic 
torque  which  bends  the  field  lines  forward  by  15°.  Restoration  of 
symmetry  by  bending  of  the  field  line  is  illustrated  in  Figure  V-2  by  the 
curves  connecting  Io  to  Jupiter.  Finally,  for  one  other  example,  Wu  et  a! . 
(1973)  assumed  the  disturbance  to  occur  at  the  sub-Io  longitude,  but  slow 
propagation  of  the  radiation  through  the  corotating  magnetosphere  prevents 
escape  until  the  sub-Io  longitude  has  increased  by  15°. 

If  the  Io-related  radiation  effectively  originates  from  longitudes  75° 
on  either  side  of  the  earth-Jupiter  line,  one  must  still  explain  why  the 
storms  are  maximized  for  the  75°  angle.  The  dominant  theory  for  the 
emission  mechanism  is  coherent  cyclotron  radiation  by  electrons  moving  in 
bunches  up  or  down  the  flux  lines.  In  the  frame  of  reference  moving  with 
the  electrons  the  radiation  is  emitted  perpendicular  to  the  field  lines. 
To  the  outside  observer  radiation  emerges  in  a  conical  sheet  about  the 
direction  of  motion  at  an  angle  near  the  required  75°.   If  the  cone  is 
beamed  toward  the  planet  it  will  be  reflected  without  loss  of  cone  angle 
by  increased  plasma  densities  at  lower  levels. 

As  a  final  step  in  using  perpendicular  beaming,  some  method  of 
modulation  must  be  invoked  to  explain  why  the  cone  of  radiation  can   be 
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received  on  earth  only  for  particular  ranges  of  sub-Io  longitude.  We 
have  seen  that  Cavaliere  and  Speranza  limit  the  emission  regions  by  the 
boundaries  of  the  non-Io-related  A  and  B  sources,  and  Goertz  believes 
sufficient  amplification  of  the  radiation  can  occur  only  when  Io  is  near 
the  north  magnetic  pole.  Goldreich  and  Lynden-Bell  theorized  that 
radiation  is  strong  only  in  parts  of  the  emission  cone.  Only  at  times 
when  the  active  portion  of  the  cone  points  toward  the  observer  will 
storms  be  received.  In  the  original  work  by  Dulk  (1965)  it  was  suggested 
that  the  orientation  of  the  magnetic  field  may  allow  intersection  of  the 
radiation  cone  and  the  ecliptic  plane  containing  the  earth  only  when  Io 
is  in  northern  magnetic  latitudes  and,  therefore,  in  the  longitudes 
surrounding  the  north  magnetic  pole. 

The  diversity  of  theories  which  share  similar  explanations  of 
beaming  geometry  stresses  the  importance  of  this  aspect  of  emission. 
It  is  possible  to  generalize  these  similarities  and  use  them  in  conjunction 
with  models  of  the  magnetic  field  in  order  to  predict  the  source  structure. 
Schatten  and  Ness  (1971)  calculated  emission  directions  for  both  the 
dipolar  magnetic  offset  from  the  center  of  Jupiter  as  suggested  by 
Warwick  (1963)  and  a  centered,  tilted  magnetic  dipole.  They  assumed  the 
radiation  emanates  in  a  cone  at  angles  between  80°  and  100°  to  the  line 
of  magnetic  flux.  Predictions  of  probability  for  30°  by  30°  intervals  in 
Io  phase  and  System  III  longitude  were  computed.  General  agreement  with 
the  positions  of  the  four  Io-related  sources  was  found,  but  the  constraints 
were  insufficient  for  determination  of  structure. 

Recent  measurements  of  the  Jovian  magnetic  field  geometry  by  two 
magnetometers  aboard  the  Pioneer  spacecraft  (Smith  et  al.,  1975;  Acuna 
and  Ness,  1976)  have  provided  a  more  definitive  basis  for  calculations. 
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Both  sets  of  investigators  described  the  magnetic  field  in  terms  of  a 
spherical  harmonic  analysis  of  the  magnetic  potential  as  represented 
below. 
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The  superscript  e  refers  to  external  contributions  to  the  field,  i  stands 

for  internal  field  sources,  a  is  the  radius  of  Jupiter,  and  k  is  the 

maximum  order  to  which  the  analysis  is  applied.  The  Pm(cose)  are  the 

Schmitt-normalized  associated  Legendre  functions,  and  Gm,Hm,gm,  and  hm 

n  n  Jn      n 

are  the  Schmitt  coefficients.  Finally,  e  and  4,  are  the  co-latitude  and 
Jovian  east  longitude,  respectively.  The  magnetic  field  can  be  determined 
by  computing 
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Usually  the  external  source  coefficients  contribute  negligibly  (see 
Smith  et  al. ,  1975)  for  points  close  to  the  planet  and  they  are 
neglected  in  the  present  analysis.  Each  group  of  investigators  carried 
out  their  computation  to  include  octopolar  moments  of  the  field  (k=3). 
The  calculated  coefficients  are  summarized  in  Table  V-l. 

Once  the  coefficients  are  known,  it  is  a  simple  matter  to  compute 
the  location  of  the  intersection  with  the  Jovian  surface  of  the  magnetic 
field  lines  passing  through  Io.  Beginning  at  Io,  the  flux  lines  are 


Table  V-l 
Octopolar  Field  Coefficients 

Smith  et  al .  (1975)  Acuna  and  Ness  (1976) 

9l                4.129  4.218 

91  -0.492  -0.440 

hl                0.531  0.563 

g2                0.042  _0.203 

g2               -0.738  -0.871 

92  °-324  0.331 

h2               -0.050  -0.037 

2 

n2               -0.381  -0.402 

9§                0.092  -0.233 

93  "0.413  -0.357 
93  0.335  0.506 
93               -0.239  -0.292 


h3  -0.084  -0.463 

h3  0.002  0.096 

h3  0.118  0.233 
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followed  in  iterative  steps  to  the  points  of  interest  at  or  near  the 
planet's  cloudtop  level.  Comparison  of  the  source  morphology  with 
sample  beaming  patterns  computed  at  selected  points  is  now  possible. 
The  radiation  will  be  assumed  to  emerge  from  the  selected  location  in 
a  conical  sheet  opening  at  an  angle  y   about  the  magnetic  flux  line.  In 
coherent  cyclotron  radiation  the  angle  y   is  determined  by  the  speed  v  of 
the  electron  along  the  magnetic  field. 

Y  =  2  "  c  (V-5) 

Inward  or  .outward  pointing  of  the  cone  axis  does  not  matter  since 
inward-directed  waves  will  be  reflected  at  lower  levels,  leaving  the 
cone  angle  unchanged.  The  effects  of  stop  zones  and  of  refraction  in 
the  plasmasphere  will  not  be  included.  They  may  be  of  importance  to  the 
beam  propagation,  but  computation  of  their  influence  is  a  complicated 
process  in  which  accuracy  depends  upon  the  validity  of  the  assumed  model 
of  the  plasma  distribution.  If  the  plasma  were  a  major  factor  in  the 
radiation  propagation,  any  instabilities  in  the  plasma  should  be 
noticeable  as  large  changes  of  source  structure.  Consistency  of  the 
source  morphology  over  many  years  suggests  that  modulation  of  emission 
direction  by  the  plasma  is  small.  Consequently,  it  will  be  neglected  in 
this  simple  analysis  of  beaming  geometry. 

Figures  V-7a  and  V-7b,  adapted  from  Goertz  (1973),  illustrate  the 
relevant  variables  for  determination  of  the  intersection  of  the  radiation 
cone  with  the  plane  of  the  ecliptic,  which  is  the  earth's  viewpoint.  The 
radiation  is  beamed  at  an  angle  u  on  either  side  of  the  projection  of 
the  magnetic  field  line  on  the  ecliptic  plane.  If  t  is  the  angle  between 
the  field  line  and  the  plane  of  the  ecliptic,  then  u  can  be  computed 
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from  the  equation 

2      2    2      2 
tan  u  =  tan  y   cos  t  -  sin  t  (V-6) 

Unless  the  emission  cone   does  not  intersect  the  ecliptic  plane,  storms 
will  be  directed  toward  the  earth  for  all  sub-Io  longitudes.  As  discussed 
previously,  some  form  of  modulation  must  be  introduced  to  limit  the 
emission  range. 

We  will  adopt  the  limit  proposed  by  Goldreich  and  Lynden-Bell , 
that  the  radiation  is  emitted  only  within  the  portion  of  the  cone  defined 
by  the  angle  6  in  Figure  V-7c.  In  the  figure,  the  cone  is  projected  onto 
a  plane  perpendicular  to  the  cone  axis.  Maximum  emission  occurs  in  the 
portion  of  the  cone  parallel  to  the  lines  of  magnetic  latitude,  repre- 
sented by  the  dashed  line.  Beyond  an  angle  6  on  either  side  of  the 
magnetic  latitude  line  radiation  is  negligible.  Reception  of  storms  will 
be  possible  only  when 

tan  t  -  tan  y   sin  <5  (V-7) 

We  will  also  assume  that  the  flux  lines  will  lead  Io  by  15°  in  longitude. 

This  is  effected  by  increasing  the  sub-Io  longitude  by  15°  in  the 

calculations.  This  is  not  the  same  as  a  warping  of  the  magnetic  field 

line,  but  the  results  should  not  be  too  different.  The  final  assumption 

requires  emission  at  the  point  where  the  electron  gyrofrequency  determined 

by  the  magnetic  field  strength  equals  the  frequency  of  the  emitted  wave. 

Since  comparison  will  be  made  with  18,  20,  and  22-MHz  source  structure, 

the  gyrofrequency  of  20  MHz  has  been  chosen. 

Figure  V-8  exhibits  the  calculated  values  of  ATTT  and  d>.  for 

111     Io 

which  radiation  can  be  received  from  the  northern  hemisphere  foot  of  Io's 
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flux  tube  with  the  given  assumptions.  The  left  diagram  corresponds  to 
use  of  the  Acuna  and  Ness  magnetic  field  model,  and  the  diagram  on  the 
right  utilizes  the  model  of  Smith  et  al .  The  points  have  been  computed 
for  20°  intervals  in  sub-Io  longitude  and  the  lines  interpolate  these 
values.  Rough  outlines  of  the  mid-decametric  sources,  based  on  diagrams 
like  Figures  V-l  and  V-4,  have  been  provided  for  reference.  To  help 
visualize  the  effects  of  different  values  for  the  variables,  three  sets 
of  parameters  have  been  used. 

The  sources  are  well -covered  by  variations  in  cone  angle  between  70° 
and  90°.  In  source  B,  especially,  the  sharp  cutoff  of  the  lower  edge  is 
followed  by  the  70°  curve.  From  equation  V-5  this  can  be  interpreted  as 
a  cutoff  of  about  30  keV  for  electron  motion  along  the  field  line.  As 
originally  suggested  by  Dulk  (1965),  a  cone  angle  of  79°  passes  through  or 
y/ery   near  the  areas  of  maximum  probability.  The  angle  of  90°,  corresponding 
to  the  least  motion  parallel  to  the  field,  covers  the  lower-probability 
upper  part  of  source  B.  As  a  consequence,  the  asymmetry  of  the  source-B 
peak  in  Io  phase  implies  a  parallel  velocity  distribution  for  the 
radiating  electrons.  The  number  of  electrons  builds  up  slowly  as  parallel 
speed  increases  and  reaches  a  maximum  for  a  cone  angle  of  79°,  followed 
by  a  rapid  decline  in  the  number  at  higher  energies.  The  reverse  would 
be  the  case  for  the  opposite  lower-frequency  asymmetries. 

From  Figure  V-8  we  see  that  for  any  given  asymmetry  in  Io  phase  for 
B,  A  and  C  should  have  opposite  skews  since  the  cone  angle  increases  with 
decreasing  Io  phase  in  those  sources.  This  opposing  symmetry  is  not  found. 
Io-A  and  Io-C  have  a  higher  mixture  of  non-Io  radiation  to  confuse  the 
identification,  but  their  direction  of  skew  is  still  unmistakeable.  One 
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possible  differentiating  influence  is  the  aforementioned  occurrence  of  A 
and  C  on  the  hotter  afternoon  side  of  the  planet  which  may  change  the 
parallel  velocity  distributions. 

For  both  magnetic  field  models  in  Figure  V-8  part  of  source  A  extends 
beyond  the  70°  cutoff,  implying  higher  speeds  for  the  electrons  along 
the  field  line,  but  this  may  be  the  confusion  of  the  non-Io  radiation 
blending  into  the  source  at  those  longitudes.  The  locations  of  the  B, 
and  B2  bulges  intimate  cone  angles  greater  than  90°  which  may  in  a  sense 
be  true.  This  discussion  will  be  more  apropos  in  the  next  section. 

The  longitudinal  width  predicted  for  the  sources  depends  upon  6. 
Inspection  of  equation  V-7  shows  confinement  of  radiation  to  the  part  of 
the  cone  determined  by  <S  to  be  the  same  as  limiting  the  activity  to  a 
symmetrical  longitude  range  about  the  north  magnetic  pole,  since  the  angle 
t  will  be  minimized  at  the  pole  longitude.  Thus,  the  theory  of  Goldreich 
and  Lynden-Bell  is  similar  in  effect  to  the  models  of  Goertz,  Cavaliere 
and  Speranza,  and  many  others,  except  in  the  case  where  y  is  close  to  90°. 
In  that  situation  6  makes  no  difference  and,  theoretically,  emission  occurs 
at  all  longitudes,  but  other  considerations  in  the  beaming  process  should 
provide  a  second  limit.  The  value  of  13°  for  5  when  Y=79°  was  suggested 
by  Goldreich  and  Lynden-Bell.  The  calculated  width  agrees  well  with 
sources  A  and  C,  but  overextends  source  B.  This  is,  of  course,  a  revisiting 
of  the  problem  of  comparing  sources  B]  and  B2  with  sources  A  and  C.  Angular 
separation  of  sources  A  and  C  in  sub-Io  longitude  is  greater  than  the 
separation  of  sources  B]  and  B2>  and  this  prevents  the  calculated  width 
of  one  of  the  source  pairs  from  matching  the  observed  extent  for  the  other. 
An  increase  of  6  from  13°  to  25°  is  necessary  to  achieve  the  same  source 
width  when  y  changes  from  79°  to  70°.  The  overextension  of  B  does  reveal 
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that  the  radiation  slants  into  the  source-D  region  on  the  low-longitude 
side.  Also,  the  Acuna-Ness  model  exhibits  the  upward  turn  often 
observed  on  the  right  side  of  source  C.  After  the  upward  turn,  the 
line  of  expected  radiation  follows  the  curve  of  the  "bridge"  observed 
between  sources  A  and  C  at  the  lower  frequencies.  This  will  be  discussed 
at  greater  length  in  the  next  section,  where  geometrical  limits  for  the 
beaming  directions  are  proposed. 

5.4  Radiation  Obscuration 

The  hypothesis  that  the  decametric  radiation  takes  place  near  the 
intersection  of  Io's  magnetic  flux  tube  with  the  Jovian  cloudtop  layer 
brings  up  the  question  of  when  these  regions  will  be  obscured  from  view 
by  the  planetary  limb.  Alexander  et  al .  (1975b)  answered  this  query  for 
the  simplest  configuration,  in  which  the  radiating  region  is  assumed  to 
be  at  the  Jovian  cloudtop  layer.  In  that  case,  the  active  areas  will  be 
visible  only  when  they  are  on  the  side  of  Jupiter  facing  the  earth.  In 
this  situation  the  longitudes  and  phases  of  Io  for  which  the  radiation 
is  obscured  are  the  shaded  sections  of  Figure  V-9,  which  has  been 
calculated  according  to  the  method  of  Alexander  et  al .  assuming  the 
declination  of  the  earth  to  be  zero.  No  warping  of  the  Io  flux  tube  or 
similar  symmetrizing  hypotheses  have  been  included.  Occultation  in  both 
the  northern  and  southern  hemispheres  is  shown.  Similar  calculations  for 
the  Smith  et  al .  (1975)  magnetic  field  model  showed  no  significant 
differences.  At  first,  the  figure  suggests  the  northern  hemisphere 
origins  of  sources  B  and  A,  while  C  would  propagate  from  the  south.  Emission 
from  cloudtop  level  is  not  realistic,  however.  If  radiation  occurs  at  the 
altitude  where  the  electron  gyrofrequency  equals  the  emitted  frequency, 
then  a  totally  different  picture  emerges. 
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Calculation  of  angles  of  obscuration  for  points  above  the  planet's 
cloud  surface  is  simplest  for  the  assumption  that  Dp  is  zero.  Figure  V-7b 
shows  an  example  in  which  all  radiation  in  the  ecliptic  plane  within 
the  angle  2a  is  cut  off.  Account  must  be  taken  of  the  rotational 
flattening  of  Jupiter,  and  this  has  been  done  through  the  equation 

R  =  R,  (1-e  cos26)  (V-8) 

where  R  is  the  planetary  radius  at  co-latitude  0,  R.  is  the  equatorial 
radius,  and  e  is  the  ellipticity  (see  Table  1-1).  The  correction  is 
important  since  the  appropriate  gyrofrequency  position  may  be  so  deep  in 
the  atmosphere  that  the  density  of  gases  makes  emission  impossible.  For 
ease  of  calculation,  the  cloudtop  layer  will  be  assumed  to  represent  this 
depth. 

Results  of  the  calculations  for  the  20-MHz  gyrofrequency  level  in 
both  the  northern  and  southern  hemispheres  are  shown  in  Figure  V-10.  The 
northern  hemisphere  shaded  regions  have  an  interesting  border  with  source 
B.  We  can  infer  that  the  high-longitude  edge  of  the  source  may  be  cut 
off  by  the  submersion  of  the  20-MHz  emission  region  in  the  cloud  layer. 
By  the  same  argument,  the  right  half  of  source  C  should  not  appear,  or  else 
we  must  again  assume  that  C  originates  in  the  southern  hemisphere.  Strict 
belief  in  the  calculated  locations  of  the  obscured  areas  is  not  advisable, 
since  the  equations  are  wery   sensitive  to  changes  in  altitude.  For  example, 
the  shaded  regions  for  the  22-MHz  level  in  Figure  V-ll  would  imply  a 
narrowing  of  source  B  by  nearly  15°,  and  the  same  computations  performed 
for  a  frequency  of  20  MHz  with  the  Smith  et  al.  (1975)  magnetic  field 
model  (Figure  V-12)  show  no  danger  of  any  source  obscuration.  We  conclude 
that  our  knowledge  of  the  magnetic  field  is  insufficient  for  accurate 
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calculations  near  the  planetary  surface.  The  Pioneer  spacecraft  did 
not  pass  close  enough  to  Jupiter  nor  remain  long  enough  for  an  accurate 
determination  of  the  magnetic  field  near  the  clouds  to  be  possible. 
The  diagrams  do  indicate,  however,  some  plausible  explanations  for  source 
morphology  in  terms  of  obscuration. 

The  more  likely  altitude  for  cutoff  of  the  decametric  radiation 
is  above  the  cloud  layer  and  dependent  on  the  plasma  density.  Refraction 
will  also  depend  on  this  factor,  and  so  the  amount  of  obscuration  may 
change  from  one  side  of  the  planet  to  the  other  as  plasma  parameters 
vary.  With  this  in  mind,  it  is  possible  to  cut  off  the  right  side  of 
source  B  by  submersion  of  the  radiating  region  without  endangering  a 
northern  hemisphere  source  C.  At  many  times  the  left  side  of  source  B 
also  exhibits  the  linear  slant  characteristic  of  a  submersion  cutoff.  If 
the  left  side  of  the  source  is  masked  at  20  MHz,  then  a  decrease  in 
gyrofrequency  to   15  MHz  would  increase  the  emission  altitude.  This 
would  allow  radiation  from  the  source  D  region  on  the  left  side  of  source  B 
as  shown  in  Figure  V-8.  At  the  same  time  removal  of  the  cutoff  would 
reveal  the  bridge  of  radiation  between  sources  A  and  C.  At  still  lower 
frequencies  no  obscuration  occurs  and  there  would  be  a  continuous  band 
of  radiation  at  Io  phases  of  90°  and  240°.  Preliminary  results  of  Desch 
and  Carr  (1977)  around  4  MHz  support  this  statement. 

Finally,  the  B]  and  B^  bulges  could  also  be  a  facet  of  source 
visibility.  An  inward-directed  radiation  cone  may  not  be  entirely 
reflected,  and  some  of  the  emission  may  escape  by  just  grazing  the  top 
of  the  reflection  layer.  In  this  sense  the  cone  angle  will  be  greater 
than  90°.  This  is  more  likely  to  occur  when  the  declination  of  the  earth 
is  positive,  and  that  is  the  time  the  B]  and  Bp  bulges  appear.  The 
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intervening  null  zone  would  still  be  an  effect  of  conditions  directly 
over  the  magnetic  pole,  such  as  the  lowered  amplification  factor 
proposed  by  Goertz.  Variation  of  emission  position  with  gyrofrequency 
would  be  a  possible  source   of  the  fluctuation  of  the  bulges  with  frequency, 

5.5  Summary 

In  the  course  of  Chapter  V  we  have  attempted  to  compare  the  second- 
order  source  characteristics  with  predictions  of  the  relevant  theories. 
Explanation  of  the  source-B  anomaly  has  led  to  speculations  concerning 
the  separate  existence  of  sources  B^  and  B2  and  their  possible  symmetry 
with  sources  A  and  C.  The  general  locations  of  the  Io-triggered  sources 
in  Io  phase  and  System  III  longitude  agree  well  with  predictions  assuming 
emission  in  a  conical  sheet  opening  at  an  angle  from  70°  to  90°  about  the 
magnetic  flux  line,  based  on  the  Pioneer  11  octopolar  magnetic  field 
models.  Beaming  directions  have  been  calculated  in  the  formalism  of  the 
Goldreich  and  Lynden-Bell  model  with  activity  occurring  at  longitudes  15° 
less  than  the  sub-Io  longitude  and  beamed  into  small  sections  of  the 
conical  sheet.  The  emission  geometry  proposed  by  many  other  perpendicular- 
beaming  theories  has  been  shown  to  be  nearly  interchangeable  with  the 
Goldreich  and  Lynden-Bell  system.  Finally,  many  of  the  second-order 
source  attributes  revealed  in  the  previous  chapters  may  be  explained  as 
occultation  of  the  emission  region  by  the  planet's  cloud  surface  or  by  an 
impenetrably  dense  plasma,  but  the  magnetic  field  is  not  yet  sufficiently 
determined  for  the  effects  to  be  computed  accurately.  Hope  for  the  near 
future  rests  in  the  Voyager  spacecraft,  one  of  which  is  intended  to 
penetrate  the  Io  flux  tube.  At  that  time  there  should  be  a  large  advance 
in  determining  the  origin  of  many  of  the  above  phenomena,  yet,  as  in 
all  past  space  ventures,  we  know  many  new  enigmas  will  arise.  In  any 
event,  Godspeed. 
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